Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913974,4490AG24.0% 89.6 / 10.3 51D44G (GAT→GGT) mukFchromosome partitioning protein MukF
Reads supporting (aligned to +/- strand):  ref base A (13/25);  new base G (12/0);  total (26/25)
Fisher's exact test for biased strand distribution p-value = 8.57e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GCTATGAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGCGTCTTTCTATG  >  NC_000913/974207‑974699
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                 
gCTATGAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTTAACGGCGAGCGGCTGGGTgggtgg                                                                                                                                                                                                                                                            >  1:385112/1‑246 (MQ=255)
gCTATGAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTga                                                                                                                                                                                                                                                             <  1:319969/249‑1 (MQ=255)
           ttACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCCAGCGTCTGGATGGGGAGATGGGTGAAggg                                                                                                                                                                                                                                                >  1:149051/1‑250 (MQ=255)
             ccTTGAAATAAAAACGGTTTATTGCCGCCAGGATCCGTATATTACCCAGTGGCGTTATACTCATGTCACCGCGCGCAATCCGCAGAGCAGGGATAAAGTATGAGTGAATTTTCCCAGCCAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTGCTCCGTCGCGCTGCCGGTAGGCCGACTCTCTTTTCTGCTGGCGGTTGCGACGCTGAACGGCGAGCGTCGGGATGGTCAGATGAGTGAAGGCg                                                                                                                                                                                                                                               <  2:16666/249‑1 (MQ=255)
                       gtaaCGCGTGATTGCCGTCATGAACCGTATAGTACCCTGGTGCGTTATATTAATGTCACCGTGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCACTGTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGAt                                                                                                                                                                                                                                    <  2:500973/249‑1 (MQ=255)
                       gAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGGCTGGGTGGTGGGATGAGTGAAGGCGGGCTGGGGGAt                                                                                                                                                                                                                                    >  1:222572/1‑251 (MQ=255)
                       aaaaCGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGAGGGTGAGATGAGTGAAGGCGAGCTGGTGGAt                                                                                                                                                                                                                                    <  1:198935/250‑1 (MQ=255)
                                tATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTTCCACGCTGAACGGCGAGCGTCTGGATGGGGGGGTGAGTGAAGGCGGGCTGGTGGAGGCATTCCGc                                                                                                                                                                                                                           >  1:709608/1‑251 (MQ=255)
                                   tGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCa                                                                                                                                                                                                                        <  1:285463/250‑1 (MQ=255)
                                       gTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGGGAAGGCGGGCTGGGGGATGCCTTCCCGCAtgtgg                                                                                                                                                                                                                   >  1:304417/1‑250 (MQ=255)
                                           ccttccctacgctctgcggatcggTTATAGTAATGTCACCGCGCGCAAACCGAAGAACAAGGATAAAGTATGAGTGAATTTTCCTATACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACGTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTTGTGGATGCATTCCGCCATgtgagtga                                                                                                                                                                                                               <  2:284330/228‑1 (MQ=255)
                                                       tACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCACCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTTCCCCGCTGTACGGCGAGCGTCTGGTCGGTGAGATGGGTGAAGGCGAGCTTGTGGGTGACTTTCGCCATGTGATGGATGCGTTTGAGc                                                                                                                                                                                                   >  1:178177/1‑250 (MQ=255)
                                                             ggggCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAAcc                                                                                                                                                                                              <  1:38567/249‑1 (MQ=255)
                                                                 cGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGCCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGGGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGaa                                                                                                                                                                                        <  2:23160/251‑1 (MQ=255)
                                                                                                      aGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACa                                                                                                                                                  <  1:28859/251‑1 (MQ=255)
                                                                                                        gATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCGGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATg                                                                                                                                                <  1:705475/251‑1 (MQ=255)
                                                                                                           aaaGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGACAAAATGACTTCTCCCTCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGGGAAGGCGCGCTGGGGGATGCATTCCGCCATGTGGGTGATGTGTTTGAGCCAACCCGCGAAACCATCGGCGTGGGCGCCCATAACACGGTAAACGACCTGGTg                                                                                                                                             >  1:414680/1‑251 (MQ=255)
                                                                                                           aaaGTATGAGTGAATTGTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGt                                                                                                                                              <  1:667060/250‑1 (MQ=255)
                                                                                                             aGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTTGACGGCGAGCGTCTGGGTGGGGAGATGAGGGAAGGCGAGCTGGGGGGTGCATTCCGCCCTGTGAGTGATGCGTTTGAGCGAACCAGCGAAAACATCGGGGTGGGCGCCCAAAAAGCGATCAACGACATGGGGCg                                                                                                                                           >  1:526374/1‑251 (MQ=255)
                                                                                                                   aGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCGCCATCTCGCTGCCGGTAGACCGACTCTCTTGTCTGCTGGCGGTTGCCACGCTGTACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTg                                                                                                                                                                                                       <  2:461781/187‑1 (MQ=255)
                                                                                                                      gAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTct                                                                                                                                  <  1:96213/251‑1 (MQ=255)
                                                                                                                             cccAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTTCCCCGCTGAACGGGGAGAGTCTGGATGGTGGGATGAGGGACAGCGAGCTAGGGGAAGGATTCCCCCCATGTAGTGTAGCGTTTGAGGCAAACCGCGCAACCCTACGCGTGCGCGCCCAATACGCGCACAACGACATAGTGCGTCaaactactactaac                                                                                                                           >  2:421545/1‑239 (MQ=255)
                                                                                                                                aGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCt                                                                                                                                                                                                                                           <  2:111567/139‑1 (MQ=255)
                                                                                                                                aGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGGCTGGGTGGTGAGATGAGTGAAGGCGAGCt                                                                                                                                                                                                                                           >  1:111567/1‑139 (MQ=255)
                                                                                                                                   cAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGGGGGATGAGTGAAGGCGGGCTGGGGGGTGCATTCCGCCCTGTGAGTGATGCGTTTGAGCAAACCAGCGGAACCATCGGCGTGGGCGCCCATAACCGGATCAACGACAAGGTGCGTCAACggcctccggaccccttt                                                                                                                     >  1:161911/1‑235 (MQ=255)
                                                                                                                                       ccccGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCTACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCa                                                                                                               <  1:480177/251‑1 (MQ=255)
                                                                                                                                        cccGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGTGAAGGCGGGCTGGGGGGTGCATTCCGCCATGTGAGTGATGCGGTTGAGCAAACCAGCGAAACCATCGGCGGGGGCGCCAAAAACGCGATCAACGACATGGTGGGTCAACGGCTGCTGAAACGCCTTAcc                                                                                                                 >  1:459694/1‑249 (MQ=255)
                                                                                                                                                ggTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAg                                                                                                       <  2:222572/250‑1 (MQ=255)
                                                                                                                                                 gTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGGAAGGCGAGCGTCTGGATGGGGAGATGGGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGGGATGCGTTTGAGGAAACCAGCGGAAACATCGGCGTGGGCGCCAATAAAGCGGTCACCGACACGGGGCGTCCAAGTCTGCTGCACCCCTTTACCCGCGAGCAggg                                                                                                     >  1:106261/1‑250 (MQ=255)
                                                                                                                                                      cTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGTGAAGGCGAGCTGGGGGATGCATTCCGCCCTGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCCATAACGCGATCAAAGACATGGTGCGGCAACGTCTGCTGAACCCCCTTACCAGCGAGCAGGCggg                                                                                                  >  1:269770/1‑248 (MQ=255)
                                                                                                                                                                  aaaTGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCCCCCTGAACCGCGGGCGTCTGGGTGGTGAGATGGGGGAAGGCGAGCTGGGGGGTGCATTCCGCCATGTGAGTGATGCGGTTGAGCAAACCCGCGAAACCATCGGCGTGCGGGCCAAAAAAGCCGTCAAACACATGGTGCGTCCACGTCCTCTGCACCCCCTTTACAGCGAGCAGGGGGGAGGGGACGCCAttt                                                                                   >  1:326676/1‑251 (MQ=255)
                                                                                                                                                                  aaaTGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAAttt                                                                                   <  1:282886/251‑1 (MQ=255)
                                                                                                                                                                     tGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTTGGTGGGGAGGTTAGGGAAGGCGAGCTGGTGGGTGCATTCCGCCCTGTGAGTGATGCGTTTGAGCAAAACAGCGAAACCATCGGCGTGGGCCCCCATAAAGCGATCCAAGACAAGGTGCGGCAACGTCCGCTGGAACCCTTTTACAGCGGGCAGGCGGAAggggaagcaaatttcc                                                                                >  1:208782/1‑238 (MQ=255)
                                                                                                                                                                          tctcCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTc                                                                             <  1:312332/249‑1 (MQ=255)
                                                                                                                                                                             ccATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGAcg                                                                        <  1:275609/251‑1 (MQ=255)
                                                                                                                                                                               aTCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCCCGCTGAACGGCGAGCGTCTGGGTGGTGGGATGAGGGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGGTTGAGCAAACCAGCGAAACCATCGGCGGGCGCGCCCAAAAAGCGGACAACGACATGGGGCGGCCACGGCTGCTTAACCGCTTTTCCAGCGAGCAGGGGGAAGGGGACGCAATTTATCGCCTGAcgcc                                                                      >  1:171110/1‑251 (MQ=255)
                                                                                                                                                                                      tGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCCCGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGGGAAGGCGAGCTGGGGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGGGCGCCCATAACGCGGACAACGAAATGGGGGGGCCACGTCTGCTGAAACGCCTTACCCAGGGGCAGGGGGGAGGGGAAGGCATTTACCGTCTGACGCCGCccggc                                                               >  1:258306/1‑251 (MQ=255)
                                                                                                                                                                                           gTAGACCGACTCTCTTTGCGGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCtcggcatcgg                                                          <  2:16726/251‑1 (MQ=255)
                                                                                                                                                                                                      ctctTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGa                                                <  1:521677/249‑1 (MQ=255)
                                                                                                                                                                                                        ctTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCCCTTTACCAGCGAGCAGGCGGGAGGGAACGCAATTTACCGTCTGACGCCGCtcggcagatcgggagagaccccc                                             >  1:162967/1‑234 (MQ=255)
                                                                                                                                                                                                        ctTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCtcggc                                                               <  2:162967/233‑1 (MQ=255)
                                                                                                                                                                                                               gctgGCGGTTGCGATGCTGAACGGCGAGCTTCTTGATGGTGAGATTAGTGAAGGTGAGATGGTGGAGGCATTCCGCCATTGGAGTGATTCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCTAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCTGAAGGGAACGCAATTTACCGGCTGACGCCGCCCGGCATCGGCATTACTGACTACTACag                                       <  2:12577/249‑2 (MQ=255)
                                                                                                                                                                                                               gctgGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATc                                      <  1:297265/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCCCGGCATCGGCATTACTGACTACTACATCCGTc                                  >  1:208629/1‑251 (MQ=255)
                                                                                                                                                                                                                    cGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGCCATGGTGCGTCAACGTCTGCTGGACCGCTTTACCAGCGCGCAGGCGGAAGGGAACGCAATTTACCGTCTGAAGCCGCTCGGCCTCGGCATTACCGACTACTTCATCCGttac                                >  2:162371/1‑248 (MQ=255)
                                                                                                                                                                                                                       ttGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCACCGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCCATTTACCGTCTGCCGCCGACCGGCCTCGGCCTTACTGACTACCACATCCGTCAgcgc                             >  2:679789/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCTGAACGGCGAGCGTCTGGATGGTGAGATGCGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCACCGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGGAGGGAACGCAATTTCCCGACTGCCGCCGCCACGCATCACACTTACTGACTACTAACTCCGTCAgacagagtttt                      >  2:705877/1‑241 (MQ=255)
                                                                                                                                                                                                                                   aCGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGGGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCTGCATTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACg                 <  2:258306/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGc             >  1:639751/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACCTCCGTCAGCGCGAGTTTTCTACGCTGCGTCt         >  1:516442/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ggATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTCCATCCGTCAGCGCGGGTTTTCTACGCTGCGTCTTTCt     >  1:243735/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  aTGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGCGTCTTTCTATg  <  1:596518/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                 
GCTATGAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCCATGTGAGTGATGCGTTTGAGCAAACCAGCGAAACCATCGGCGTGCGCGCCAATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTGAACCGCTTTACCAGCGAGCAGGCGGAAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCATTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGCGTCTTTCTATG  >  NC_000913/974207‑974699

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