Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_001416 47,317 N→T ?87I (NTC→ATC)  lambdap78 ← putative envelope protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00141647,3170NT100.0% 103.5 / NA 34?87I (NTC→ATC) lambdap78putative envelope protein
Reads supporting (aligned to +/- strand):  ref base N (0/0);  new base T (17/17);  total (17/17)

ACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGANGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTG  >  NC_001416/47268‑47363
                                                 |                                               
aCCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGAt                                                 <  1:3877/50‑1 (MQ=255)
 ccTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATg                                                <  1:3331/50‑1 (MQ=255)
  cTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGt                                              >  1:1013/1‑50 (MQ=255)
  cTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGt                                              >  1:2856/1‑50 (MQ=255)
  cTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCAAACACGATGt                                              <  1:252/50‑1 (MQ=255)
    ttttAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGa                                            <  1:2664/50‑1 (MQ=255)
     tttAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGaa                                           >  1:157/1‑50 (MQ=255)
       tAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGCTGAat                                          >  1:4029/1‑50 (MQ=255)
       tAAGTCCTTGTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAata                                         >  1:918/1‑50 (MQ=255)
        aaGTCATTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAatat                                        <  1:3161/50‑1 (MQ=255)
         aGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAAtatt                                       <  1:3080/50‑1 (MQ=255)
           tCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGATGTGAAtattat                                     <  1:5000/50‑1 (MQ=255)
                   ttaACACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCt                             <  1:3790/50‑1 (MQ=255)
                     aaCACGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAg                           <  1:322/50‑1 (MQ=255)
                      acacGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGa                          >  1:3171/1‑50 (MQ=255)
                       cacGGTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGAt                         <  1:4434/50‑1 (MQ=255)
                          ggTGTTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGt                      >  1:4390/1‑50 (MQ=255)
                             gtTATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTaaa                   <  1:4027/50‑1 (MQ=255)
                              ttATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAat                  >  1:2762/1‑50 (MQ=255)
                              ttATCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAat                  >  1:4512/1‑50 (MQ=255)
                               tATCGTTTTTTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAata                 <  1:207/50‑1 (MQ=255)
                                aTCGTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAatat                <  1:3093/50‑1 (MQ=255)
                                 tCGTTTTCCAACACGATGTGAATATTATCTGTGGCTAGATAGTAAAtata               <  1:2994/50‑1 (MQ=255)
                                  cGTTTTCTAACACGATGTGAATATTATTTGTGGCTAGATAGTAAATATaa              <  1:2244/50‑1 (MQ=255)
                                   gTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATaag             >  1:1903/1‑49 (MQ=255)
                                   gTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAAt             >  1:95/1‑50 (MQ=255)
                                   gTTTTCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAAt             <  1:887/50‑1 (MQ=255)
                                    ttttCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAAtg            >  1:451/1‑50 (MQ=255)
                                    ttttCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAAtg            >  1:2264/1‑50 (MQ=255)
                                    ttttCTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAAtg            >  1:738/1‑50 (MQ=255)
                                        cTAACACGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTgag        >  1:545/1‑50 (MQ=255)
                                            cacGATGTGACTATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGt    <  1:3141/50‑1 (MQ=255)
                                             acGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGtt   >  1:4583/1‑50 (MQ=255)
                                             acGATGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGtt   >  1:2972/1‑50 (MQ=255)
                                              cGATGTGAATATTATCTGTGGCTAGATAGTAAATATTATGTGAGACGTtg  <  1:623/50‑1 (MQ=255)
                                                 |                                               
ACCTTTTTAAGTCCTTTTATTAACACGGTGTTATCGTTTTCTAACACGANGTGAATATTATCTGTGGCTAGATAGTAAATATAATGTGAGACGTTG  >  NC_001416/47268‑47363

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: