Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913458,1890AG22.6% 105.3 / 11.3 54D255G (GAT→GGT) clpXATP‑dependent Clp protease ATP‑binding subunit ClpX
Reads supporting (aligned to +/- strand):  ref base A (11/30);  new base G (12/0);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 5.07e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGC  >  NC_000913/457940‑458431
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
aTGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGAt                                                                                                                                                                                                                                                           <  1:334030/251‑1 (MQ=255)
  gTCCAGAAAGCTCAGCGTGGTAGTGTCTTCTTCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATa                                                                                                                                                                                                                                                          <  2:590669/250‑1 (MQ=255)
  gTCCAGAAAGCACCGCGTTGTATTGACTACATCGAGGAAATCGACAAGATTTCTCGTAAGGCAGACAACCCTTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGAGCGAAGGTACGGTAGCTGCTTTTCCGCCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATGCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCTGGTCTGGATa                                                                                                                                                                                                                                                          <  2:309240/250‑1 (MQ=255)
       gAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATTCCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGa                                                                                                                                                                                                                                                   <  1:327016/251‑1 (MQ=255)
          aGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACTGTTCCGGTGAAGTCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGAttt                                                                                                                                                                                                                                                <  2:674395/251‑1 (MQ=255)
                  gTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCgt                                                                                                                                                                                                                                       <  1:442378/251‑1 (MQ=255)
                    ggTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTTCCGGGCTGGGTAAA‑GTGATTTTCCCAcgtgt                                                                                                                                                                                                                                      >  1:310300/1‑247 (MQ=255)
                          gTCTACAGCGATGAAACCGACAAGATGTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGAGCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑c                                                                                                                                                                                                                           <  2:299559/251‑1 (MQ=255)
                              aCATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCCGGGCTGGGTAAA‑GTGATTT‑TCCACCGGG‑TAAG‑AAA‑cggc                                                                                                                                                                                                                        >  1:504786/1‑248 (MQ=255)
                                      gAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTTATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCAt                                                                                                                                                                                                               <  1:654719/251‑1 (MQ=255)
                                                       tcGTAAGTCAGACAACCCGTCCATTCACCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTTTTGAAACTGAGCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCTTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCTACg                                                                                                                                                                                               <  2:607674/250‑1 (MQ=255)
                                                                   cAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCCGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGGCGGTATGGGGAAA‑TGGATTT‑CCCACCGTGCTG‑G‑AAA‑CGGCCTCCG‑CATTGGTGTTgtggcgtcggcaaaagagaac                                                                                                                                                                                   >  2:249182/1‑230 (MQ=255)
                                                                    aaCCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGTCGGGCTTGGTAAA‑AGGATTT‑TCCCCCGTG‑TA‑G‑AAA‑CAGGCTCCCGCATTGGtttttgtgggcaggtaacagagaagtc                                                                                                                                                                                 >  2:567434/1‑228 (MQ=255)
                                                                       ccGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGa                                                                                                                                                                              <  1:567434/251‑1 (MQ=255)
                                                                        cGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGACACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGAAGGGCGTAAACATCCGCAGCAGGAATTCTTGAAGGTTGATACCTCTAAGACCATGTTTATTTGTGGCGGTGCCGTTGCCGGTCGGGGTAAA‑GAGATTT‑CCCGCCCTT‑AA‑G‑AAA‑CATGGTCCCGCATTGGTTTTGGAGCTACgggtaaagagaaacacgac                                                                                                                                                                             >  1:80252/1‑234 (MQ=255)
                                                                         gTCCATTACCCGAGACGTTTCCGGTGCAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGCTCCTGTTTATTTGTGGCGGGTCGTTTGCCGGTCTGGGAAAA‑ATGAATT‑CCCCCCCCG‑GA‑G‑AAAACCGCCTC‑GGCCCTTGTTTTGGTGCGCCGGTAAAAGCGAAgagcggc                                                                                                                                                                             >  1:573517/1‑244 (MQ=255)
                                                                           ccATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGGTAAA‑GGGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGGCGGTAAAAGCGAAGTCCGACaa                                                                                                                                                                           >  1:255493/1‑250 (MQ=255)
                                                                                            tCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGCTCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGGCCGGATAAA‑GTGATTT‑TCCCCCGTG‑TT‑G‑AAA‑CCGGCTCCGGCCTTTGTTTTGGCGCGACGGTAAAAACGGAGGCCGACAAAACAAGCGGAGGCgcgcg                                                                                                                                                         >  1:161514/1‑247 (MQ=255)
                                                                                              cGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGACGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTTCCCGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCCTTGGTTTTTGCGCGACGGTAAAAACGAAGTCCGACAAAGCAAGCGAAGGCGAgcagg                                                                                                                                                       >  1:271503/1‑248 (MQ=255)
                                                                                              cGGTGAAGGCGTACAGCAGGCACTGGTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCGTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgc                                                                                                                                                          <  2:271503/248‑1 (MQ=255)
                                                                                               ggTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAgctgct                                                                                                                                                      <  1:668033/251‑1 (MQ=255)
                                                                                                             gcagGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCa                                                                                                                                                <  2:15173/243‑1 (MQ=255)
                                                                                                             gcagGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGCTTGCCGGTCTGGGTAAA‑GGGATTT‑CCCCCCGTG‑TT‑G‑AAA‑CCGGCTCCGGCCTTGGTTTTTTCGCGGCGGTAAAAGCGAGGGCCGGCAAAGCAAGCGAAAGGGAGCTGCTggggcgagatcggg                                                                                                                                        >  1:15173/1‑239 (MQ=255)
                                                                                                              cagGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAATGGGTTTT‑CCC‑CCCTG‑TA‑G‑AAA‑ACGGCGCCGGCGTTGGTTTTGGCGCGACGGGAAAAACGAAGTCCGCCACAACAAACGAAGGGGAGGTGCTGGCGCAGGTTGaaac                                                                                                                                       >  1:530834/1‑249 (MQ=255)
                                                                                                              cagGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGCCGGTCTGGGTAAA‑GGGGTTT‑TCCCCCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGGCGGGAAAAGCGAAGTCCGACAAAGGAAACGAAGGGGAGCGACTGGGGCAGGGTGAAcc                                                                                                                                       >  1:83562/1‑251 (MQ=255)
                                                                                                                           aaCTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATc                                                                                                                          <  1:210564/251‑1 (MQ=255)
                                                                                                                           aaCTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATc                                                                                                                          <  1:44004/251‑1 (MQ=255)
                                                                                                                               gATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAg                                                                                                                       <  1:510590/250‑1 (MQ=255)
                                                                                                                                   gAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTgg                                                                                                                 <  1:615103/251‑1 (MQ=255)
                                                                                                                                   gAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTTCCGGGCTGGGTAAA‑ATGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTTGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGGAAGCGGGCTGCTGGGGCAGGTTGGACCGGAAGATCTGATCAAG‑TTTgg                                                                                                                 >  1:621461/1‑251 (MQ=255)
                                                                                                                                       gtcaggTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCtt                                                                                                             <  2:253267/247‑1 (MQ=255)
                                                                                                                                       gTACGGTAGCTGCTGTTCCACCGCAAGGTGTGCGTAAACATCCGCCGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGTGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAACGCGAAGTCCGCCAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCtt                                                                                                             >  2:328470/1‑251 (MQ=255)
                                                                                                                                        tccGGTAGCTGCTGTTCCACCGCAAGGAGGACGTAAACATCCGCAGCCGGACTTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGGGTTTTCCGGCCTTGGTAAA‑GTCATTT‑CCCCCCGTG‑TT‑G‑ACA‑CCGGCTCCGGCGTTCGTTTTGGCGGGGCGGGAAAACCGAAATCCGGCCAAGCAAGGGGAGGCGAGGTTCTGGTGCCGGGTGAACCGGAAGATCTGATccaatttcggcttt                                                                                                             >  1:621939/3‑238 (MQ=255)
                                                                                                                                             aaGCTGGTGGTCCACCGCAAGGTGGGCTTAAACATCCGCAGCAGGAATTCTGGCCGGTTGATAACGCTAAGTTCGTGGTTATTTGTGGCGGTGCGTTTGTCGGGCTGGATAAA‑GTGAGTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCTGCATTGGCTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAATGCTAGCTGCTGGCGCAGGTCGAACCGGAAGTGCTGATCAAG‑TTTGGTCTTATCCCt                                                                                                       <  2:690166/250‑1 (MQ=255)
                                                                                                                                                ctgctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCCTTTGACGGTCTGGAAAAA‑GGGATTT‑CCCAACCTG‑GA‑G‑CAA‑ACGGGTCAGGCATTGGGTTTAGCGCGAAGGGAAAACCGAGGACCGAAAAAGGAAGCCGAGGAGAGCTGCTGGGGCAGGTTGAAACGGCAGATCCGGTCAAG‑ATTGGTCTTTTCCCTg                                                                                                      >  1:31045/1‑249 (MQ=255)
                                                                                                                                                   ctgTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGgc                                                                                                                                                                                                     <  1:709903/152‑1 (MQ=255)
                                                                                                                                                        ccACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGGCGGTCTGGGTAAA‑GTGATTT‑CCCACCGTG‑GA‑G‑AAA‑CCGGCTCCGGCCTTGGTTTTTGCGCGACGGGAAAAACGAAGACCGACAAAGCAAGCGAAGGCGCGCTGCTGGCGCAGGTTGAACCGGAAGCGCTGATAAAG‑TTTTGTCTTTTCCCTGCG‑TTTTTTgg                                                                                           >  1:191579/1‑251 (MQ=255)
                                                                                                                                                           ccGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGgtcg                                                                                        <  1:349675/251‑1 (MQ=255)
                                                                                                                                                                aGGTGGGAGTACACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTCTATTTGTGGCGGCGGGGTTGCCGGGCTGGATCAA‑GAGATTT‑CCCACCATG‑GA‑GGAAA‑CTGG‑TCCGGCACTGGTGTTTGCGCGGCCGTAAAAGCGGAGTACGACGAAACAAGCGAAGGGGAGCTGGTGGCGCAGGGTGAACCCGAAGATATGATCATGATTTGGTCT‑ATACCTGAG‑TTTAttagttgttttg                                                                                   >  2:307389/1‑241 (MQ=255)
                                                                                                                                                                                ccGCAGCAGGAATTCGTGCTGGTTGATACCTCTAAGTTCCTGTGTATTGGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGAGTT‑CCCCCCGTG‑GA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCGGAG‑TTTATTGGTCGTCTGCCTGTTGTCGCAACGt                                                                    <  2:319164/250‑1 (MQ=255)
                                                                                                                                                                                                  gaGGTTGAGACCTCTAAGAGCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGAGTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACTGTAAAAGCGAAGTCCGACAAAGCAAGCGAATGCGAGCTGCGGGCGCAGGTTGAACCGGAAGATCGGATCAAG‑GTTGGTCTTAGCCCTGAG‑TTTATTTGTCGGCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCgaaga                                                 <  2:310300/250‑1 (MQ=255)
                                                                                                                                                                                                    ggTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGGTTGGCGGTCTGGATAAA‑ATGGTTT‑CCCACCGTG‑GA‑G‑GAA‑CCGGCTCCGGCCTTGGTTTTTGCGCGCCGGTAAAAACGAAGTCCGACAAAGCAAGCGAAGGCGGGCTGCTGGCGCCGGGTGGACCGGAAGCACTGATAAAG‑TTTGGTTTTTTCCCTGGGGTTTTTTGACGGTCT‑CCGGTGGTCGCAACGTTTAATGAAcaggggggaaa                                                 >  2:449428/1‑239 (MQ=255)
                                                                                                                                                                                                      ttGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGCTAAA‑GTGATTT‑CCCACCGTG‑GA‑G‑AAA‑CCGGCTCCCGCATTGGGTTTTGCGCGACGGGAAAAGCGGAGTCCGACAAAGCAAGCGAAGGCGAGGTGCTGGGGCCGGTTGAACCGGAAGATATGATAAAG‑TTTGGTCTTATTCCTGAG‑GTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATTAACTGAGGgaagaaact                                             >  1:449946/1‑248 (MQ=255)
                                                                                                                                                                                                       tGATCCCTCTAAGTGCCTGTTTATTTGTGTCGGTGCGTTTGCCGGGCTGGATAAA‑GTGAGTT‑CCCTCCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TGTATTGGTCGGCTGCCTGTTGTCGCAACGCTGAATGAACTGAGCGAAGAAGct                                             <  2:263709/250‑1 (MQ=255)
                                                                                                                                                                                                               ctAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTc                                      <  1:171804/249‑1 (MQ=255)
                                                                                                                                                                                                                tAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGa                                                                                                                                                            <  1:194107/132‑1 (MQ=255)
                                                                                                                                                                                                                 aaGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGag                                  <  2:255493/251‑2 (MQ=255)
                                                                                                                                                                                                                                 gtgGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGaaa                  <  1:449428/251‑1 (MQ=255)
                                                                                                                                                                                                                                  tgGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGaaa                  >  1:241068/1‑250 (MQ=255)
                                                                                                                                                                                                                                      ggTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGc             >  1:380525/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCTGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGc  <  1:335849/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGc  >  1:607869/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   gTCTTGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGTAAGATc                                                                                                                               <  1:82222/126‑1 (MQ=255)
                                                                                                                                                                                                                                                   gTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑CCA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAACCGCGAAGTCCGACAAAGCACGCGAAGGCGAGCTGCTGGCGCAGGTTGACCCGGAAGATc                                                                                                                               >  2:82222/1‑126 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          
ATGTCCAGAAAGCACAGCGTGGTATTGTCTACATCGATGAAATCGACAAGATTTCTCGTAAGTCAGACAACCCGTCCATTACCCGAGACGTTTCCGGTGAAGGCGTACAGCAGGCACTGTTGAAACTGATCGAAGGTACGGTAGCTGCTGTTCCACCGCAAGGTGGGCGTAAACATCCGCAGCAGGAATTCTTGCAGGTTGATACCTCTAAGATCCTGTTTATTTGTGGCGGTGCGTTTGCCGGTCTGGATAAA‑GTGATTT‑CCCACCGTG‑TA‑G‑AAA‑CCGGCTCCGGCATTGGTTTTGGCGCGACGGTAAAAGCGAAGTCCGACAAAGCAAGCGAAGGCGAGCTGCTGGCGCAGGTTGAACCGGAAGATCTGATCAAG‑TTTGGTCTTATCCCTGAG‑TTTATTGGTCGTCTGCCGGTTGTCGCAACGTTGAATGAACTGAGCGAAGAAGCTCTGATTCAGATCCTCAAAGAGCCGAAAAACGCCCTGACCAAGC  >  NC_000913/457940‑458431

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