Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913691,8830CT26.5% 88.2 / 11.0 50E153E (GAG→GAA) ybeYendoribonuclease YbeY
Reads supporting (aligned to +/- strand):  ref base C (12/24);  new base T (12/1);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 2.87e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGG  >  NC_000913/691634‑692131
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          
cgATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTAtt‑ctttc                                                                                                                                                                                                                                                            <  1:331464/250‑1 (MQ=255)
 gATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGCAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTCCTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGGTAGGT‑TGTCCGG‑CCAATGACTGGCCAATTATT‑CTTtctc                                                                                                                                                                                                                                                          >  2:280775/1‑251 (MQ=255)
 gATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTtctc                                                                                                                                                                                                                                                          <  1:174684/251‑1 (MQ=255)
 gATAAGGTCGTGCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGGGTGAATTGTCGTCGCGCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTTGCAAATTATT‑CTTtctc                                                                                                                                                                                                                                                          <  1:495428/251‑1 (MQ=255)
    aaGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCgg                                                                                                                                                                                                                                                        <  1:349535/250‑1 (MQ=255)
    aaGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGc                                                                                                                                                                                                                                                       <  1:479756/251‑1 (MQ=255)
    aaGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGATGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCgg                                                                                                                                                                                                                                                        <  1:82045/250‑1 (MQ=255)
             cTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑GCGTCTGCCAGT‑CAAAAGCCCGGCCAATTATT‑TTTTTTCGGCCATGTACgg                                                                                                                                                                                                                                              >  2:556599/1‑251 (MQ=255)
               gCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACgg                                                                                                                                                                                                                                              <  1:206207/249‑1 (MQ=255)
               gCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACg                                                                                                                                                                                                                                               <  1:20072/248‑1 (MQ=255)
               gCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAt                                                                                                                                                                                                                                            <  1:209650/251‑1 (MQ=255)
               gCCCGGAATCACGGATCAGCGCCAGCAGTGCGTCACGGTTTGTCGGTTCACCGTGGAAAAATTGGTGGAGTAACAGGGAGAAAAATCCCTTCTTTTTGCGTAGGTTGGCACTACGGTGTGAATTGTCGTCGCTCATGGCTTCGGATGTGTTCTCATGTTAGTGAATGAAAACGCCGTCGGTGAGCACCAACTGCAGGC‑GCGTCGGCCAGT‑CAAAGTCGTGGTAAATGATG‑CTTTCTCTGCAATGTACgg                                                                                                                                                                                                                                              <  2:618586/249‑1 (MQ=255)
                        cACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑GCGTCTGCCCGT‑CCAATGACTGGCAAATTATT‑CTTTTTCCGCAATGTACGGAT‑CCT‑CCTAGc                                                                                                                                                                                                                                 >  1:243891/1‑251 (MQ=255)
                                     cagcagTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCATAGCAAGCAt                                                                                                                                                                                                                    <  1:549517/251‑1 (MQ=255)
                                              gTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGt                                                                                                                                                                                                           <  1:99198/251‑1 (MQ=255)
                                                       tttCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAggg                                                                                                                                                                                                   <  1:32518/250‑1 (MQ=255)
                                                       tttCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAggg                                                                                                                                                                                                   <  1:280775/250‑1 (MQ=255)
                                                        ttCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTTTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGGCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATGATT‑CTTTCTCGGCTATGTGCGAAG‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGc                                                                                                                                                                                                  <  1:326852/250‑1 (MQ=255)
                                                           ggTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGGCTGCCAGT‑CAAATGCCTGGCAAATTATT‑TTTTTTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCCTAATCTCTTTTTCCAGGGCTTcc                                                                                                                                                                                              >  1:212349/1‑251 (MQ=255)
                                                                  cGTGGAAAAGTTGGCTGAGTATCAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCGTCCATttctt                                                                                                                                                                                       <  1:172931/251‑1 (MQ=255)
                                                                      gAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTAATGTGTGAATTGTCGTCGCTAATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑TCGCCCTCCCGTACAAA‑GCCCGGCACAATATT‑TTTTTTCGGCCACGTTCGGAT‑CCTACACA‑CCCAGAGCAAGCATACTCTCTTTTTCGAGGGCTTTCATttctttttc                                                                                                                                                                                   >  2:78150/1‑247 (MQ=255)
                                                                        aaaGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑GCGTCCGCCATT‑TAAATATCCGGCATACTATT‑CTTTTTTGGCGACTTACGGAT‑CCT‑CCAAGCCCAGAGCAAGCACAATCGCTGTTTCGAGGGGTTCCCTTTCTTCTGCtt                                                                                                                                                                                 >  2:380395/1‑251 (MQ=255)
                                                                              ggCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCCGT‑CCAATTCCTGGCCAATTATT‑CTTTTTCGGCAATGTACGGGT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTTTTCTGCTTCGTCAt                                                                                                                                                                           >  1:373255/1‑251 (MQ=255)
                                                                                    gtaacagggaccacactCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑GCGTCTGCCCGT‑CCAATGCCTGGCGCATTTTT‑TTTTCTCGGCGATGTAACGAT‑CCT‑CCTAACCCCGAGCCAGCCTAATCTTTTTTTCGAGGGCTTCCCTTTCTTCTGCTTCTTCCTCtttcc                                                                                                                                                                     >  2:276706/17‑248 (MQ=255)
                                                                                       aCAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑GCGGCTGCCAGG‑CAAAAGCCTGGCAAATTATT‑TTTTCTCGGCAATGTACGGAC‑CCT‑CCTAGCCCAGAGCAAGCATAATCTCTGTTTTGAGGGCTTCCATTTCTTCTGCCTCGTCATTTTCggtgg                                                                                                                                                                  >  1:312329/1‑247 (MQ=255)
                                                                                           ggAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑TTTTTTCGGCCATGTACGGAT‑CCT‑CATAGCCCAGAGGCAGCATAATCTCTTTTTCGAGGGCTTCCATTTCTTTTGCTTCGTCATCTTCGATGTGATc                                                                                                                                                              >  1:28759/1‑251 (MQ=255)
                                                                                                    tCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCCGT‑CAAATGCCTGGCAAATTTTT‑CTTTTTCCGCCATGTTCGGCT‑CCT‑CCTAGCCCCGAGCAAACATAATTTTTGTTTCGAGGGCTTCCATTTTTTCTGCTTCGTCATCTTCGATGTGGTCGTAACCtt                                                                                                                                                      >  1:185157/1‑249 (MQ=255)
                                                                                                       cttcttGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCCTGGAGTCGTATGGGTTCTCATGTTAGTTCATCAACACGCAGTCGTTAATCACCAAAGGCGGGG‑CAGTCTGCAAGA‑CCAAAACCAGGCCAATTATTTCTTT‑TCCGCACTGTACGGGTACCT‑‑ATAGCCCACAGGAAGCATCATCTCCGTTTCGAGGGCTTCCCTTTCTTTTCCTTCGCCATCTTCGATGTGATCGTAACCTaaaa                                                                                                                                                   >  1:285638/1‑248 (MQ=255)
                                                                                                          cttGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑GCGTCTGCCCGT‑CAAAAGCCTGGGAAATTATT‑CTTTTTTGGCAATGTACGGAT‑CCT‑CCTAGCCCAGAGCAAGCATAATCTCTTTTTTGAGGGCTTCCCTTTTTTTTTCTTCGTCCTCTTCTATGTGATCGTAACCTTACAAAtt                                                                                                                                               >  1:53121/1‑250 (MQ=255)
                                                                                                                  ttATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑GCGTCTGCCCGT‑CAAATGCCTGGCAAATTATT‑TTTTTTCGGCCATGTACGGAT‑CCT‑CCTAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCCTTTCTTCTTCTTCGTCATCTTCGTTTTGGTCGTAACCTAAccaaatcagaattc                                                                                                                                       >  2:9810/1‑238 (MQ=255)
                                                                                                                       gtgtCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGTCAAT‑CAAATGCCTGGCCAATTATT‑TTTTTTCGGCCAAGGACGGAT‑CCT‑CATTGCCCAGAGCAAACAAAATTTCTGCTTCGAGGGGTTCCATTTTTTCTGCTTCGTCATCTTCGTTGTTATCGTAACCTAACAAAATCAGACTGCCGTGc                                                                                                                                  >  2:190173/1‑251 (MQ=255)
                                                                                                                            actactGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑TTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCCGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCCTTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAAAAAATGCAGACTGCCGTGcaccac                                                                                                                             >  1:202719/1‑251 (MQ=255)
                                                                                                                                              gtcgCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTgcgc                                                                                                           <  2:28759/251‑1 (MQ=255)
                                                                                                                                                    cATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTcc                                                                                                       <  1:463338/249‑1 (MQ=255)
                                                                                                                                                    cATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTGCTGCTTCGTCATCTTCGATGTGATCGTAACCGAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAg                                                                                                     <  1:619042/251‑1 (MQ=255)
                                                                                                                                                     aTGGCGTCGTATGGGTTCTCATGATAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGCG‑AACCGTGCCAGG‑CAAAAGCCCGGCAAAATAAT‑CTTTTTCGGAAATGTACGGAT‑TCT‑CATAGTCCAGAAGAAGCGTAATCTCTGTTTTTAGGGCTTCCATTTCTTCTCCTTCGTCCTCTTCCGTTTTATCGCAACCCAACAAAACAAAACTGCCGTACACCACCCTATGCGCCCAGCGGGCCCCCAGt                                                                                                    >  1:544569/1‑251 (MQ=255)
                                                                                                                                                       ggcatcgtATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTGCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGTGCCTCCAGTg                                                                                                   <  1:445153/246‑1 (MQ=255)
                                                                                                                                                                   gTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAAGGTACGGAT‑CCT‑GATAGCCCAGAGCAAGCATAATCTCTGTGTCGAGGGCTTGCATTTCTTCTTCTTCGTCAGCTTCGATGTGATCGTAACCGAACAAATGCAGGGTGCCGGGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTc                                                                                      <  1:304643/251‑1 (MQ=255)
                                                                                                                                                                                  aaTCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTtctc                                                                       <  2:373255/251‑1 (MQ=255)
                                                                                                                                                                                  aaTCAAAACGCCGTCGTGAATCACCAACGGCGGGG‑ACGTCGGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCAGTTCTTCTGCGTCGTCATCTTCGATGTGATCGTAACCTAACAAATGc                                                                                                                                              <  2:245624/180‑1 (MQ=255)
                                                                                                                                                                                      aaaaCGCCGTCGTTAATCACCAACGGCGGCG‑GCGTCTGCCAATACAAA‑GACCGGCAAATTATT‑CTTTCTACGCCATGTAACGTC‑CCT‑CACATCCCAAAGCAAACATAATCTCTGTTTCGAGGGCTTCTTTTTTTTCTGCTTCGTCATCTTCTAAGGGATCAGCACCTAACAAAAGCAGACTGCCCTTCACCACCCTATGCGCCCAGTGCCCCACCCGTGGTTTGCCTTGCTCCTGGGcttacctcccaacc                                                                   >  1:676792/1‑240 (MQ=255)
                                                                                                                                                                                           gCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGcttccttc                                                                         <  1:549933/240‑1 (MQ=255)
                                                                                                                                                                                            ccGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGG‑CCAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTTGAGGGCTTCCCTTTCTTTTGCTTCGTCATCTTCGCTGTGATCGGAACCTTACCAATTCAGACTGCCGTGCACCCCCACATGCCCCCCGTGCGCCTCCAGTGGGTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTg                                                              >  2:209403/1‑250 (MQ=255)
                                                                                                                                                                                             cgtcgtTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGAc                                                            <  2:53121/251‑1 (MQ=255)
                                                                                                                                                                                              gtcgtTAATCACCAACGGCGGGG‑ACGTCTCCCCGT‑CCAATGCCTGGCAAAATTTT‑TTTTCTCGGCAAAGTACGGGT‑ACT‑CATAACCCCGAGCAAGCCTAATCTCTGTTTCGAGGGCTTCCCTTTCTTCTGCTTCGTCATCTTCGCTGTGATCGTAACCTAACAAAAGCAGACTGCCGTGCACCACCCTATGCGCCCAGTGCGCCTCCAGGGGTTTGCCCTGCTCCTGAGCTTCCTTCTCAACCACCTGACg                                                           >  1:217343/1‑251 (MQ=255)
                                                                                                                                                                                                cgtTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑GCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACgg                                                          <  1:380395/250‑1 (MQ=255)
                                                                                                                                                                                                             cggcggGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTGCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGGTTGTCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATc                                            <  1:652309/251‑1 (MQ=255)
                                                                                                                                                                                                                      aCGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGt                                     >  2:173577/1‑250 (MQ=255)
                                                                                                                                                                                                                      aCGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGt                                     >  2:96927/1‑250 (MQ=255)
                                                                                                                                                                                                                          cTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCg                                 >  2:124255/1‑250 (MQ=255)
                                                                                                                                                                                                                             ccAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCCGTAGCGACAt                             >  2:585277/1‑251 (MQ=255)
                                                                                                                                                                                                                                       tGCCTGGCAACTTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGCTCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGcc                    >  2:268827/1‑251 (MQ=255)
                                                                                                                                                                                                                                         ccTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCAAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGcc                    >  1:535238/1‑249 (MQ=255)
                                                                                                                                                                                                                                             gCAAATTATT‑CGTGTTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCAGAATCTCTGTTTCGAGGGCTTCCATTTCGTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAggcgg              <  2:189363/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           tctcGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACgg  <  1:227409/250‑1 (MQ=255)
                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                          
CGATAAGGTCGTTCTGCCCGGAATCACGGATCAGCGCCAGCAGTTCGTCACGGTTTTTCGGTTCACCGTGGAAAAGTTGGCTGAGTAACAGGGAGAAAAATCCCTTCTTGTTGCTTATCGTGTCACTACTGTGTGAATTGTCGTCGCTCATGGCGTCGTATGGGTTCTCATGTTAGTTAATCAAAACGCCGTCGTTAATCACCAACGGCGGGG‑ACGTCTGCCAGT‑CAAATGCCTGGCAAATTATT‑CTTTCTCGGCAATGTACGGAT‑CCT‑CATAGCCCAGAGCAAGCATAATCTCTGTTTCGAGGGCTTCCATTTCTTCTGCTTCGTCATCTTCGATGTGATCGTAACCTAACAAATGCAGACTGCCGTGCACCACCATATGCGCCCAGTGCGCCTCCAGTGGTTTGCCTTGCTCCTGAGCTTCCTTCTCAACCACCTGACGGCAGATAACCAGATCGCCCAGTAGCGACATTTCCATGCCAGGCGGCACTTCAAACGG  >  NC_000913/691634‑692131

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: ‑