Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,706,017 C→T 31.3% intergenic (‑227/+81) 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,706,0170CT31.3% 441.5 / 224.1 298intergenic (‑227/+81)dppB/dppAdipeptide/heme ABC transporter permease/dipeptide/heme ABC transporter periplasmic binding protein; dipeptide chemotaxis receptor
Reads supporting (aligned to +/- strand):  ref base C (101/86);  new base T (83/2);  total (210/88)
Fisher's exact test for biased strand distribution p-value = 2.43e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

Only 100 of 300 total aligned reads displayed.

GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGT  >  NC_000913/3705918‑3706114
                                                                                                    |                                                                                                          
  gcCACCGGTTATACAGATTGTTGACAAAGTGCCCTTTGTTTATGCCGGATGTGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATAc                                                                                                            <  2:3355360/98‑1 (MQ=255)
  gcCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAt                                                                                                          >  2:2606558/1‑100 (MQ=255)
   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                                                         >  2:2841549/1‑100 (MQ=255)
   ccACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACAtt                                                                                                         >  2:3938411/1‑100 (MQ=255)
       cGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑a                                                                                                   <  1:169121/99‑1 (MQ=255)
         ggTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑Gaa                                                                                               >  2:2171039/1‑100 (MQ=255)
            ataCAGATGGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAAtt‑‑‑t                                                                                         <  1:341159/100‑1 (MQ=255)
               cAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑g                                                                                    <  2:2995149/100‑1 (MQ=255)
                 gATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTa                                                                                  <  1:3150194/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAgg                                                                                <  1:247495/100‑1 (MQ=255)
                   ttgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAgg                                                                                >  2:340911/1‑100 (MQ=255)
                    tgttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGc                                                                               >  2:4266395/1‑100 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAG‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGcc                                                                              >  1:3088648/1‑100 (MQ=255)
                     gttgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGcc                                                                              <  2:715622/100‑1 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTGACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTg                                                                            <  1:2210155/100‑1 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTg                                                                            >  1:1683117/1‑100 (MQ=255)
                       tgACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTg                                                                            <  2:1470226/100‑1 (MQ=255)
                       aaaCAAAGCGCACTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTg                                                                            >  2:1834330/3‑100 (MQ=255)
                        gACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGa                                                                           >  1:3744466/1‑100 (MQ=255)
                            aagcgcacTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAg                                                                       >  2:2574175/8‑100 (MQ=255)
                             aGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGc                                                                      <  1:1976622/100‑1 (MQ=255)
                                  gcTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAgc                                                                 <  1:2371245/100‑1 (MQ=255)
                                    tttttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTACAAAAATCTTGCCAATTCAATATATTG‑C‑‑A‑GGACC‑‑‑AT‑‑‑GTAGGCCTGATAAGCGTAgcgc                                                               >  2:209628/10‑100 (MQ=18)
                                     ttGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCa                                                              <  2:1735579/100‑1 (MQ=255)
                                     ctttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTG‑C‑‑A‑GTTTA‑‑‑TTG‑‑GTAGGCCTGATAAgcgcagcgca                                                              >  2:2700622/9‑100 (MQ=21)
                                      tttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTG‑C‑‑A‑GTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCAt                                                             >  2:1117994/8‑100 (MQ=21)
                                      tGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCAt                                                             >  2:3772392/1‑100 (MQ=255)
                                      cttgtaggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑C‑‑A‑TGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAgcgc                                                               >  2:945490/8‑98 (MQ=17)
                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAATTC‑TGCAAATTCAACAAATTG‑C‑‑‑‑GAATCC‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                           >  2:3553759/7‑100 (MQ=18)
                                       ttgtaggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑C‑‑‑‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATc                                                            >  2:4343869/7‑100 (MQ=18)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                            >  1:141880/1‑100 (MQ=255)
                                       gTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                            <  2:3531107/100‑1 (MQ=255)
                                        tttATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCa                                                           >  1:2461817/1‑100 (MQ=255)
                                        tgtaggCCGGATAAGGCGTTCACGCCGCATCAGGCCTAC‑AAAGACCTGCAAATTCAATATATTG‑C‑‑A‑GGATA‑‑‑CGT‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                           >  2:2981499/6‑100 (MQ=17)
                                        tgtaggCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATc                                                            >  2:2264423/6‑99 (MQ=17)
                                         gtaggCCGGATAAGGCGTTCACGCTGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑C‑‑A‑CCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAg                                                          >  2:645012/5‑100 (MQ=18)
                                         gtaggCCGGATAAGGCGCTTACGCCGCATCCGGCCTAC‑AAGAACGTGCAAATTCAACAAATTGTG‑‑ATGCATC‑‑‑CT‑‑‑GTAGGCCTGATAAGCGTAGCGCATCa                                                           >  2:515598/5‑100 (MQ=18)
                                          taggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑C‑‑A‑GAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAgg                                                         >  1:1483394/4‑100 (MQ=17)
                                          taggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTG‑C‑‑‑‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAgg                                                         >  1:3253062/4‑100 (MQ=18)
                                          tATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAgg                                                         <  1:3801484/100‑1 (MQ=255)
                                           aggCCGGATAAGGCGTTTACGCCGCATCCGGCATAC‑AAAAGCGTGCAAATTCAATACATTG‑C‑‑A‑TGGGC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  1:4033404/3‑100 (MQ=17)
                                           aggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑C‑‑A‑TGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGc                                                        >  1:2863329/3‑100 (MQ=17)
                                           aTGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGc                                                        >  2:380383/1‑100 (MQ=255)
                                            tGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                       >  2:685377/1‑100 (MQ=255)
                                            ggctggataaggcgTTCACGCCGCATCCAGCCTAC‑AAAATTGTGCAAAATCAATGGATTG‑C‑‑A‑CAGCT‑‑‑AGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:2329378/11‑100 (MQ=25)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCACTATATTG‑C‑‑A‑GAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:732366/2‑100 (MQ=17)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑C‑‑A‑GAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  1:966751/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑C‑‑A‑ATTCT‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGct                                                       >  2:868667/2‑99 (MQ=18)
                                            ggCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAAATGTGCAAATTCAATAAATTG‑C‑‑‑‑‑AATTCAACTT‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGc                                                        >  2:1218092/2‑100 (MQ=21)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑T‑‑G‑GGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCa                                                       >  2:2075891/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAACCTTGCTAATTCAATATATTG‑C‑‑AGGGACT‑‑‑AT‑‑‑GTAGGCCTGATAAGCATAGCGCATCAGGCa                                                       >  2:2470596/2‑100 (MQ=18)
                                            ggCCGGATAAGGCGTTCACGCCGCATCCGGCATAC‑‑AAACCGCGCAAATTCAATATATTG‑CGGA‑GAAAA‑‑‑T‑‑‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                      >  2:3562998/2‑98 (MQ=21)
                                             gtcggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTG‑C‑‑A‑GGAAC‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                      >  2:3266383/10‑99 (MQ=21)
                                             gCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAt                                                      >  1:3624096/1‑100 (MQ=255)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑C‑‑‑‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                      >  1:4109146/1‑98 (MQ=18)
                                             gCCGGATAAGGCGTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAACATATTG‑C‑‑‑‑‑AATT‑‑‑CTCTTGTAGGCCTGATAAGCGTAGCGCATCAGGctg                                                      >  2:3930721/1‑98 (MQ=18)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑C‑‑A‑TGAGC‑‑‑CAC‑‑GTAGTCCTGATAAGCGTAGCGCATCAGGCag                                                      >  1:1008699/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAGATCGTGCAAATTCAATATATTG‑C‑‑A‑TGAGC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGACa                                                       >  2:1010024/1‑99 (MQ=17)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAATCGCTTAAATTCAATATATTG‑C‑‑A‑CCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaa                                                      >  2:168247/1‑99 (MQ=21)
                                             gCCGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTG‑C‑‑A‑GGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCag                                                      >  2:2227271/1‑99 (MQ=18)
                                             gCCGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCa                                                       >  1:3517417/1‑99 (MQ=17)
                                              ccggataaggcgTTTACGCCGCATCCGGCCTAC‑AAAAGCGTGCAAATTCAATACATTG‑C‑‑A‑TGGGC‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                     >  2:1252887/9‑98 (MQ=17)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAGTCGTATAAATTCAATGTATTG‑C‑‑A‑GTTTA‑‑‑TTG‑‑GTAGGCCTGATAAGCGCAGCGCATCAGGCaat                                                     <  1:1117991/92‑3 (MQ=25)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCATGAAAATTCAATACATTG‑C‑‑A‑AGATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCacg                                                     >  2:4129596/9‑98 (MQ=17)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCATGAAAATTCAATATATTG‑C‑‑A‑GGAGC‑‑‑TGC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCagt                                                     >  1:1203728/9‑98 (MQ=18)
                                              ccggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAAGCACGCAAACTCAATATATTG‑C‑‑A‑GAGAT‑‑‑CAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaat                                                     >  2:4052065/9‑98 (MQ=14)
                                              ccGGATAAGGCGTTCACGCCGCATCCGGCCTAC‑‑AAATCATGCAAATCCAAAATATTG‑C‑‑A‑GGAGT‑‑‑TAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAttt                                                    >  1:1279166/1‑100 (MQ=21)
                                              ccGGATAAGGCGTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAtt                                                     >  2:4174774/1‑100 (MQ=18)
                                               cggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAACTTTGCAAATTCAATATATTG‑C‑‑‑‑‑AATC‑‑‑TCCGTGTAGGCCTGATAAGCGTAGCGCATCAGGCaatt                                                    >  2:2769333/8‑97 (MQ=18)
                                               cggataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttt                                                    >  2:2418941/8‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑GAAATCGTGCTAATTCAATATATTG‑T‑‑G‑GGCTC‑‑‑CAC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgttt                                                   >  1:3320088/7‑100 (MQ=18)
                                                ggataaggcgTTCACGCCGCATCCGGCCTAC‑AAAATCACTTAAATTCAATATATTG‑C‑‑A‑CCATC‑‑‑CTC‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAattc                                                   >  1:2987483/7‑99 (MQ=21)
                                                ggataaggcgTTCACGCCGCATCCGACCCAC‑AAAGTCTTGCAAAATCAATAAATTG‑C‑‑A‑GGAAC‑‑‑GAT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCaacta                                                   >  1:1517992/7‑96 (MQ=21)
                                                 gataaggcgTTAACGCCGCATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCATGATAAGCGGAGCGCATCAGGCAttttt                                                  >  1:1375457/6‑100 (MQ=18)
                                                   taaggcgTTTACGCCGCATCCGGCCTAC‑AAAATCGTGCTAATTCAATATATTG‑C‑‑A‑GAAAC‑‑‑CTT‑‑GTAGGCCTGATAAGCGTAGCGCATCAGGCAgtttt                                                  >  1:1881978/4‑98 (MQ=18)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttt                                                  >  1:2104398/1‑98 (MQ=255)
                                                   tGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgc                                                >  2:2022169/1‑100 (MQ=255)
                                                    gcggcgTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTtttgcc                                               <  1:1834323/100‑1 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGCGCGCATCAGGCATTttttccttttc                                          >  1:1315856/1‑96 (MQ=255)
                                                           tAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCAttttttttttttttt                                        >  1:3867332/1‑90 (MQ=255)
                                                             aCGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATc                                      <  1:3283167/100‑1 (MQ=255)
                                                              cGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCa                                     >  1:2268796/1‑100 (MQ=255)
                                                                    aTCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATAAGGCATTTTTGCCTTTGCCATCAGTCtt                                >  1:571339/1‑99 (MQ=255)
                                                                      cctgCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             >  2:1970570/4‑100 (MQ=255)
                                                                      ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             <  1:846601/100‑1 (MQ=255)
                                                                      ccGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTa                             >  2:47104/1‑100 (MQ=255)
                                                                       cGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGTCTTTGCTTTTAGTCGTGTac                            >  1:4230648/1‑99 (MQ=18)
                                                                       cGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTAt                            >  1:2680948/1‑100 (MQ=255)
                                                                        ggCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATg                           >  1:531758/1‑100 (MQ=255)
                                                                            tAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                       <  1:2293729/100‑1 (MQ=255)
                                                                            tAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtt                       >  1:3912872/1‑100 (MQ=255)
                                                                              c‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCtttt                     >  1:1793546/1‑100 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGTCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                  <  2:442052/100‑1 (MQ=255)
                                                                                  aaTCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAat                  <  2:1883434/100‑1 (MQ=255)
                                                                                    tCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAatta                <  2:2710519/100‑1 (MQ=255)
                                                                                      gTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAattatt              <  2:3121656/100‑1 (MQ=255)
                                                                                          aaaTTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGAt          <  2:1347561/100‑1 (MQ=255)
                                                                                           aaTTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATa         <  1:3185528/100‑1 (MQ=255)
                                                                                               cAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAgaga     <  2:1372329/100‑1 (MQ=255)
                                                                                                  tACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTACTTATTCGATCGAGACGt  >  1:4062190/1‑100 (MQ=255)
                                                                                                    |                                                                                                          
GCGCCACCGGGTATACAGATTGTTGACAAAGTGCGCTTTGTTTATGCCGGATGCGGCGTTAACGCCTTATCCGGCCTAC‑AAAATCGTGCAAATTCAATACATTG‑C‑‑A‑GAATT‑‑‑TTC‑‑GTAGGCCTGATAAGCGGAGCGCATCAGGCATTTTTGCCTTTGCCATCAGTCTTGTATGGCTTTTAATTATTCGATAGAGACGT  >  NC_000913/3705918‑3706114

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