Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,133,3690CT24.0% 111.5 / 15.0 50E197K (GAA→AAA) yghSputative ATP‑binding protein YghS
Reads supporting (aligned to +/- strand):  ref base C (4/34);  new base T (12/0);  total (16/34)
Fisher's exact test for biased strand distribution p-value = 1.50e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCTT  >  NC_000913/3133122‑3133608
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
cATCACCCCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTg                                                                                                                                                                                                                                                  <  2:91693/251‑1 (MQ=255)
   cacaCCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑aa                                                                                                                                                                                                                                              <  1:453727/251‑1 (MQ=255)
              accGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGTTCCAGGATCTTTTCGCCCTTTAATG‑TCAGTTGCGGCGGAACGGGGATTTTTTTCTGT‑AACGCTGCCAGtt                                                                                                                                                                                                                                   >  1:58668/1‑251 (MQ=255)
                    agcAGCCCCGGAATGGAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTggtgg‑t                                                                                                                                                                                                                            <  2:498182/249‑1 (MQ=255)
                      cAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑Tc                                                                                                                                                                                                                           <  2:232647/250‑1 (MQ=255)
                      cAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGGGCAGCCATTCGTTAAATCAGCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑Tc                                                                                                                                                                                                                           <  2:242938/250‑1 (MQ=255)
                                                   cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTtcatc                                                                                                                                                                                            <  2:58668/251‑1 (MQ=255)
                                                   cgGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTTGCCCCTCCCGTTCAAGGATCTTTGGGCCCTTGTATT‑TCCGTTGGGGCGTAACGGCGGTTTTTTTCTGT‑TACGCCTCCAGTTTGTGG‑GTGGGTTTACGCGCAAACCCGG‑tgtttttatt                                                                                                                                                                                            >  2:4809/1‑242 (MQ=255)
                                                          aCTATTAACAGTTTGAGTTTGTAGTGGATCCATTCGTTAAATCATCCCCTATCCATTTTGGTTATCAGGAAAGCGCGGCGTGAATTCCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTTG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCa                                                                                                                                                                                     <  2:119514/251‑1 (MQ=255)
                                                           cgtttAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCaa                                                                                                                                                                                    <  2:35786/248‑1 (MQ=255)
                                                             aTTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTTCGCCATTGAATG‑TCAGTTGCGGCGGAACGGCGATTTTTTTCTGT‑AAAGCTGCCCGTTGGTGG‑TCAGGGTTACGCGCAAACGAGG‑GTTTTTCATCAATGCCa                                                                                                                                                                                     >  1:316896/1‑248 (MQ=255)
                                                                aaCAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGAGCTTTTCGCCTTTGAATG‑TCAGTTGCGTCGTAACGGCGATTTTTGCCTGT‑AACGCTTCCAGTTGGTGG‑TCAGTTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACg                                                                                                                                                                                <  1:596656/250‑1 (MQ=255)
                                                                    ggttgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCTTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCTGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTTGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCATGACGAA‑Tca                                                                                                                                                                          <  1:379007/249‑1 (MQ=255)
                                                                       tgagtttgTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCATACAGTTCCAGGAGTTTTTCGGCATTTAAAT‑TCAGGTGCGGCGTAAAGGCGATTTTTTTCTTT‑AAAGCTGGCCGTTGGGGG‑GTAGGTGTACGCGCCAAAACGG‑TTTTGTCATTAATGCCAAGACGAa‑acaac                                                                                                                                                                        >  1:130809/1‑250 (MQ=255)
                                                                              gTAGTGCTTCCATTCGTTAAATCATCCCCGATCCAGTTTGGTGATCAGGAAAGCGCGGCGTGTGTTGCGCGTAGTGACGCTTGCAGAATTTCAGCGGCGGGGTGTCGCCCATCCAGTTCAAGGAGTTTTGCGCCAGTGAATG‑TCAGTTGCGGCGTAACGGCGATTGTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATTCggg                                                                                                                                                                <  2:427068/251‑1 (MQ=255)
                                                                                   gCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCaaa                                                                                                                                                            <  1:575347/250‑1 (MQ=255)
                                                                                         aTTCGTCAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCAGTTCAAGGATCTTTGCGCCCTTGAATG‑TCAGATGCGGCGGAACGGCGATTTTTTTCTTT‑AACGCTGTCAGTTTGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCCCGACGAA‑TCAACAATAAGGGGAAATAAAATg                                                                                                                                                     >  1:154438/1‑251 (MQ=255)
                                                                                            cGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCAGTTTAAGGAACTTTGCGCCATTTGATG‑TCAGTTGCGGCGGAAAGGGGAATTTTTTCTGT‑AAAGCTGGCCGTTGGGGGGTC‑GGTTTAAGCGCCAAAGCGG‑GTTGTTCATTAATTCCAAGACGAA‑ATAACAAAACAGGCAAATAAGATGcc                                                                                                                                                   >  1:550542/1‑250 (MQ=255)
                                                                                             gTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGGCGGCCATCCAGTTCAAGGATCTTTGCGCCCTTTGATG‑TCAGTTGCGGCGTAACGGCGATTTTTTTCTTT‑AACCCTGCCCGGTGGGGG‑TCAGGTTTTCGCGCAAACCCGG‑GTTGTTTTTCCATTCCAAGACGAAATCAAAAAAACGGGaaaaaaaaagcca                                                                                                                                                   >  1:402061/1‑241 (MQ=255)
                                                                                               tAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGccatc                                                                                                                                                <  1:597797/250‑1 (MQ=255)
                                                                                                 aaTCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGccatccat                                                                                                                                             <  1:224822/251‑1 (MQ=255)
                                                                                                  atcatcCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTAGCCCATCCCGTTCAAGGGAATTTTCGCCCCTTAATG‑TTAGTTGCGGCGTAAAGGCGATTTTTTTCTGT‑AACGCTGCCAGTTGGGGG‑GCAGGGTTATGCGCAAACGCGG‑GTTTTTTATCAATGCCCAGACGAA‑ACAAAAAAACGGGCAAATAAAATGCCACCCAtt                                                                                                                                            >  1:585825/1‑251 (MQ=255)
                                                                                                       ccccTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTAc                                                                                                                                       <  1:51872/251‑1 (MQ=255)
                                                                                                                aTTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCATCCCGTTCAAAGATCTTTGCGCCATTGATTG‑GCAGTTTCGGGGTAACGGGGCTTTTTTTCTGT‑AACGCTTCCAGTTGGGGG‑GTAGGTTTACGCGCCAACGCGG‑GTTGTTCATCAATGCCCAGAAGAA‑ACAACAAAAAGGGGAAAAAAAAAGCCATCCATTGGTaaacctttaac                                                                                                                              >  1:105964/1‑241 (MQ=255)
                                                                                                                  tttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGCCCCTCCCGGTCAAGGGTCTTTGCGGCATTGAATG‑TCAGGTGCGGCGGAACGGCGGTTTTTTTCTGT‑AACGCTGCCCGGTGGTGG‑TGAGGTTTAAGCGCAAAAGCGG‑GTTGTTTATTAATGGCAAAAAGGA‑TCAACAAAAAGGGGAAATAAAATGGCATCCATTGGTACAACtttcgctc                                                                                                                            >  1:107704/1‑245 (MQ=255)
                                                                                                                   ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTc                                                                                                                            <  1:402422/250‑1 (MQ=255)
                                                                                                                   ttGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGGCGCCCCTCCCGCTCAAGAATTTTTTCGCCCCTGAATGGTC‑GTTGCGGCGGAACGGCGGTTTTTTTCTGTCAAC‑CTGGCCGCTGGTGG‑TCAGGGTTCCGCGCAAAACCGG‑GTTGTTTATCCAAGCACCGACACA‑ACTCCAATACCGGCAAATAAGAAGCCCTCCCTTGGGACAGCTTCAGCTcc                                                                                                                           >  1:587326/1‑250 (MQ=255)
                                                                                                                        gATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTg                                                                                                                      <  1:407943/251‑1 (MQ=255)
                                                                                                                           cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCTGCGTAACGGCGATTTTTTCCTGT‑AACGCGGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑GTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGTTACAGCTGCAGCTCGCGCTGCCt                                                                                                                   <  2:591420/251‑1 (MQ=255)
                                                                                                                           cAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGGGTCGGCCCCCCCGGTCCAGGGTCTTTTCGCCATTGTATG‑TCAGTTTCGGCGTAACGGCGGTTTTTTTCTTT‑TACGCTGCCAGTTGGGGG‑GCGGGGTTTCGCGCCAAAGCGGGTTTTTTATTCAAT‑CCAAGACGGA‑TCAACAATACGGGCAAATAAGATGCCATCCCTTTGTACAGCTTcctctcccgccgcct                                                                                                                   >  2:248612/1‑237 (MQ=255)
                                                                                                                              gAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTaa                                                                                                                <  1:228921/251‑1 (MQ=255)
                                                                                                                                            aaTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAAc                                                                                                  <  1:367201/251‑1 (MQ=255)
                                                                                                                                             aTTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAAcc                                                                                                 <  1:408878/251‑1 (MQ=255)
                                                                                                                                                gtgcggaggggcgctatcAGAATTGCATCGGCGGGGTGTCGCCCATCCAGTTCAATGATGTTTGCGCCATTGAATG‑TTAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTTGTGT‑TGCGTGTTACGCGCAATCGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGTTGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAGCATTTTTATCCAACCGt                                                                                               <  2:453556/233‑1 (MQ=255)
                                                                                                                                                gcgcGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGtt                                                                                              <  1:105678/251‑1 (MQ=255)
                                                                                                                                                   cgTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTAcc                                                                                           <  1:527705/251‑1 (MQ=255)
                                                                                                                                                       gTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGccc                                                                                       <  1:248612/251‑1 (MQ=255)
                                                                                                                                                                             ggcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGc                                                                 <  2:414839/251‑1 (MQ=255)
                                                                                                                                                                              gcggGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGcc                                                                <  1:42970/251‑1 (MQ=255)
                                                                                                                                                                                                           cTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTCATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCAATTCAACCt                                   >  2:259759/1‑251 (MQ=255)
                                                                                                                                                                                                             ttGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATGTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCaa                                      <  1:126453/246‑1 (MQ=255)
                                                                                                                                                                                                              tGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCACCAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGACTTAACATTTTTATCCAACCGTTACCACCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCg                                >  1:490307/1‑251 (MQ=255)
                                                                                                                                                                                                                    aTTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGc                          >  1:453822/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGtgatga              <  1:14210/250‑1 (MQ=255)
                                                                                                                                                                                                                                   gcggtgTAACGGCGATTTTTGCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTGTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGc            <  2:130809/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cggcgTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCa           <  1:4809/251‑1 (MQ=255)
                                                                                                                                                                                                                                          aaCGGCGATTTGTTCCTGT‑AACGCTGCCAGGTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTGCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATGTTTATCCAACCGTTCCCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAg      <  2:556504/250‑1 (MQ=255)
                                                                                                                                                                                                                                             ggCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTCATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCtt  <  1:259759/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCtt  >  1:436031/1‑248 (MQ=255)
                                                                                                                                                                                                                                                gATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCtt  <  2:436031/248‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
CATCACACCCGACCACCGCAATCAGCCCCGGAATGTAGTTGGGCTGCGGCGCGGTTGTACTATTAACAGTTTGAGTTTGTAGTGCATCCATTCGTTAAATCATCCCCTATCCATTTTGGTGATCAGGAAAGCGCGGCGTGAATTGCGCGTAGTGACGCTTGCAGAATTTCATCGGCGGGGTGTCGCCCATCCAGTTCAAGGATCTTTGCGCCATTGAATG‑TCAGTTGCGGCGTAACGGCGATTTTTTCCTGT‑AACGCTGCCAGTTGGTGG‑TCAGGTTTACGCGCAAACGCGG‑TTTGTTCATCAATGCCAAGACGAA‑TCAACAATACGGGCAAATAAGATGCCATCCATTGGTACAGCTTCAGCTCGCGCTGCCTTAACATTTTTATCCAACCGTTACCGCCCGTGGTTTTTGCCAATTGCGGGCCATCAAAGCGAAACCCCGGCACTTCAACCTGCGGGTAGCGGTCGGTGATGAGCAGAAAGCCTT  >  NC_000913/3133122‑3133608

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