Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,208,1320GA67.9% 34.0 / 45.8 54intergenic (‑600/‑15)purH/rrsEbifunctional AICAR transformylase/IMP cyclohydrolase/16S ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base G (9/8);  new base A (24/12);  total (34/20)
Fisher's exact test for biased strand distribution p-value = 3.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGC  >  NC_000913/4207893‑4208367
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           
aaaaCAGCCGGTTAGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAGGATGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCTGACAATCTGTGTGGGCACTAGAAGATAAGGATTCTTAACGTCGCCAGACGAAAACTGAATAACAAGTCTCAAGAGTGAACACGAAATTCATTACGAAATTTAATTCTTACAGCATAAAACt                                                                                                                                                                                                                                  >  2:167867/2‑251 (MQ=11)
  aaCAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACttt                                                                                                                                                                                                                                <  2:42282/251‑1 (MQ=38)
             cggttGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTTAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAag                                                                                                                                                                                                                      <  2:17261/250‑1 (MQ=255)
                        aaaaaTCCTGAAATTCAGGGGTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAAGCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATg                                                                                                                                                                                                          <  1:306725/251‑1 (MQ=255)
                             tCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTc                                                                                                                                                                                                      >  1:251334/1‑250 (MQ=25)
                                         gggTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTg                                                                                                                                                                                         <  1:345405/251‑1 (MQ=35)
                                                     aaaGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCt                                                                                                                                                                              >  1:54803/1‑250 (MQ=25)
                                                         agGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAacac                                                                                                                                                                         >  1:12380/1‑251 (MQ=25)
                                                            aaaGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATg                                                                                                                                                                      >  2:386867/1‑251 (MQ=25)
                                                               gCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACGGCTCGTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTGTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCaa                                                                                                                                                                   <  1:44696/251‑1 (MQ=255)
                                                                cGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAg                                                                                                                                                                  <  1:342442/251‑1 (MQ=25)
                                                                 gTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACGTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGt                                                                                                                                                                 <  2:97476/251‑1 (MQ=25)
                                                                      ataCGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAAc                                                                                                                                                            <  1:264199/251‑1 (MQ=25)
                                                                       taCGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACg                                                                                                                                                           >  1:456054/1‑251 (MQ=25)
                                                                        aCGCCACCTCGCGTCAGTGCGCTAAAGCGCTTCGCAACTGCTCTTTAACAATTTATCAGACAATCTATGTGGGCACTCGAGGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAAAACTTGATAGTCGAACgt                                                                                                                                                          >  2:356203/1‑250 (MQ=11)
                                                                        aCGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACgg                                                                                                                                                          >  1:285422/1‑251 (MQ=25)
                                                                               cTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAgg                                                                                                                                                   >  1:110591/1‑251 (MQ=25)
                                                                                               aaaGCGCGTCGCAACAGCTCTTAAACACTTTCTCAGACAATCTGTGTGGGCCCTCGAAGATACGGATTCTTAACGTCGACAGACGACAAATGAATACCAAGTCTCCTGAGTGAACCCGTAATTCATTACGAAGTTTAATTCTTTACGCATCAAACTTTTCACTTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCACGTCGCACGGTAACAGGAAGAAGCTTGcttctt                                                                                                                                   >  2:201288/1‑251 (MQ=11)
                                                                                                             cTGCTCTCTAACAATCTATCAGACCATCTGTGAGGGCCCTCGAAGCTACGGATTCTTAACCACGCAAGACGCAAAATGAATACCAAGCCTCAAGAGTGAACACGTAATTCATTACGACGTTTAATTCTTTGAGCATCAAACTTTTAAATTGCAGAGTTTGATCATGGCTCAGCTTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGATGACGAGTGGc                                                                                                                     >  2:279404/1‑251 (MQ=11)
                                                                                                               gCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCg                                                                                                                    <  2:111030/250‑1 (MQ=14)
                                                                                                                  cggtAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCGTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTTTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTGGCTGACGAGTGGCAGAc                                                                                                                 <  2:94853/247‑1 (MQ=9)
                                                                                                                     tAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACggg                                                                                                              <  2:324120/250‑1 (MQ=31)
                                                                                                                     tAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACggg                                                                                                              <  1:386867/250‑1 (MQ=14)
                                                                                                                        cAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTg                                                                                                            >  1:128804/1‑249 (MQ=14)
                                                                                                                          aTTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGt                                                                                                         >  1:97801/1‑250 (MQ=31)
                                                                                                                                 aGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTg                                                                                                 <  1:158166/251‑1 (MQ=31)
                                                                                                                                     aaTCTGTGTGGGCACTCGAAGATACGGATTCTTCACGTCGCAAGACGAAAAATGAATACCACGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGACGAGTTTGATCATGGCGCAGATTGAACGCTGTCGGCAGGACTAACACATGCAAGCCGCACGGTAACAGGAAGAAGCTTGCTTCATTGCTGCCGAGGGACGGACGGGGGGACAATGTCTGGGa                                                                                              >  2:78200/1‑250 (MQ=11)
                                                                                                                                          gtgtgGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTg                                                                                         <  2:430884/250‑1 (MQ=14)
                                                                                                                                          gtgtgGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGc                                                                                        >  1:192996/1‑251 (MQ=14)
                                                                                                                                          gtgtgGGCACTCGAAGATACGGATACTTACCGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGCACACGTAATTCATTACGAAGTTTAATTATTTGAGCATCACACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGAGTAATGTCTGGGAAACTGc                                                                                        >  2:232977/1‑251 (MQ=255)
                                                                                                                                                  actcgaatatctggATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGAGAAAACGTAATTAATTACGAAGATTAATTCTTTGAGCATCAAACTTTTAACTTGAAGAGTTTGATCATGGTTCTGATTGAACGATGGCGGCAGGCCTCACACACGACAGTCGAACGGTAACAGGAAGAAGCTTGCTGCTTAGAATAAGAGTGGCGGACGGGAGAGTCATGCATGGGaaagtggctggttg                                                                                 >  2:64013/13‑239 (MQ=2)
                                                                                                                                                  aCTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGa                                                                                <  1:312901/251‑1 (MQ=31)
                                                                                                                                                       aaGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAggggg                                                                           >  1:422216/1‑251 (MQ=31)
                                                                                                                                                               gATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAAAACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAactac                                                                   >  1:312544/1‑251 (MQ=14)
                                                                                                                                                                 ttCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGCAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGAAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGTGAAACTGCCTGATGGAGGGGGATAACTACTg                                                                 >  1:278232/1‑251 (MQ=9)
                                                                                                                                                                                    aaaaaTGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTAAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGGGAGTAATGGCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCg                                              >  2:40278/1‑251 (MQ=14)
                                                                                                                                                                                      aaaTGAATACCAAGTCTCCAAAGTGAACACGTAATTCATTACGACGTTAAATTCTTTAACCATCCAACTTTTAAATTGAAGAGTTTGATCATAGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGAAGGGTGAGTAATGTCTGGGAAACTGCCTGATGGGGGGGGATAACTACTGGAAACGGTAGATCATACCGc                                             >  2:179115/1‑250 (MQ=11)
                                                                                                                                                                                        aTGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGTGGGATAACTACTGGAAACGGTAGCTAATACCGCATa                                          <  2:98499/251‑1 (MQ=14)
                                                                                                                                                                                          gAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATAAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAAc                                        <  1:332625/251‑1 (MQ=14)
                                                                                                                                                                                               acaagtaccacgagTGAACACGTAATTCATTACGACGTTTACTTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTAATGGAAACGGTAGCTAATACAGCATAACGTCGc                                   >  2:65610/12‑251 (MQ=11)
                                                                                                                                                                                                      tcAAGCGTGAACACGTAATTCATTACGACGTTTAATTCTTTGAGCATCAAACTTTTACATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGCAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTCCTGGAACCGGTAGCTAATACCGCATAACGGCGCAAGAcc                             >  2:457691/1‑250 (MQ=11)
                                                                                                                                                                                                      tcAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAggggggt                                                                         >  2:140752/1‑204 (MQ=14)
                                                                                                                                                                                                      tcAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGAt                                                                         <  1:140752/206‑1 (MQ=14)
                                                                                                                                                                                                       cAAGAGTGCACACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGCAACTGCCTGATGGAGGGGGATAACTCCTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCa                            >  2:395595/1‑250 (MQ=9)
                                                                                                                                                                                                        aaGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCaa                           >  1:359148/1‑250 (MQ=31)
                                                                                                                                                                                                             tGAACTCGTAATTCATTACGACGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTAAACGCTGGCGGCAGGCCTAACACCTGCAAGTCGAACGGTCACAGGCAGCAGCTTGCTTCTTTGCTGACGAGTGGCGGCAGGGTGAGTAATGTCTGGGAAACTGACTGCTGGAGGGGGATAGCTACGGGAAACGGTAGCTCATACCGCATAACGTAGCacaaccaatggtgg                     >  2:63199/1‑238 (MQ=11)
                                                                                                                                                                                                              gAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAgggg                    >  1:289554/1‑251 (MQ=31)
                                                                                                                                                                                                               aaCACGTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGCGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAggggg                   >  1:300895/1‑251 (MQ=14)
                                                                                                                                                                                                                    gTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCtt              <  2:262639/251‑1 (MQ=31)
                                                                                                                                                                                                                    gTAATTCATTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCtt              >  2:367769/1‑251 (MQ=14)
                                                                                                                                                                                                                           aTTACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGCAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTACAAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGCCGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATACAGTCGCAAGACCAAAGTGGGGGACCTTCGGGcac       >  2:320771/1‑249 (MQ=11)
                                                                                                                                                                                                                             tACGAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGCTCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGCACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGAGCGTAATGTCTGGGCAACTACCTGATGGAGGGGGATATCTACTGGAAACGGTAGCTAATACCGCATAACGTCCCCAGACCAAAGCGGGGGACCTTCGGGCctct     >  2:105925/1‑251 (MQ=11)
                                                                                                                                                                                                                              aCGAAGTTTAATTCTTTGAGCATCAAACTTGTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCtt    <  2:78886/251‑1 (MQ=14)
                                                                                                                                                                                                                                gAAGTTTAATTCTTTAAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGc  >  1:107388/1‑251 (MQ=11)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                           
CAAACAGCCGGTTCGGTTGAAGAGAAAAATCCTGAAATTCAGGGTTGACTCTGAAAGAGGAAAGCGTAATATACGCCACCTCGCGACAGTGCGCTAAAGCGCGTCGCAACTGCTCTTTAACAATTTATCAGACAATCTGTGTGGGCACTCGAAGATACGGATTCTTAACGTCGCAAGACGAAAAATGAATACCAAGTCTCAAGAGTGAACACGTAATTCATTACGAAGTTTAATTCTTTGAGCATCAAACTTTTAAATTGAAGAGTTTGATCATGGCTCAGATTGAACGCTGGCGGCAGGCCTAACACATGCAAGTCGAACGGTAACAGGAAGAAGCTTGCTTCTTTGCTGACGAGTGGCGGACGGGTGAGTAATGTCTGGGAAACTGCCTGATGGAGGGGGATAACTACTGGAAACGGTAGCTAATACCGCATAACGTCGCAAGACCAAAGAGGGGGACCTTCGGGCCTCTTGC  >  NC_000913/4207893‑4208367

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: