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.

GCTATGAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCA‑TTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGCGTCTTTCTATG  >  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)
                                   tGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑a                                                                                                                                                                                                                        <  1:285463/250‑1 (MQ=255)
                                       gTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGGGAAGGCGGGCTGGGGGATGCCTTCCCGC‑Atgtgg                                                                                                                                                                                                                   >  1:304417/1‑250 (MQ=255)
                                           ccttccctacgctctgcggatcggTTATAGTAATGTCACCGCGCGCAAACCGAAGAACAAGGATAAAGTATGAGTGAATTTTCCTATACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACGTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTTGTGGATGCATTCCGCC‑ATgtgagtga                                                                                                                                                                                                               <  2:284330/228‑1 (MQ=255)
                                                       tACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCACCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTTCCCCGCTGTACGGCGAGCGTCTGGTCGGTGAGATGGGTGAAGGCGAGCTTGTGGGTGACTTTCGCC‑ATGTGATGGATGCGTTTGAG‑c                                                                                                                                                                                                   >  1:178177/1‑250 (MQ=255)
                                                             ggggCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAAcc                                                                                                                                                                                              <  1:38567/249‑1 (MQ=255)
                                                                 cGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGCCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGGGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGaa                                                                                                                                                                                        <  2:23160/251‑1 (MQ=255)
                                                                                                      aGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACa                                                                                                                                                  <  1:28859/251‑1 (MQ=255)
                                                                                                        gATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCGGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATg                                                                                                                                                <  1:705475/251‑1 (MQ=255)
                                                                                                           aaaGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGACAAAATGACTTCTCCCTCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGGGAAGGCGCGCTGGGGGATGCATTCCGCC‑ATGTGGGTGATGTGTTTGAG‑CCAACCCGCGAAACCATCGGCGTGGGCGCC‑CATAACACGGTAAACGACCTGGTg                                                                                                                                             >  1:414680/1‑251 (MQ=255)
                                                                                                           aaaGTATGAGTGAATTGTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGt                                                                                                                                              <  1:667060/250‑1 (MQ=255)
                                                                                                             aGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTTGACGGCGAGCGTCTGGGTGGGGAGATGAGGGAAGGCGAGCTGGGGGGTGCATTCCGCC‑CTGTGAGTGATGCGTTTGAG‑CGAACCAGCGAAAACATCGGGGTGGGCGCC‑CAAAAAGCGATCAACGACATGGGGCg                                                                                                                                           >  1:526374/1‑251 (MQ=255)
                                                                                                                   aGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCGCCATCTCGCTGCCGGTAGACCGACTCTCTTGTCTGCTGGCGGTTGCCACGCTGTACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTg                                                                                                                                                                                                       <  2:461781/187‑1 (MQ=255)
                                                                                                                      gAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTct                                                                                                                                  <  1:96213/251‑1 (MQ=255)
                                                                                                                             cccAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTTCCCCGCTGAACGGGGAGAGTCTGGATGGTGGGATGAGGGACAGCGAGCTAGGGGAAGGATTCCCCCCATGT‑AGTGTAGCGTTTGAGGCAAACC‑GCGCAACCCTACGCGTGCGCGCCCAATA‑CGCGCACAACGACATAGTGCGTCaaactactactaac                                                                                                                           >  2:421545/1‑239 (MQ=255)
                                                                                                                                aGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCt                                                                                                                                                                                                                                           <  2:111567/139‑1 (MQ=255)
                                                                                                                                aGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGGCTGGGTGGTGAGATGAGTGAAGGCGAGCt                                                                                                                                                                                                                                           >  1:111567/1‑139 (MQ=255)
                                                                                                                                   cAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGGGGGATGAGTGAAGGCGGGCTGGGGGGTGCATTCCGCC‑CTGTGAGTGATGCGTTTGAG‑CAAACCAGCGGAACCATCGGCGTGGGCGCC‑CATAACCGGATCAACGACAAGGTGCGTCAACggcctccggaccccttt                                                                                                                     >  1:161911/1‑235 (MQ=255)
                                                                                                                                       ccccGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCTACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑a                                                                                                               <  1:480177/251‑1 (MQ=255)
                                                                                                                                        cccGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGTGAAGGCGGGCTGGGGGGTGCATTCCGCC‑ATGTGAGTGATGCGGTTGAG‑CAAACCAGCGAAACCATCGGCGGGGGCGCC‑AAAAACGCGATCAACGACATGGTGGGTCAACGGCTGCTG‑AAACGCCTTAcc                                                                                                                 >  1:459694/1‑249 (MQ=255)
                                                                                                                                                ggTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAg                                                                                                       <  2:222572/250‑1 (MQ=255)
                                                                                                                                                 gTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGGAAGGCGAGCGTCTGGATGGGGAGATGGGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGGGATGCGTTTGAG‑GAAACCAGCGGAAACATCGGCGTGGGCGCC‑AATAAAGCGGTCACCGACACGGGGCGTCCAAGTCTGCTG‑CACCCCTTTACC‑CGCGAGCAggg                                                                                                     >  1:106261/1‑250 (MQ=255)
                                                                                                                                                      cTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGTGAAGGCGAGCTGGGGGATGCATTCCGCC‑CTGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑CATAACGCGATCAAAGACATGGTGCGGCAACGTCTGCTG‑AACCCCCTTACC‑AGCGAGCAGGCggg                                                                                                  >  1:269770/1‑248 (MQ=255)
                                                                                                                                                                  aaaTGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCCCCCTGAACCGCGGGCGTCTGGGTGGTGAGATGGGGGAAGGCGAGCTGGGGGGTGCATTCCGCC‑ATGTGAGTGATGCGGTTGAG‑CAAACCCGCGAAACCATCGGCGTGCGGGCC‑AAAAAAGCCGTCAAACACATGGTGCGTCCACGTCCTCTGCACCCCCTTTAC‑‑AGCGAGCAGGGGG‑GAGGGGACGCCAttt                                                                                   >  1:326676/1‑251 (MQ=255)
                                                                                                                                                                  aaaTGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAAttt                                                                                   <  1:282886/251‑1 (MQ=255)
                                                                                                                                                                     tGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTTGGTGGGGAGGTTAGGGAAGGCGAGCTGGTGGGTGCATTCCGCC‑CTGTGAGTGATGCGTTTGAG‑CAAAACAGCGAAACCATCGGCGTGGGCCCC‑CATAAAGCGATCCAAGACAAGGTGCGGCAACGTCCGCTGGAACC‑CTTTTAC‑AGCGGGCAGGCGG‑AAggggaagcaaatttcc                                                                                >  1:208782/1‑238 (MQ=255)
                                                                                                                                                                          tctcCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTc                                                                             <  1:312332/249‑1 (MQ=255)
                                                                                                                                                                             ccATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGAcg                                                                        <  1:275609/251‑1 (MQ=255)
                                                                                                                                                                               aTCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCCCGCTGAACGGCGAGCGTCTGGGTGGTGGGATGAGGGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGGTTGAG‑CAAACCAGCGAAACCATCGGCGGGCGCGCC‑CAAAAAGCGGACAACGACATGGGGCGGCCACGGCTGCTT‑AACCGCTTTTCC‑AGCGAGCAGGGGG‑AAGGGGACGCAATTTATCGCCTGAcgcc                                                                      >  1:171110/1‑251 (MQ=255)
                                                                                                                                                                                      tGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCCCGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGGGAAGGCGAGCTGGGGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGGGCGCC‑CATAACGCGGACAACGAAATGGGGGGGCCACGTCTGCTG‑AAACGCCTTACCCAG‑GGGCAGGGGGGAGGGGAAGGCA‑TTTACCGTCTGACGCCGCccggc                                                               >  1:258306/1‑251 (MQ=255)
                                                                                                                                                                                           gTAGACCGACTCTCTTTGCGGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCtcggcatcgg                                                          <  2:16726/251‑1 (MQ=255)
                                                                                                                                                                                                      ctctTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCA‑TTACTGa                                                <  1:521677/249‑1 (MQ=255)
                                                                                                                                                                                                        ctTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGGTGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCCCTTTACC‑AGCGAGCAGGCGG‑GAGGGAACGCAATTTACCGTCTGACGCCGCtcggcagatcgggagagaccccc                                             >  1:162967/1‑234 (MQ=255)
                                                                                                                                                                                                        ctTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCtcggc                                                               <  2:162967/233‑1 (MQ=255)
                                                                                                                                                                                                               gctgGCGGTTGCGATGCTGAACGGCGAGCTTCTTGATGGTGAGATTAGTGAAGGTGAGATGGTGGAGGCATTCCGCC‑ATTGGAGTGATTCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCT‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCTG‑AAGGGAACGCAATTTACCGGCTGACGCCGCCCGGCATCGGCA‑TTACTGACTACTACag                                       <  2:12577/249‑2 (MQ=255)
                                                                                                                                                                                                               gctgGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCA‑TTACTGACTACTACATc                                      <  1:297265/250‑1 (MQ=255)
                                                                                                                                                                                                                  ggCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCCCGGCATCGGCA‑TTACTGACTACTACATCCGTc                                  >  1:208629/1‑251 (MQ=255)
                                                                                                                                                                                                                    cGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGCCATGGTGCGTCAACGTCTGCTG‑GACCGCTTTACC‑AGCGCGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGAAGCCGCTCGGCCTCGGCA‑TTACCGACTACTTCATCCGttac                                >  2:162371/1‑248 (MQ=255)
                                                                                                                                                                                                                       ttGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCACCGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCCATTTACCGTCTGCCGCCGACCGGCCTCGGCC‑TTACTGACTACCACATCCGTCAgcgc                             >  2:679789/1‑251 (MQ=255)
                                                                                                                                                                                                                              gCTGAACGGCGAGCGTCTGGATGGTGAGATGCGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCACCGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑GAGGGAACGCAATTTCCCGACTGCCGCCGCCACGCATCA‑CACTTACTGACTACTAACTCCGTCAgacagagtttt                      >  2:705877/1‑241 (MQ=255)
                                                                                                                                                                                                                                   aCGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGGGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCTGCA‑TTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACg                 <  2:258306/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCA‑TTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGc             >  1:639751/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCA‑TTACTGACTACTACCTCCGTCAGCGCGAGTTTTCTACGCTGCGTCt         >  1:516442/1‑251 (MQ=255)
                                                                                                                                                                                                                                                ggATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCA‑TTACTGACTACTCCATCCGTCAGCGCGGGTTTTCTACGCTGCGTCTTTCt     >  1:243735/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  aTGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCA‑TTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGCGTCTTTCTATg  <  1:596518/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                 
GCTATGAAGCATTACTTGAATTAGAAACGCGTTATTGCCGTCAGGAACCGTATATTACCCTGGGGCGTTATATTCATGTCACCGCGCGCAAACCGCAGAGCAAGGATAAAGTATGAGTGAATTTTCCCAGACAGTCCCCGAACTGGTTGCCTGGGCCAGAAAAAATGACTTCTCCATCTCGCTGCCGGTAGACCGACTCTCTTTTCTGCTGGCGGTTGCCACGCTGAACGGCGAGCGTCTGGATGGTGAGATGAGTGAAGGCGAGCTGGTGGATGCATTCCGCC‑ATGTGAGTGATGCGTTTGAG‑CAAACCAGCGAAACCATCGGCGTGCGCGCC‑AATAACGCGATCAACGACATGGTGCGTCAACGTCTGCTG‑AACCGCTTTACC‑AGCGAGCAGGCGG‑AAGGGAACGCAATTTACCGTCTGACGCCGCTCGGCATCGGCA‑TTACTGACTACTACATCCGTCAGCGCGAGTTTTCTACGCTGCGTCTTTCTATG  >  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: ‑