Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,705,984 T→G 39.6% intergenic (‑194/+114) dppB ← / ← dppA dipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,705,9840TG39.6% 416.8 / inf 332intergenic (‑194/+114)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base T (106/94);  new base G (126/5);  total (233/99)
Fisher's exact test for biased strand distribution p-value = 9.46e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-03

Only 100 of 338 total aligned reads displayed.

TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083
                                                                                                   |                                                                                                            
   ctgcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTAt                                                                                                           >  2:2558601/1‑100 (MQ=255)
     gcCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATcc                                                                                                         >  2:2894380/1‑100 (MQ=255)
        ggATGGCGCTGTACGCAGAGGCACAGCGCCACCGTGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                      <  1:1011624/100‑1 (MQ=25)
        ggATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGc                                                                                                      <  2:522833/100‑1 (MQ=37)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGcc                                                                                                     <  1:62179/100‑1 (MQ=255)
         gATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATAAGGCGTTAACGCCGCATCCGGca                                                                                                     <  2:1316653/100‑2 (MQ=37)
             gcgcTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑a                                                                                                >  1:4255298/1‑100 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCTGGATGCGGCGTTAACGCCTTATACGGCCTAC‑AAAAt                                                                                            >  1:1692412/1‑100 (MQ=255)
                 tGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTAGGCCGGATAAGGCGTTAACGCCGCATCCGGcataaacaaa                                                                                             <  2:4276910/100‑10 (MQ=21)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAAt                                                                                            <  2:4015599/99‑1 (MQ=255)
                  gCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATc                                                                                           >  1:437807/1‑100 (MQ=255)
                   cACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCg                                                                                          <  2:3665069/100‑1 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGCCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTCATCCGGCCTAC‑AAAATCGTg                                                                                        >  2:2690664/1‑100 (MQ=255)
                     cGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTg                                                                                        <  1:1842167/100‑1 (MQ=255)
                          agGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAAt                                                                                   >  1:2859126/1‑100 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                                <  2:3580955/100‑1 (MQ=255)
                             cacaGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCa                                                                                >  2:2124495/1‑100 (MQ=255)
                                  cgcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACa                                                                           <  2:1835161/100‑1 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                         >  1:1619533/1‑100 (MQ=255)
                                    ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                         >  2:922712/1‑100 (MQ=255)
                                     cACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTg                                                                        <  2:2727898/100‑1 (MQ=255)
                                     cACCGAGTATACAGATTGTTGACAAAGTGCGCTTTGTCTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                         <  1:1357764/99‑1 (MQ=255)
                                          ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga‑a                                                                 >  1:1493618/1‑100 (MQ=255)
                                          ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑Ga‑a                                                                 >  2:2473613/1‑100 (MQ=255)
                                                cAGATTTTTGACAAAGTGCGCTTTGTTTATGCCGGATGCTGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑At‑t‑‑‑tt                                                         <  1:3551303/98‑1 (MQ=255)
                                                  gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTa                                                   <  1:3733186/100‑1 (MQ=255)
                                                    gcataaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAgg                                                 >  2:975526/7‑100 (MQ=18)
                                                     tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGc                                                >  2:2514801/1‑100 (MQ=255)
                                                     tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGc                                                <  2:3929540/100‑1 (MQ=255)
                                                      gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGcc                                               <  1:867488/100‑1 (MQ=255)
                                                        tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTg                                             <  2:4063820/100‑1 (MQ=255)
                                                         gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCACATTCACTACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCt                                              >  1:1551668/1‑98 (MQ=255)
                                                          aCAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGAt                                           <  1:259877/100‑1 (MQ=255)
                                                              aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGaaaa                                         >  1:1657215/1‑98 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                     <  2:2403523/100‑1 (MQ=255)
                                                                tGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCgg                                     <  1:360453/100‑1 (MQ=255)
                                                                  cgcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAg                                   <  1:190388/100‑1 (MQ=255)
                                                                    cTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgcg                                 >  1:532550/1‑100 (MQ=255)
                                                                      ttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGC‑‑‑A‑AT‑C‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCa                               >  2:3031945/9‑100 (MQ=18)
                                                                      ctttgtgggtcGGATAAGGCGTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCA‑GG‑AA‑C‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCa                               >  1:2015075/11‑100 (MQ=21)
                                                                       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCAt                              >  2:449426/8‑100 (MQ=21)
                                                                       tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATc                             >  2:1251666/8‑100 (MQ=21)
                                                                       tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                              >  2:3691175/1‑100 (MQ=255)
                                                                       tGTTTATGCCGGAGGCTGCGTTAACGACTTACCCGGCATAC‑AAAATCGTGCCAATTCACTACATTGCA‑GA‑AT‑T‑‑‑TGC‑‑GTAGGCCTGATAAGCGGAgcg                                 <  1:3117850/97‑1 (MQ=17)
                                                                       cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGTTCGTGCAAATTCAATAAATTGCA‑‑A‑TA‑T‑‑‑GAC‑‑GTAGGCCTGATAAGCGTAGCGCATc                             >  2:4157011/8‑100 (MQ=18)
                                                                        gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                             >  2:1471074/1‑100 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            <  1:754853/100‑1 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            >  2:1127174/1‑100 (MQ=255)
                                                                         tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                            <  2:3407641/100‑1 (MQ=255)
                                                                          gtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAACATATTGCA‑‑A‑TTCT‑‑‑CTT‑‑GTAGGCATGATAAGAGTAACGCATCAg                           >  1:2275138/5‑100 (MQ=12)
                                                                          gtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AG‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCa                            >  2:3738793/5‑99 (MQ=17)
                                                                           taggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  1:4184301/4‑100 (MQ=17)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATTGTGTAAATTCAATAAATTGAT‑CG‑AT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  2:153646/4‑100 (MQ=21)
                                                                           taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTGCA‑GG‑AG‑C‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                          >  1:2083808/4‑100 (MQ=17)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                         >  1:1750264/3‑100 (MQ=14)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1335584/3‑100 (MQ=21)
                                                                            aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTGCA‑GG‑AG‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:1998899/3‑100 (MQ=18)
                                                                             tGCCGGATGCGGCGTTAACTACATATCCGGCCTAC‑AAAATCGAGCAAATTCAATACATTGCA‑GA‑AC‑T‑‑‑TTC‑‑GTACGCCTGATAAGCGAAGCGCAACAGGc                         >  1:3194132/1‑99 (MQ=17)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  2:3325062/1‑100 (MQ=255)
                                                                             tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                        >  1:1295317/1‑100 (MQ=255)
                                                                             ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTGCA‑‑A‑TTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                        >  1:2915855/2‑99 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:3323310/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  1:1321910/2‑100 (MQ=18)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTGCA‑GT‑TT‑A‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCa                        >  1:657558/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                        >  2:753226/2‑100 (MQ=21)
                                                                             ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTGCA‑AT‑CT‑C‑‑‑CGT‑‑GTAGTCCTGATAAGCGTAGCGCATCAgg                          >  1:4079169/2‑98 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:2798223/1‑99 (MQ=17)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  1:3923365/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                       >  2:665223/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaa                       >  2:4168318/1‑99 (MQ=18)
                                                                              gCCGGATAAGGCGTTCACTCCGCATCCGGCCTAC‑AAAAGCACGCAATTTCAATATATTGCA‑GA‑GA‑T‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:3056236/1‑99 (MQ=14)
                                                                              gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTGCA‑CC‑AT‑C‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                       >  1:1085725/1‑99 (MQ=21)
                                                                               cgggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:2978717/9‑98 (MQ=17)
                                                                               ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  2:4344735/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTGCA‑TG‑AT‑C‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAACCCATCAAGCag                       >  2:2871569/9‑98 (MQ=17)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGCAAATTCAATAGATTGCA‑GA‑GA‑T‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  1:2341855/9‑98 (MQ=14)
                                                                               ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAACATGCAAATTCAATAGATTGCA‑GAGAT‑T‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                      >  2:1942664/9‑98 (MQ=17)
                                                                               ccggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                      >  2:1555057/9‑100 (MQ=18)
                                                                               ccGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                      >  1:2759637/1‑98 (MQ=18)
                                                                                cggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCATGCAAATTCAATATATTGCA‑GA‑GT‑T‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtt                     >  1:222894/8‑100 (MQ=17)
                                                                                cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTGTGCAAATTCAATATATTGCA‑GG‑AA‑A‑‑‑CAC‑‑GTAGGCCTGATAAGCGAAGCGCATCAGGCagtt                     >  2:3195072/8‑97 (MQ=18)
                                                                                aggataaggcgTTCACGCCGAATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GA‑GA‑T‑‑‑GGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                     >  1:2682450/8‑100 (MQ=14)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTGCA‑GA‑AA‑C‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:2073193/7‑100 (MQ=18)
                                                                                 ggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTGC‑‑‑A‑AT‑TCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaatt                     >  2:2369838/7‑97 (MQ=21)
                                                                                 ggataaggcgTTCACGCCGCATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GA‑GA‑T‑‑‑GGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAtttt                    >  1:3475694/7‑100 (MQ=14)
                                                                                 ggataaggcgTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTGCA‑GG‑AC‑C‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                    >  2:3068104/7‑100 (MQ=18)
                                                                                 ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTGCAGGA‑AC‑G‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                    >  2:2282018/7‑96 (MQ=25)
                                                                                  gataaggcgTTTACGCCGCATCCGGCCTAC‑AAAACAGTGCAAATTCAATATATTGCA‑AT‑CT‑T‑‑‑CTG‑‑GTAGGCCTGATAAGCATCGCGCATCAGGCaagctt                   >  2:1402029/6‑95 (MQ=21)
                                                                                  gataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgctc                   >  1:4111492/6‑95 (MQ=17)
                                                                                    taaggcgTTCACGCCGCATCCTGCCTAC‑AAAAGCACGCAAATTCAATATATTGCA‑GA‑GA‑T‑‑‑GGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCATTtttgc                 >  2:4205153/4‑100 (MQ=14)
                                                                                      cggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcct               >  1:3437667/1‑100 (MQ=255)
                                                                                       ggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTGCA‑AG‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacgctcgcatg              >  1:326761/1‑90 (MQ=18)
                                                                                       ggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctt              >  2:224091/1‑100 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgc           >  1:2547126/1‑100 (MQ=255)
                                                                                          gTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttg            >  2:1078257/1‑99 (MQ=255)
                                                                                           ttAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcctttgcc          <  1:1043851/100‑1 (MQ=255)
                                                                                            tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCa         >  2:1740001/1‑100 (MQ=255)
                                                                                              aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc       <  2:1296512/100‑1 (MQ=255)
                                                                                                  cTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTc   <  2:100552/100‑1 (MQ=255)
                                                                                                   ttATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCt  <  2:2639548/100‑1 (MQ=255)
                                                                                                   |                                                                                                            
TTGCTGCCGGATGGCGCTGCACGCAGAGGCACAGCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTGCA‑GA‑AT‑T‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCT  >  NC_000913/3705885‑3706083

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑