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

GGCTGCTCTGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTAC  >  NC_000913/3978637‑3979088
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
ggCTGCTCTGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTTGTGCAGGCGTTGGTTTCTGCGGGCTGCTTCTACGCGGTTTTCTATTTCACCt                                                                                                                                                                                                                                  >  2:493133/1‑228 (MQ=255)
ggCTGCTCTGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGTCGTTTCTTTCTCCGGGCTGCTCCTACCCGGTTTACTATTTccccc                                                                                                                                                                                                                                  >  2:626897/1‑227 (MQ=255)
  ctgctCTGCACGCTGTTTATTGCCACGCTTGCCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTa                                                                                                                                                                                                                                 <  1:287220/227‑1 (MQ=255)
        tGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTTCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAAccc                                                                                                                                                                                                                         >  1:498528/1‑229 (MQ=255)
        tGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGGTTCGCGCCGCCAGAAATCTTGTTGCAGGCGGTGCTTTCTGGGGGCTGCTTCTACGCGGTTTCCTATGTCACCTACAAAAc                                                                                                                                                                                                                           >  2:287220/1‑227 (MQ=255)
            cGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCTGCCTGCTTCTACGCGGGTTACTCTGTCCCCTACAAAACCCGCCt                                                                                                                                                                                                                     >  1:25506/1‑229 (MQ=255)
             gCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGACGGCGTTGCTTTCTGCGGGCTCCTCCTACGCGGATTACTATGTAACCTACAAAAccccccca                                                                                                                                                                                                                    >  2:335736/1‑227 (MQ=255)
                         cACGTTGCCCTGGTTTGAGTTTCGCCGTGTGCTCTTTATCTTCAATGTCTTCTTTTCATTGCTTTTTTTGCTCGCCTTTTGCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCCCCTACAAAACCCGCCTACGCAAACGCGt                                                                                                                                                                                                         <  1:687244/228‑1 (MQ=255)
                         cACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCCCCTACAAAACCCGCCTACGCAAACGCGt                                                                                                                                                                                                         >  2:129646/1‑228 (MQ=255)
                         cACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAATACCCGCCTACGCAACCGCGt                                                                                                                                                                                                         >  1:262206/1‑228 (MQ=255)
                           cGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTTTCGCGCCGCCAGCAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTCCTACGCGCTTTACTATTCCCCCTACAAAACCCCCCTACGCaacaccgttgg                                                                                                                                                                                                      >  2:285091/1‑220 (MQ=255)
                           cGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTCCGCGGTTTACTATGTCACCTCCAAAACCCGCCTACTCAAACGCGTTg                                                                                                                                                                                                       >  1:269938/1‑228 (MQ=255)
                           cGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCCCCTACAAAACCCCCCTACGACAACCCGtgg                                                                                                                                                                                                       >  2:172442/1‑226 (MQ=255)
                            gCTGACCTGGTTTGAGTTTCGCCGTGTTCCCTTTAACTCCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTg                                                                                                                                                                                                    <  1:85428/229‑1 (MQ=255)
                                                                 ttCAATGTCTTCTTTGCATTGCTGTTTTTGCTCATCTTTGTCTTCGGCTCCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTTCGCGGTTTACTATGTCCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACcgc                                                                                                                                                                 <  1:234694/228‑1 (MQ=255)
                                                                        tttttttttCATTCTGTTTTTTGTTCCTCTTTTTCTTCGGCTCCCCGTGGCCCAGCGTGCTGGTCTTTCGTTTTGTTGTGGGGGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCGGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGAc                                                                                                                                                         <  2:682372/224‑1 (MQ=255)
                                                                        tttttttgtcgtttttGTTTTTGCTTTCTTTTTTTTTCGGCTTCCCTGTGGCAATGGTGCTGGTGTTCGGCTTTGATGTGGGTGTCGCGCCGCAAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTATTATGTCCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGAc                                                                                                                                                         <  2:243263/214‑1 (MQ=255)
                                                                        tcttctTTTCATGGCTGTTTTTGCTCTCGTTTTTCTGCGGCTTCCCGCTGCCCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGAc                                                                                                                                                         <  1:626898/229‑1 (MQ=255)
                                                                             ttttCATTGCTGTTTTTGCTCCCCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGTTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATc                                                                                                                                                    <  1:435770/229‑1 (MQ=255)
                                                                                  aTTTCTTTTTTTTCTCCCCTTTTGTTGCGGCTTCCCGTTGGCCGGCGTGCTGGTTTTTCGTTTGGATGTTGGTGTCGCGCCGCCAGAAATTTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACg                                                                                                                                               <  2:491236/229‑1 (MQ=255)
                                                                                    tGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACAAATCTTActgg                                                                                                                                             >  1:660404/1‑226 (MQ=255)
                                                                                       gtttttttGCTCCCCTTTTTGTGGGGCTTCCCGCTGCCCAGGGGGGTGGTATTCGGTTTGTGTGTGGGTGCGGCGCCGTCAGAAATGTTGTTGCGGTCGTTGCTGTCTGCGGGCTGCTTCTGCGCGGTTTACTATGTCCCCTAAAAAACCCGCCTACGCAAACGCGGTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTggg                                                                                                                                           <  2:642719/226‑1 (MQ=255)
                                                                                       cgttttttGCTCACTTTTTCCTTCGGCTTCCCGCTGGCCAGCGTGCTTGTATTTCGCTTTGATGTTGGTTGCGCGCCGCCAGAAATCTTTTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTggg                                                                                                                                           <  1:194141/226‑1 (MQ=255)
                                                                                                 tctccTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGAt                                                                                                                                 <  1:172442/226‑1 (MQ=255)
                                                                                                      tttttcgcgggCTTCCCGTTGCCCAGGGGGCTGGTATTTCGCTTTGATTTTGGTGCCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTa                                                                                                                            <  1:516121/220‑1 (MQ=255)
                                                                                                           cttcGGCTTCCCGCTGCCCGGGTGGCTGGTATTTGGTTTTGATGTTGGTGCCCCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATcgcg                                                                                                                       <  1:584746/229‑1 (MQ=255)
                                                                                                           cttcGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATcgc                                                                                                                        <  1:129646/228‑1 (MQ=255)
                                                                                                           ctccggtctccccgtgcccAGCGTGCTGGTTTTTCGCTTTGATGTTGGGTTCGCGCCGCAAGAAATCTTGTTGCAGTCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATcgcg                                                                                                                       <  1:632857/212‑1 (MQ=255)
                                                                                                               ggCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTCCGCAAACGCGTTGCTGATGTCCCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGGACCCAATCTTACGGGGGGGATCCTGTTGGGTATcggccggg                                                                                                                   >  2:308652/1‑223 (MQ=255)
                                                                                                                tcTTCCCCGTGCCCGGTGTGCTGTTTTGTTGGTTTGATGTTGGGGTCGCGCCGCAAGAACTCGTGTTGGAGGGGTTGTTTCCGGGGGGCTGCTTCTCCGGGGTTTACTATGTCCCCTACAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCACGCTGGt                                                                                                                  <  2:22670/228‑1 (MQ=255)
                                                                                                                gttccccgtttcacgggtggCTGTTTTTTCGTTTTGAGTTGTGGGTCGCGCCGCCAGAATCTTTGTTGCCGGGTTTGTTTTCGGGGGGTGGCTTCTGCGCGTTTTATTATGTCCCCTTCAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGt                                                                                                                  <  2:230278/210‑1 (MQ=255)
                                                                                                                 cTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGTGGCTCCTGATGGGTATCGCGCTGGt                                                                                                                  >  1:462123/1‑228 (MQ=255)
                                                                                                                    cccGCTGACCCGCGTGCGGTGATTTCGCTGTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGGGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATTCCCCTCCAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCg                                                                                                             <  1:28731/229‑1 (MQ=255)
                                                                                                                     ccGCTGACCTGCGTGTTGGTTTTTCGTTTTGATGGTGGTGGCGCGCCGCAAGAAATCTTGTTGGCGGCGTTGGTTCTGGGGGGCTTCTTCGACGCGGTTTACTATGGCCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTATGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCg                                                                                                             <  2:225014/229‑1 (MQ=255)
                                                                                                                              agttgttgggTTTTTCGTTTTGTTGTGGGTTGCGCGCCGCCAAAAATCTTGTGGGAGGCTTTGTTTTCTGCGGGCTTCTGTTACGGGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCAtct                                                                                                    <  2:498528/221‑1 (MQ=255)
                                                                                                                                cGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCCtcttc                                                                                                  >  1:20703/1‑229 (MQ=255)
                                                                                                                                  tGCTGGTATTTCGCTTTGATGTTGGTGTCTCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTCCTGCTGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGAc                                                                                                                                                         >  1:645534/1‑172 (MQ=255)
                                                                                                                                     tGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCTTCCGCTGTTTACCATGAACCGCCTGGGGACCACTCTTACGTGGGGGATCCTGATGGGGATCGCCCTGGTAAGCGTCGGCACCttcttc                                                                                               >  2:560048/1‑227 (MQ=255)
                                                                                                                                                 ttGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCTTTGCTGATGTACCGCGCCGTCCCCTGTTTACCATGAACCGCGTGGATACCCATCTTCCTTGGGTGATCCTGATGGGTAACAGCCTGTACACCGCCGGCACCTTCTTCATGCCCAACGGCt                                                                                  >  2:104673/1‑228 (MQ=255)
                                                                                                                                                 ttGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCt                                                                                  >  1:542092/1‑228 (MQ=255)
                                                                                                                                                      ttttGGGGTCGCCGCGACAAAAATCTTGTTGCAGGCGTTGCTTTGTGCGGGCTGTTTCTACGCGGTTTACTAGGTCACCTACAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTg                                                                             <  1:672087/227‑1 (MQ=255)
                                                                                                                                                          ggTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTg                                                                          >  2:526304/1‑228 (MQ=255)
                                                                                                                                                             tgtCGCGCCGCCAGAAGTTTTGTTGCGGTCGTTGCTTTCTGCGGGCGCGTTGTACGCGTTTTACTACTCCCCCTACAAATCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTc                                                                       <  2:479309/229‑1 (MQ=255)
                                                                                                                                                                   gccgccAGAATTTTTTTGGCAAGCTTTGCTTTGCCCGGGCTGCTTCTACGCGGTTTACGTGTCCCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCtt                                                                 <  2:269938/229‑1 (MQ=255)
                                                                                                                                                                   gccgccAGAAATTTTGGTGCCGGTGTTGCTTTCCGCGGGCTGCTTCCACGGGGTTTACTATGTCCCTTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCtt                                                                 <  1:216723/229‑1 (MQ=255)
                                                                                                                                                                     cgcaAGAATTTTTGTTGCAGGCGGTGCTTTCGGCGGGTTGCTTTTACGCGGTTTATTATGTCCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGACTTaa                                                               <  1:335736/229‑1 (MQ=255)
                                                                                                                                                                          taaaTCTGTTTGCAGGCGTTGTTTTCTTGGGTCTTCTTCTCCGCGGTTTGATATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGGGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCt                                                          <  1:505574/228‑1 (MQ=255)
                                                                                                                                                                          agatttttgttgCAGGCGTGGCTTTCTGGGGGTCGCTTCTACGCGGTTTACTATTCCACCTCCAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCt                                                          <  1:652222/223‑1 (MQ=255)
                                                                                                                                                                           aaaTTTTGTTGGAGGCGTTGCTTTCTGCGGGCTGTTTCTACGGGGTTTCCTAGTCCCCCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTa                                                         <  1:676496/229‑1 (MQ=255)
                                                                                                                                                                               ttttttggcggCGTTTTTTTTTGCGGGCTTTTTCCACGGGGTTTACTATGTCCCCTACAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGt                                                     <  2:20703/220‑1 (MQ=255)
                                                                                                                                                                                  gttgCAGGCGTTGCTTTCTGGGGGCTGATTCTACGCGGTTTACTATGTCACCTCCAAACCCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAg                                                  <  1:287034/229‑1 (MQ=255)
                                                                                                                                                                                        ggCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGCGCGATCCTGATGGGAATCGCGCTGGTAAGCGTCGGCACCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTAAAGTCAGATCttt                                            >  1:479396/1‑229 (MQ=255)
                                                                                                                                                                                             tGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGGGATCCTGATGGGTATCGCGCTGGTAAGCGTCGCCATCTTCTTCATGCACACCGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGCCAGATCTTTTCCAg                                       >  2:227946/1‑229 (MQ=255)
                                                                                                                                                                                                              ttCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCggcgtg                       >  1:17146/1‑228 (MQ=255)
                                                                                                                                                                                                                       gTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGCGACCACTCTTACGTGGGTGATCCTCATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACACTCCGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAAc             >  1:508358/1‑229 (MQ=255)
                                                                                                                                                                                                                        tttACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGATGATCCTGATGGGTATCGCGCGGGTAAGCGTCGGCATCTTCTTCATGCCCAACGGCTTTTTGCTGTTCCGGCTTAACCCCTACAGTACGATCTTTTCCGGTGAACACTCCGGGGGGGCGGTaaaaac            >  2:87327/1‑227 (MQ=255)
                                                                                                                                                                                                                          tACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGCGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCACCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCCAATCTTTTCCCGTGAGTTCTCCGGGCTGGCATTacaaacat          >  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)
                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
GGCTGCTCTGCACGCTGTTTATTGCCACGCTGACCTGGTTTGAGTTTCGCCGTGTGCGCTTTAACTTCAATGTCTTCTTTTCATTGCTGTTTTTGCTCACCTTTTTCTTCGGCTTCCCGCTGACCAGCGTGCTGGTATTTCGCTTTGATGTTGGTGTCGCGCCGCCAGAAATCTTGTTGCAGGCGTTGCTTTCTGCGGGCTGCTTCTACGCGGTTTACTATGTCACCTACAAAACCCGCCTACGCAAACGCGTTGCTGATGTACCGCGCCGTCCGCTGTTTACCATGAACCGCGTGGAGACCAATCTTACGTGGGTGATCCTGATGGGTATCGCGCTGGTAAGCGTCGGCATCTTCTTCATGCACAACGGCTTTTTGCTGTTCCGGCTTAACTCCTACAGTCAGATCTTTTCCAGTGAAGTCTCCGGCGTGGCGTTAAAACGCTTCTTTTAC  >  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: