Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913458,1890AG20.6% 142.4 / 10.5 63D255G (GAT→GGT) clpXATP‑dependent Clp protease ATP‑binding subunit ClpX
Reads supporting (aligned to +/- strand):  ref base A (16/34);  new base G (13/0);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 6.48e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGC  >  NC_000913/457950‑458431
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
aGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑GTTGCCGGGCTGGAT‑AAA‑GTGAttt                                                                                                                                                                                                                                           >  1:121650/1‑251 (MQ=255)
 gCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTc                                                                                                                                                                                                                                          <  2:237851/251‑1 (MQ=255)
      gCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTTCCGGTCTGGGT‑AAA‑GTGATTTcccccc                                                                                                                                                                                                                                     >  1:63994/1‑248 (MQ=255)
         tGGTATTGTCTACATCGATGAACTCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGGCGGTCTGGGT‑AAA‑GTGATTTTCCAACgtg                                                                                                                                                                                                                                  >  2:273696/1‑251 (MQ=255)
             aTTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGAAAGCCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCGGGCGCTGTTGAAACTGAGCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGTCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATTCCTCTAATATCCTGGTTATTTGTGGCGGTGCG‑TTTGCCGGTTTGGAT‑AAA‑GTGATTTCCCACCGTGTAg                                                                                                                                                                                                                               <  1:332172/250‑1 (MQ=255)
               tGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGaa                                                                                                                                                                                                                             <  1:371349/250‑1 (MQ=255)
                      aTCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTTCCGGTCTGGAT‑AAA‑GTGGTTTCCCCCCGGGTAGAAAACGGCTc                                                                                                                                                                                                                     >  1:134925/1‑251 (MQ=255)
                          aTGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGACACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTATAAGATCCTGTTTATTTGTGGCGGTGGG‑TTTGCCGGGCTGGATCAAA‑TGGATTTCC‑ACCGTGTAGGAACAGGCTCCGGc                                                                                                                                                                                                                 >  2:354758/1‑251 (MQ=255)
                                ccGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACGGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGt                                                                                                                                                                                                           <  1:153830/250‑1 (MQ=255)
                                              cGTAAGGCAGACAACCCGTCCATTACCCGAGACGTTGCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCGAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGGGGAGAAACCGGCTCCGTCATTGGTTTTGGCTCGAc                                                                                                                                                                                                <  2:326486/248‑1 (MQ=255)
                                                 aaGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGGGCGACGG‑Taa                                                                                                                                                                                          >  1:9415/1‑250 (MQ=255)
                                                   gTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTTCCGGGCTGGAT‑AAA‑GTGATTTTCCACCGCGGAGAAACCGGCTCCGGCACTGGGTTTTGCGCGACGG‑Taaa                                                                                                                                                                                         >  1:136395/1‑249 (MQ=255)
                                                             ccGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGa                                                                                                                                                                             <  1:280783/251‑1 (MQ=255)
                                                               gTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCGCTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGAGCG‑GTTGCCGGGCTGGGA‑AAA‑GTGATTTTCCACCGTGTAGAAACCGGCTCCGGCATTTGTTTTGGCGCGGCGG‑Taaaaaagaaagacgcc                                                                                                                                                                            >  1:268972/1‑238 (MQ=255)
                                                               gTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATGCCTCTAAGATCCTGGTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACa                                                                                                                                                                           <  2:360688/251‑1 (MQ=255)
                                                                         cGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTTCCGGGCTGGGT‑AAA‑GGGATTTTCCACCGGGTAGAAACCGGCTCCGGCATTGGTTTTGGGGCGGCGG‑TAAAAGCGAAGACCGGCAaaacaaacg                                                                                                                                                                  >  2:132989/1‑247 (MQ=255)
                                                                           agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAACCTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TGTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAaagcaa                                                                                                                                                                     <  1:90400/245‑1 (MQ=255)
                                                                           agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTTCCGGTCTGGGT‑AAA‑GTGATTTTCCACCGTGTAGAAAACGGCGCCCGCCTTGGGTTTGGCGCGGCGG‑GAAAAGCGAAGGCCGGCaaacaaaggcg                                                                                                                                                                 >  2:21783/1‑241 (MQ=255)
                                                                           agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCAGGATGGGT‑AAAAGGGGTTTCCC‑CCGTGGAAAAAACGGATCCGGCGTTGGTTTTGGCGGGGAGG‑TAAAAACGAAGTACGGCaaacgaaaaaagg                                                                                                                                                               >  2:90400/1‑241 (MQ=255)
                                                                             ccGTTGCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGc                                                                                                                                                            <  2:56596/250‑1 (MQ=255)
                                                                             aCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑GTTGTCGGGCTGGGT‑AAA‑GTGGTTTCCCACCGGGTAGAAACCGGCTCCGGCGTTGGTTTTGGGGCGGCGGGTAAAA‑AGGAGTCCCGCAAAGCAAGC‑GAAggg                                                                                                                                                            >  1:279263/1‑250 (MQ=255)
                                                                              cGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCg                                                                                                                                                           <  1:280467/251‑1 (MQ=255)
                                                                                tttCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑GTTGCCGGTCTGGAT‑AAA‑GGGATTTCCCACCGGGTAGAAACCGGACCCGGCATTGGGTTTTGCGCGACGG‑GCAAAGCGAAGTGCGGCAaagtaagcagaggcg                                                                                                                                                            >  1:175627/1‑242 (MQ=255)
                                                                                 ttCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAg                                                                                                                                                         <  2:252750/250‑1 (MQ=255)
                                                                                    cGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAgctgc                                                                                                                                                     <  1:9380/251‑1 (MQ=255)
                                                                                     ggTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAgctgct                                                                                                                                                    <  1:227526/251‑1 (MQ=255)
                                                                                        gAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGgc                                                                                                                                                 <  1:160994/251‑1 (MQ=255)
                                                                                                 cagcagGCACTGTTGAAACTGATCGAAGGTACGGTCGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGGCAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGa                                                                                                                                        <  2:125030/251‑1 (MQ=255)
                                                                                                  agcagGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGaa                                                                                                                                       <  1:425762/251‑1 (MQ=255)
                                                                                                     agGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTTGC‑GGGCATGAT‑AAA‑ATGGTTTTCCACCGGGTAGAAAACGGCTCCGGCATTGGTTTTGGCGCGACGG‑TGAAAGCGAAGTCCGACCAACCAAGC‑GAAGGCGAGCTGCCGGCGCAGGGTGAAcc                                                                                                                                     >  2:273737/1‑250 (MQ=255)
                                                                                                     agGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCg                                                                                                                                    <  2:139748/251‑1 (MQ=255)
                                                                                                         aCTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAg                                                                                                                                <  1:177926/251‑1 (MQ=255)
                                                                                                             ttgcaaCTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTTGC‑GGTCTGGGT‑AAC‑ATGATTTCCCCCCGGGGAGAAAACGGCTCCGCCATTGGGTTTGTTGCGACCG‑GAAAAACGAAGCACGAAACAGCAACG‑GAACGCGCGCTGCTGGCGCCGGttacacaggaagagct                                                                                                                            >  2:152904/5‑237 (MQ=255)
                                                                                                               gAAAATGAACGAAGGTAAGGTCGCGGCTGTTCAACCGAAAGGTGGGCGTACATATCCGCAGCGGGAATTCTTGCAGGTTGCTAGCTCTGCGAGCCGGTTTATTTGTGGCGGTGTG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAAGCGGCTCCGGCATTGGTCGTGGCGCGACGG‑TAAAAGCGAAGTGCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGACCCGGAAGATCtta                                                                                                                          <  2:180002/251‑3 (MQ=255)
                                                                                                                     gATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCaa                                                                                                                      <  2:207120/249‑1 (MQ=255)
                                                                                                                        cGAAGGTACGTTAGCTGCTGTTCCACAGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGAGCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTg                                                                                                                 <  2:373046/251‑1 (MQ=255)
                                                                                                                             gTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCtt                                                                                                            <  1:132989/251‑1 (MQ=255)
                                                                                                                               aCGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGTCGGTCTGGGT‑AAA‑GTGATTTCCCCCCGTGGTGAAAACGGCACCCGCATTGGTTTTGGCGCGGCGG‑TTAAAGCGAAATCACACCAAGAAAGCCGAAG‑CGGGCTGCTGGCGCAGGGTGAACCGGACGAACTGCTTAAATTTGGTCTTAt                                                                                                          >  2:386758/1‑251 (MQ=255)
                                                                                                                                     gctgctGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑GAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAAACGGAAGATCTGATCAAGGTTGGTCTTATCCCTg                                                                                                     >  1:284725/1‑250 (MQ=255)
                                                                                                                                        gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGt                                                                                                  <  2:329522/250‑1 (MQ=255)
                                                                                                                                        gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGt                                                                                                  <  2:279263/250‑1 (MQ=255)
                                                                                                                                                aCCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGGT‑AAA‑GGGATTTCCCACCGGGGAGAAAACGGCTCCGGCATTGGGTTTGGCGGGACGG‑GAAAAGCGAAGGCCGGCAAAGCAAGC‑GAAGGGGAGCTGCTGGCGCCGGGTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGgtc                                                                                         >  1:23070/1‑251 (MQ=255)
                                                                                                                                                          gggctgaattctcagcagcaGGAATGCTTGCAGGTTGATACCTCTAAGATCCTGTTTATGTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCGGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTGATTGGTCGCCTGCCgg                                                                                <  2:366167/236‑1 (MQ=255)
                                                                                                                                                          gggCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCCCCGTGTAGAAACCGGCTCCGGCATTGGGTTTTGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCCGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCggg                                                                               >  1:27347/1‑250 (MQ=255)
                                                                                                                                                           ggCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGG‑GTTGCCGGTCTGGGT‑AAA‑GTGATTTTCCACCGTGTAGAAAACGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCGGCTGGCGCAGGTTGAAACGGAAGATGTGGTCAAGGTTGGTCTTATCCCCGAGGTTATTGGTCGTGTGCCGGtt                                                                              >  1:60987/1‑251 (MQ=255)
                                                                                                                                                                   cATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGGT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTTGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGGCTTATCCCTGAGTTTTTTGGTCGGCTGCCGGGTGTCGCAAc                                                                      >  1:49279/1‑251 (MQ=255)
                                                                                                                                                                   cATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑GTTGCCCGTCTGGAT‑AAA‑GTGATTTTCCACCGTGTAGAAACCGGCCCCGGCATTGGTTTTGGCGCGCCGG‑GAAAAGCGGAGTCCGACAAAAGAAGC‑GAAGGCGAGCTGCTGGCGCCGGGTGAAACGGAAGATCTGATCAAGTTTGGTCTTTTCCCTGAGGTTTTTGGGCGGCTGCCGGGTGTCGCAAc                                                                      >  1:184897/1‑251 (MQ=255)
                                                                                                                                                                      ccGCAGCAGGAATTCTTGCAGGTTGATACCGCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGtt                                                                   <  2:87392/251‑1 (MQ=255)
                                                                                                                                                                       cGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTAGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTtg                                                                  <  1:259463/251‑1 (MQ=255)
                                                                                                                                                                        gcagcaGGAATTCTTGCAGGTTGATTCCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTtga                                                                 <  1:105212/251‑1 (MQ=255)
                                                                                                                                                                           gcagGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCTGGGCTGGAA‑AAA‑GTGATTTCCCACCGTGTAGAAAACCGCACCCGCATTGGTTTTGGCGCGACGG‑GAAAAACGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGAAGGGTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTAATGGTCGTCTGCCGGTTGTCGCAACGTtgaagg                                                              >  2:239118/1‑249 (MQ=255)
                                                                                                                                                                           gcagGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCgg                                                                                                                                                                                                                  <  2:265534/105‑1 (MQ=255)
                                                                                                                                                                                 aTTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTATGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCCGCCCCGGCATTGGTTTTGGCGCGGCGG‑TAAAAGCGAAGTCCGACAAAACAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGGACCGGAAGATCTGGTCAAGTTTGGTCTTATCCCTGAGGTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTg                                                         >  1:398919/1‑250 (MQ=255)
                                                                                                                                                                                        cAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGGTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCgaaga                                                 <  2:127672/251‑1 (MQ=255)
                                                                                                                                                                                          ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGGCCGGGT‑GAA‑GGGATTTTCCACCCGGTAGAAACCGGCACCGGCATTGGTTTTGGCGCGGCGG‑GAAAAACGAAGTCCCGCAAAGCAAAC‑GAAGGCGAGCGGCTGGGGGAGGGTGAACCGGGAAAACTGGACAAGTTTGGTCTTATCCCTGAGGTTTTTGGTCGGCTGCCGGGTGTGGCAACGATGGATGAACTGAGCgacgaa                                                >  1:6876/1‑250 (MQ=255)
                                                                                                                                                                                           gTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑GTTGCCCGTCTGGAT‑AAA‑AGGATTTCCCACCGGGTAGAAACCGGCCCCGGCATTGGTTTTGGCGCGACGG‑TGAAAGCGAAGGCCGACAAAGCAAGC‑GAAGGGGAGCCGCTGGCGCCGGTTGAACCGGCAGAGCTGATCAAGTTTGGTCTTCTCCCCGGGGTTATTTGTCGTCTGCCGGTTGTCGCCACGTTGGATGAACTGGGCgcagcagc                                              >  2:70175/1‑244 (MQ=255)
                                                                                                                                                                                                                  taattggggCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGCTTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCTTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAagag                        <  2:268972/243‑1 (MQ=255)
                                                                                                                                                                                                                   aTTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGc                       <  2:252874/250‑1 (MQ=255)
                                                                                                                                                                                                                      tgtgGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAaagcaagc                                                                                                                                                                   <  1:121531/108‑1 (MQ=255)
                                                                                                                                                                                                                       gtgGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAATAGCCGaaa                  <  1:406081/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTAAGATCCTCAAAGAGCCGAAAAACGc             >  2:165618/1‑251 (MQ=255)
                                                                                                                                                                                                                             gTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGcc            >  2:286866/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGc  >  2:169723/1‑249 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
AGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCG‑TTTGCCGGTCTGGAT‑AAA‑GTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGG‑TAAAAGCGAAGTCCGACAAAGCAAGC‑GAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGC  >  NC_000913/457950‑458431

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