Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,978,861 A→C 26.2% T88P (ACC→CCC)  wzyE → putative ECA polysaccharide chain elongation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,978,8610AC26.2% 88.6 / 3.6 62T88P (ACC→CCC) wzyEputative ECA polysaccharide chain elongation protein
Reads supporting (aligned to +/- strand):  ref base A (27/16);  new base C (4/15);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 5.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCTGCTCTGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTAC  >  NC_000913/3978637‑3979088
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
ggCTGCTCTGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTTGTGCAGGCGTTGGTTTCTGCGGGCTGCTTCTACGCGGTTTTCTATTT‑CACCt                                                                                                                                                                                                                                  >  2:493133/1‑228 (MQ=255)
ggCTGCTCTGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGTCGTTTCTTTCTCCGGGCTGCTCCTACCCGGTTTACTATTT‑ccccc                                                                                                                                                                                                                                  >  2:626897/1‑227 (MQ=255)
  ctgctCTGCACGCTGTTTATTGCCACGCTTGCCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTa                                                                                                                                                                                                                                 <  1:287220/227‑1 (MQ=255)
        tGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTTCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAAccc                                                                                                                                                                                                                         >  1:498528/1‑229 (MQ=255)
        tGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GGTTCGCGCCGCCAGAAATCTTGTTGCAGGCGGTGCTTTCTGGGGGCTGCTTCTACGCGGTTTCCTATGT‑CACCTACAAAAc                                                                                                                                                                                                                           >  2:287220/1‑227 (MQ=255)
            cGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCTGCCTGCTTCTACGCGGGTTACTCTGT‑CCCCTACAAAACCCGCCt                                                                                                                                                                                                                     >  1:25506/1‑229 (MQ=255)
             gCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGACGGCGTTGCTTTCTGCGGGCTCCTCCTACGCGGATTACTATGT‑AACCTACAAAAccccccca                                                                                                                                                                                                                    >  2:335736/1‑227 (MQ=255)
                         cACGTTGCCCTGGTTTGAGTTTCGCCGTGTGCTCTTTATCTTCAATGTCTTCTTTTCATTGCTT‑TTTTTGCTCGCCTTTTGCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CCCCTACAAAACCCGCCTACGCAAACGCGt                                                                                                                                                                                                         <  1:687244/228‑1 (MQ=255)
                         cACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CCCCTACAAAACCCGCCTACGCAAACGCGt                                                                                                                                                                                                         >  2:129646/1‑228 (MQ=255)
                         cACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAATACCCGCCTACGCAACCGCGt                                                                                                                                                                                                         >  1:262206/1‑228 (MQ=255)
                           cGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTTTCGCGCCGCCAGCAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTCCTACGCGCTTTACTATTC‑CCCCTACAAAACCCCCCTACGCaacaccgttgg                                                                                                                                                                                                      >  2:285091/1‑220 (MQ=255)
                           cGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTCCGCGGTTTACTATGT‑CACCTCCAAAACCCGCCTACTCAAACGCGTTg                                                                                                                                                                                                       >  1:269938/1‑228 (MQ=255)
                           cGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CCCCTACAAAACCCCCCTACGACAACCCGtgg                                                                                                                                                                                                       >  2:172442/1‑226 (MQ=255)
                            gCTGACCTGGTTTGAGTTTCGCCGTGTTCCCTTTAACTCCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGC‑TTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTg                                                                                                                                                                                                    <  1:85428/229‑1 (MQ=255)
                                                                 ttCAATGTCTTCTTTGCATTGCTG‑TTTTTGCTCATCTTTGTCTTCGGCTCCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTTCGCGGTTTACTATGT‑CCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACcgc                                                                                                                                                                 <  1:234694/228‑1 (MQ=255)
                                                                        tttttttttCATT‑CTGTTTTTTGTTCCTCTTTTTCTTCGGCTCCCCGTGGCCCAGCGTGCTGGTCTTTCGTTTTG‑TTGTGG‑GGGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCGGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGAc                                                                                                                                                         <  2:682372/224‑1 (MQ=255)
                                                                        tttttttgtcgtttttG‑TTTTTGCTTTCTTTTTTTTTCGGCTTCCCTGTGGCAATGGTGCTGGTGTTCGGCTTTG‑ATGTGG‑GTGTCGCGCCGCAAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTATTATGT‑CCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGAc                                                                                                                                                         <  2:243263/214‑1 (MQ=255)
                                                                        tcttctTTTCATGGCTG‑TTTTTGCTCTCGTTTTTCTGCGGCTTCCCGCTGCCCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGAc                                                                                                                                                         <  1:626898/229‑1 (MQ=255)
                                                                             ttttCATTGCTG‑TTTTTGCTCCCCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGTTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATc                                                                                                                                                    <  1:435770/229‑1 (MQ=255)
                                                                                  aTTTCTT‑TTTTTTCTCCCCTTTTGTTGCGGCTTCCCGTTGGCCGGCGTGCTGGTTTTTCGTTTGG‑ATGTTG‑GTGTCGCGCCGCCAGAAATTTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACg                                                                                                                                               <  2:491236/229‑1 (MQ=255)
                                                                                    tGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACAAATCTTActgg                                                                                                                                             >  1:660404/1‑226 (MQ=255)
                                                                                       gtttttttGCTCCCCTTTTTGTGGGGCTTCCCGCTGCCCAGGGGGGTGGTA‑TTCGGTTTGTGTGTGG‑GTGCGGCGCCGTCAGAAATGTTGTTGCGGTCGTTGCTGTCTGCGGGCTGCTTCTGCGCGGTTTACTATGT‑CCCCTAAAAAACCCGCCTACGCAAACGCGGTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTggg                                                                                                                                           <  2:642719/226‑1 (MQ=255)
                                                                                       cgttttttGCTCACTTTTTCCTTCGGCTTCCCGCTGGCCAGCGTGCTTGTATTTCGCTTTG‑ATGTTG‑GTTGCGCGCCGCCAGAAATCTTTTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTggg                                                                                                                                           <  1:194141/226‑1 (MQ=255)
                                                                                                 tctccTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGAt                                                                                                                                 <  1:172442/226‑1 (MQ=255)
                                                                                                      tttttcgcgggCTTCCCGTTGCCCAGGGGGCTGGTATTTCGCTTTG‑ATTTTG‑GTGCCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTa                                                                                                                            <  1:516121/220‑1 (MQ=255)
                                                                                                           cttcGGCTTCCCGCTGCCCGGGTGGCTGGTATTTGGTTTTG‑ATGTTG‑GTGCCCCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATcgcg                                                                                                                       <  1:584746/229‑1 (MQ=255)
                                                                                                           cttcGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATcgc                                                                                                                        <  1:129646/228‑1 (MQ=255)
                                                                                                           ctccggtctccccgtgcccAGCGTGCTGGTTTTTCGCTTTG‑ATGTTG‑GGTTCGCGCCGCAAGAAATCTTGTTGCAGTCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATcgcg                                                                                                                       <  1:632857/212‑1 (MQ=255)
                                                                                                               ggCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTCCGCAAACGCGTTGCTGATGTCCCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGGACCCAATCTTACGGGGGGGATCCTGTTGGGTATcggccggg                                                                                                                   >  2:308652/1‑223 (MQ=255)
                                                                                                                tcTTCCCCGTGCCCGGTGTGCTGTTTTGTTGGTTTG‑ATGTTG‑GGGTCGCGCCGCAAGAACTCGTGTTGGAGGGGTTGTTTCCGGGGGGCTGCTTCTCCGGGGTTTACTATGT‑CCCCTACAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCACGCTGGt                                                                                                                  <  2:22670/228‑1 (MQ=255)
                                                                                                                gttccccgtttcacgggtggCTGTTTTTTCGTTTTG‑A‑GTTGTGGGTCGCGCCGCCAGAATCTTTGTTGCCGGGTTTGTTTTCGGGGGGTGGCTTCTGCGCGTTTTATTATGT‑CCCCTTCAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGt                                                                                                                  <  2:230278/210‑1 (MQ=255)
                                                                                                                 cTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGTGGCTCCTGATGGGTATCGCGCTGGt                                                                                                                  >  1:462123/1‑228 (MQ=255)
                                                                                                                    cccGCTGACCCGCGTGCGGTGATTTCGCTGTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGGGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTAT‑T‑CCCCTCCAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCg                                                                                                             <  1:28731/229‑1 (MQ=255)
                                                                                                                     ccGCTGACCTGCGTGTTGGTTTTTCGTTTTG‑ATGGTG‑GTGGCGCGCCGCAAGAAATCTTGTTGGCGGCGTTGGTTCTGGGGGGCTTCTTCGACGCGGTTTACTATGG‑CCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTATGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCg                                                                                                             <  2:225014/229‑1 (MQ=255)
                                                                                                                              agttgttgggTTTTTCGTTTTG‑TTGTGG‑GTTGCGCGCCGCCAAAAATCTTGTGGGAGGCTTTGTTTTCTGCGGGCTTCTGTTACGGGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtct                                                                                                    <  2:498528/221‑1 (MQ=255)
                                                                                                                                cGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCCtcttc                                                                                                  >  1:20703/1‑229 (MQ=255)
                                                                                                                                  tGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCTCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTCCTGCTGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGAc                                                                                                                                                         >  1:645534/1‑172 (MQ=255)
                                                                                                                                     tGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCTTCCGCTGTTTACCATGAACCGCCTGGGGACCACTCTTACGTGGGGGATCCTGATGGGGATCGCCCTGGTAAGCGTCGGCACCttcttc                                                                                               >  2:560048/1‑227 (MQ=255)
                                                                                                                                                 ttG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCTTTGCTGATGTACCGCGCCGTCCCCTGTTTACCATGAACCGCGTGGATACCCATCTTCCTTGGGTGATCCTGATGGGTAACAGCCTGTACACCGCCGGCACCTTCTTCATGCCCAACGGCt                                                                                  >  2:104673/1‑228 (MQ=255)
                                                                                                                                                 ttG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCt                                                                                  >  1:542092/1‑228 (MQ=255)
                                                                                                                                                      ttttG‑GGGTCGCCGCGACAAAAATCTTGTTGCAGGCGTTGCTTTGTGCGGGCTGTTTCTACGCGGTTTACTAGGT‑CACCTACAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTg                                                                             <  1:672087/227‑1 (MQ=255)
                                                                                                                                                          g‑gTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTg                                                                          >  2:526304/1‑228 (MQ=255)
                                                                                                                                                             tgtCGCGCCGCCAGAAGTTTTGTTGCGGTCGTTGCTTTCTGCGGGCGCGTTGTACGCGTTTTACTACTC‑CCCCTACAAATCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTc                                                                       <  2:479309/229‑1 (MQ=255)
                                                                                                                                                                   gccgccAGAATTTTTTTGGCAAGCTTTGCTTTGCCCGGGCTGCTTCTACGCGGTTTAC‑GTGTCCCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCtt                                                                 <  2:269938/229‑1 (MQ=255)
                                                                                                                                                                   gccgccAGAAATTTTGGTGCCGGTGTTGCTTTCCGCGGGCTGCTTCCACGGGGTTTACTATGT‑CCCTTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCtt                                                                 <  1:216723/229‑1 (MQ=255)
                                                                                                                                                                     cgcaAGAATTTTTGTTGCAGGCGGTGCTTTCGGCGGGTTGCTTTTACGCGGTTTATTATGT‑CCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGACTTaa                                                               <  1:335736/229‑1 (MQ=255)
                                                                                                                                                                          taaaTCTGTTTGCAGGCGTTGTTTTCTTGGGTCTTCTTCTCCGCGGTTTGATATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGGGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCt                                                          <  1:505574/228‑1 (MQ=255)
                                                                                                                                                                          agatttttgttgCAGGCGTGGCTTTCTGGGGGTCGCTTCTACGCGGTTTACTATTC‑CACCTCCAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCt                                                          <  1:652222/223‑1 (MQ=255)
                                                                                                                                                                           aaaTTTTGTTGGAGGCGTTGCTTTCTGCGGGCTGTTTCTACGGGGTTTCCTAGTC‑CCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTa                                                         <  1:676496/229‑1 (MQ=255)
                                                                                                                                                                               ttttttggcggCGTTTTTTTTTGCGGGCTTTTTCCACGGGGTTTACTATGT‑CCCCTACAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGt                                                     <  2:20703/220‑1 (MQ=255)
                                                                                                                                                                                  gttgCAGGCGTTGCTTTCTGGGGGCTGATTCTACGCGGTTTACTATGT‑CACCTCCAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAg                                                  <  1:287034/229‑1 (MQ=255)
                                                                                                                                                                                        ggCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGCGCGATCCTGATGGGAATCGCGCTGGTAAGCGTCGGCACCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTAAAGTCAGATCttt                                            >  1:479396/1‑229 (MQ=255)
                                                                                                                                                                                             tGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGGGATCCTGATGGGTATCGCGCTGGTAAGCGTCGCCATCTTCTTCATGCACACCGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGCCAGATCTTTTCCAg                                       >  2:227946/1‑229 (MQ=255)
                                                                                                                                                                                                              ttCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtg                       >  1:17146/1‑228 (MQ=255)
                                                                                                                                                                                                                       gTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGCGACCACTCTTACGTGGGTGATCCTCATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACACTCCGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAAc             >  1:508358/1‑229 (MQ=255)
                                                                                                                                                                                                                        tttACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGATGATCCTGATGGGTATCGCGCGGGTAAGCGTCGGCATCTTCTTCATGCCCAACGGCTTTTTGCTGTTCCGGCTTAACCCCTACAGTACGATCTTTTCCGGTGAACACTCCGGGGGGGCGGTaaaaac            >  2:87327/1‑227 (MQ=255)
                                                                                                                                                                                                                          tACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGCGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCACCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCCAATCTTTTCCCGTGAGTTCTCCGGGCTGGCATTacaaacat          >  2:224354/1‑225 (MQ=255)
                                                                                                                                                                                                                                   cACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCCCAACGGCTTTTTGCTGTTCCGGCTTAACTCCTCCAGTCAGACCTTTTCCATTGAAGTCTCCGGAGTGGCGTTAAAAAGCTTCTTTTa   >  2:348482/1‑228 (MQ=255)
                                                                                                                                                                                                                                   cACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTCCGGTCAGATCTTTTCCAGTGAAGTTTCCGGCGTGCCGTTAAAACGCTTCTTTTa   >  1:307219/1‑228 (MQ=255)
                                                                                                                                                                                                                                   cACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGAACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACCGTCAGATCTTTTCCAGTGAATGCTCCGGCTTGGCGTTAAAAAGCTTCTTTTa   >  1:297743/1‑228 (MQ=255)
                                                                                                                                                                                                                                   cACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGAATCTTCTTCATGCACAACGGCTTTTTGCGGTTACGGCTTAACTCCTACAGTCAGATATTTCCCAGTGACGTATCTGGCGTGGCGTTAAAACCCTTCTTTTat  >  2:267501/1‑228 (MQ=255)
                                                                                                                                                                                                                                   cACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGCGATCCTGACGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACATTCAGACATTTTCCAGTGAAGTCTCGTGCGCGGCGTTAAAACGCTTCTTTTat  >  2:133639/1‑228 (MQ=255)
                                                                                                                                                                                                                                   cACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCCCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGTCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACCCCTACAGTCCGATCCTTTCCTGTGAAGTTCCCGGCGGGGCTATAAAACGCGTCtttttg  >  2:581126/1‑227 (MQ=255)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
GGCTGCTCTGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTG‑TTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTG‑ATGTTG‑GTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGT‑CACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTAC  >  NC_000913/3978637‑3979088

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