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.

AGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGC  >  NC_000913/457950‑458431
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
aGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGCCGGGCTGGATAAAGTGAttt                                                                                                                                                                                                                                           >  1:121650/1‑251 (MQ=255)
 gCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTc                                                                                                                                                                                                                                          <  2:237851/251‑1 (MQ=255)
      gCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCCGGTCTGGGTAAAGTGATTTcccccc                                                                                                                                                                                                                                     >  1:63994/1‑248 (MQ=255)
         tGGTATTGTCTACATCGATGAACTCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGGCGGTCTGGGTAAAGTGATTTTCCAACgtg                                                                                                                                                                                                                                  >  2:273696/1‑251 (MQ=255)
             aTTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGAAAGCCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCGGGCGCTGTTGAAACTGAGCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGTCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATTCCTCTAATATCCTGGTTATTTGTGGCGGTGCGTTTGCCGGTTTGGATAAAGTGATTTCCCACCGTGTAg                                                                                                                                                                                                                               <  1:332172/250‑1 (MQ=255)
               tGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGaa                                                                                                                                                                                                                             <  1:371349/250‑1 (MQ=255)
                      aTCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCCGGTCTGGATAAAGTGGTTTCCCCCCGGGTAGAAAACGGCTc                                                                                                                                                                                                                     >  1:134925/1‑251 (MQ=255)
                          aTGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGACACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTATAAGATCCTGTTTATTTGTGGCGGTGGGTTTGCCGGGCTGGATCAAATGGATTTCCACCGTGTAGGAACAGGCTCCGGc                                                                                                                                                                                                                 >  2:354758/1‑251 (MQ=255)
                                ccGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACGGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGt                                                                                                                                                                                                           <  1:153830/250‑1 (MQ=255)
                                              cGTAAGGCAGACAACCCGTCCATTACCCGAGACGTTGCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCGAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGGGGAGAAACCGGCTCCGTCATTGGTTTTGGCTCGAc                                                                                                                                                                                                <  2:326486/248‑1 (MQ=255)
                                                 aaGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGGGCGACGGTaa                                                                                                                                                                                          >  1:9415/1‑250 (MQ=255)
                                                   gTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCCGGGCTGGATAAAGTGATTTTCCACCGCGGAGAAACCGGCTCCGGCACTGGGTTTTGCGCGACGGTaaa                                                                                                                                                                                         >  1:136395/1‑249 (MQ=255)
                                                             ccGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGa                                                                                                                                                                             <  1:280783/251‑1 (MQ=255)
                                                               gTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCGCTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGAGCGGTTGCCGGGCTGGGAAAAGTGATTTTCCACCGTGTAGAAACCGGCTCCGGCATTTGTTTTGGCGCGGCGGTaaaaaagaaagacgcc                                                                                                                                                                            >  1:268972/1‑238 (MQ=255)
                                                               gTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATGCCTCTAAGATCCTGGTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACa                                                                                                                                                                           <  2:360688/251‑1 (MQ=255)
                                                                         cGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCCGGGCTGGGTAAAGGGATTTTCCACCGGGTAGAAACCGGCTCCGGCATTGGTTTTGGGGCGGCGGTAAAAGCGAAGACCGGCAaaacaaacg                                                                                                                                                                  >  2:132989/1‑247 (MQ=255)
                                                                           agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAACCTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTGTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAaagcaa                                                                                                                                                                     <  1:90400/245‑1 (MQ=255)
                                                                           agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCCGGTCTGGGTAAAGTGATTTTCCACCGTGTAGAAAACGGCGCCCGCCTTGGGTTTGGCGCGGCGGGAAAAGCGAAGGCCGGCaaacaaaggcg                                                                                                                                                                 >  2:21783/1‑241 (MQ=255)
                                                                           agaCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCAGGATGGGTAAAAGGGGTTTCCCCCGTGGAAAAAACGGATCCGGCGTTGGTTTTGGCGGGGAGGTAAAAACGAAGTACGGCaaacgaaaaaagg                                                                                                                                                               >  2:90400/1‑241 (MQ=255)
                                                                             ccGTTGCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGc                                                                                                                                                            <  2:56596/250‑1 (MQ=255)
                                                                             aCGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGTCGGGCTGGGTAAAGTGGTTTCCCACCGGGTAGAAACCGGCTCCGGCGTTGGTTTTGGGGCGGCGGGTAAAAAGGAGTCCCGCAAAGCAAGCGAAggg                                                                                                                                                            >  1:279263/1‑250 (MQ=255)
                                                                              cGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCg                                                                                                                                                           <  1:280467/251‑1 (MQ=255)
                                                                                tttCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGCCGGTCTGGATAAAGGGATTTCCCACCGGGTAGAAACCGGACCCGGCATTGGGTTTTGCGCGACGGGCAAAGCGAAGTGCGGCAaagtaagcagaggcg                                                                                                                                                            >  1:175627/1‑242 (MQ=255)
                                                                                 ttCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAg                                                                                                                                                         <  2:252750/250‑1 (MQ=255)
                                                                                    cGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctgc                                                                                                                                                     <  1:9380/251‑1 (MQ=255)
                                                                                     ggTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctgct                                                                                                                                                    <  1:227526/251‑1 (MQ=255)
                                                                                        gAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGgc                                                                                                                                                 <  1:160994/251‑1 (MQ=255)
                                                                                                 cagcagGCACTGTTGAAACTGATCGAAGGTACGGTCGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGGCAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGa                                                                                                                                        <  2:125030/251‑1 (MQ=255)
                                                                                                  agcagGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGaa                                                                                                                                       <  1:425762/251‑1 (MQ=255)
                                                                                                     agGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTTGCGGGCATGATAAAATGGTTTTCCACCGGGTAGAAAACGGCTCCGGCATTGGTTTTGGCGCGACGGTGAAAGCGAAGTCCGACCAACCAAGCGAAGGCGAGCTGCCGGCGCAGGGTGAAcc                                                                                                                                     >  2:273737/1‑250 (MQ=255)
                                                                                                     agGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCg                                                                                                                                    <  2:139748/251‑1 (MQ=255)
                                                                                                         aCTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAg                                                                                                                                <  1:177926/251‑1 (MQ=255)
                                                                                                             ttgcaaCTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTTGCGGTCTGGGTAACATGATTTCCCCCCGGGGAGAAAACGGCTCCGCCATTGGGTTTGTTGCGACCGGAAAAACGAAGCACGAAACAGCAACGGAACGCGCGCTGCTGGCGCCGGttacacaggaagagct                                                                                                                            >  2:152904/5‑237 (MQ=255)
                                                                                                               gAAAATGAACGAAGGTAAGGTCGCGGCTGTTCAACCGAAAGGTGGGCGTACATATCCGCAGCGGGAATTCTTGCAGGTTGCTAGCTCTGCGAGCCGGTTTATTTGTGGCGGTGTGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAAGCGGCTCCGGCATTGGTCGTGGCGCGACGGTAAAAGCGAAGTGCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGACCCGGAAGATCtta                                                                                                                          <  2:180002/251‑3 (MQ=255)
                                                                                                                     gATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCaa                                                                                                                      <  2:207120/249‑1 (MQ=255)
                                                                                                                        cGAAGGTACGTTAGCTGCTGTTCCACAGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGAGCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTg                                                                                                                 <  2:373046/251‑1 (MQ=255)
                                                                                                                             gTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCtt                                                                                                            <  1:132989/251‑1 (MQ=255)
                                                                                                                               aCGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGTCGGTCTGGGTAAAGTGATTTCCCCCCGTGGTGAAAACGGCACCCGCATTGGTTTTGGCGCGGCGGTTAAAGCGAAATCACACCAAGAAAGCCGAAGCGGGCTGCTGGCGCAGGGTGAACCGGACGAACTGCTTAAATTTGGTCTTAt                                                                                                          >  2:386758/1‑251 (MQ=255)
                                                                                                                                     gctgctGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGGAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAAACGGAAGATCTGATCAAGGTTGGTCTTATCCCTg                                                                                                     >  1:284725/1‑250 (MQ=255)
                                                                                                                                        gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGt                                                                                                  <  2:329522/250‑1 (MQ=255)
                                                                                                                                        gctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGt                                                                                                  <  2:279263/250‑1 (MQ=255)
                                                                                                                                                aCCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAAGGGATTTCCCACCGGGGAGAAAACGGCTCCGGCATTGGGTTTGGCGGGACGGGAAAAGCGAAGGCCGGCAAAGCAAGCGAAGGGGAGCTGCTGGCGCCGGGTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGgtc                                                                                         >  1:23070/1‑251 (MQ=255)
                                                                                                                                                          gggctgaattctcagcagcaGGAATGCTTGCAGGTTGATACCTCTAAGATCCTGTTTATGTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCGGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTGATTGGTCGCCTGCCgg                                                                                <  2:366167/236‑1 (MQ=255)
                                                                                                                                                          gggCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCCCCGTGTAGAAACCGGCTCCGGCATTGGGTTTTGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCCGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCggg                                                                               >  1:27347/1‑250 (MQ=255)
                                                                                                                                                           ggCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGGTTGCCGGTCTGGGTAAAGTGATTTTCCACCGTGTAGAAAACGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCGGCTGGCGCAGGTTGAAACGGAAGATGTGGTCAAGGTTGGTCTTATCCCCGAGGTTATTGGTCGTGTGCCGGtt                                                                              >  1:60987/1‑251 (MQ=255)
                                                                                                                                                                   cATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTTGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGGCTTATCCCTGAGTTTTTTGGTCGGCTGCCGGGTGTCGCAAc                                                                      >  1:49279/1‑251 (MQ=255)
                                                                                                                                                                   cATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGCCCGTCTGGATAAAGTGATTTTCCACCGTGTAGAAACCGGCCCCGGCATTGGTTTTGGCGCGCCGGGAAAAGCGGAGTCCGACAAAAGAAGCGAAGGCGAGCTGCTGGCGCCGGGTGAAACGGAAGATCTGATCAAGTTTGGTCTTTTCCCTGAGGTTTTTGGGCGGCTGCCGGGTGTCGCAAc                                                                      >  1:184897/1‑251 (MQ=255)
                                                                                                                                                                      ccGCAGCAGGAATTCTTGCAGGTTGATACCGCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGtt                                                                   <  2:87392/251‑1 (MQ=255)
                                                                                                                                                                       cGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTAGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTtg                                                                  <  1:259463/251‑1 (MQ=255)
                                                                                                                                                                        gcagcaGGAATTCTTGCAGGTTGATTCCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTtga                                                                 <  1:105212/251‑1 (MQ=255)
                                                                                                                                                                           gcagGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCTGGGCTGGAAAAAGTGATTTCCCACCGTGTAGAAAACCGCACCCGCATTGGTTTTGGCGCGACGGGAAAAACGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGAAGGGTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTAATGGTCGTCTGCCGGTTGTCGCAACGTtgaagg                                                              >  2:239118/1‑249 (MQ=255)
                                                                                                                                                                           gcagGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCgg                                                                                                                                                                                                                  <  2:265534/105‑1 (MQ=255)
                                                                                                                                                                                 aTTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTATGGATAAAGTGATTTCCCACCGTGTAGAAACCCGCCCCGGCATTGGTTTTGGCGCGGCGGTAAAAGCGAAGTCCGACAAAACAAGCGAAGGCGAGCTGCTGGCGCAGGTTGGACCGGAAGATCTGGTCAAGTTTGGTCTTATCCCTGAGGTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTg                                                         >  1:398919/1‑250 (MQ=255)
                                                                                                                                                                                        cAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGGTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCgaaga                                                 <  2:127672/251‑1 (MQ=255)
                                                                                                                                                                                          ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGGCCGGGTGAAGGGATTTTCCACCCGGTAGAAACCGGCACCGGCATTGGTTTTGGCGCGGCGGGAAAAACGAAGTCCCGCAAAGCAAACGAAGGCGAGCGGCTGGGGGAGGGTGAACCGGGAAAACTGGACAAGTTTGGTCTTATCCCTGAGGTTTTTGGTCGGCTGCCGGGTGTGGCAACGATGGATGAACTGAGCgacgaa                                                >  1:6876/1‑250 (MQ=255)
                                                                                                                                                                                           gTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGCCCGTCTGGATAAAAGGATTTCCCACCGGGTAGAAACCGGCCCCGGCATTGGTTTTGGCGCGACGGTGAAAGCGAAGGCCGACAAAGCAAGCGAAGGGGAGCCGCTGGCGCCGGTTGAACCGGCAGAGCTGATCAAGTTTGGTCTTCTCCCCGGGGTTATTTGTCGTCTGCCGGTTGTCGCCACGTTGGATGAACTGGGCgcagcagc                                              >  2:70175/1‑244 (MQ=255)
                                                                                                                                                                                                                  taattggggCGGTGCGTTTGCCGGTCTGGATAAAGTGCTTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCTTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAagag                        <  2:268972/243‑1 (MQ=255)
                                                                                                                                                                                                                   aTTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGc                       <  2:252874/250‑1 (MQ=255)
                                                                                                                                                                                                                      tgtgGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAaagcaagc                                                                                                                                                                   <  1:121531/108‑1 (MQ=255)
                                                                                                                                                                                                                       gtgGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAATAGCCGaaa                  <  1:406081/251‑1 (MQ=255)
                                                                                                                                                                                                                            ggTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTAAGATCCTCAAAGAGCCGAAAAACGc             >  2:165618/1‑251 (MQ=255)
                                                                                                                                                                                                                             gTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGcc            >  2:286866/1‑251 (MQ=255)
                                                                                                                                                                                                                                          gTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGc  >  2:169723/1‑249 (MQ=255)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                                      
AGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAAGTGATTTCCCACCGTGTAGAAACCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAGTTTGGTCTTATCCCTGAGTTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGC  >  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: