Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,132,787 A→C 100% I204S (ATC→AGC)  wcaA ← putative colanic acid biosynthesis glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,132,7870AC100.0% 61.8 / NA 24I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (11/13);  total (11/13)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCT  >  NC_000913/2132648‑2132922
                                                                                                                                           |                                                                                                                                       
gagcTGGTAGAGGGTATACAGCTTATATTTTTGGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGTAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGctct                                                                                                                                      <  2:677353/140‑1 (MQ=255)
gATCTGGTATAGGGTAAACAGCTGATCCTTGTTGCTGGCGCGGTCGGATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGctct                                                                                                                                      <  2:40261/143‑1 (MQ=255)
  tCTGGTAGAGGGCACACAGTTGTGAGTTTTGGTGGGAGCCGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGCACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                                    <  2:23570/143‑1 (MQ=255)
  gcTGGTAGAGGGTTAACAGCTTATTTGTTTTGTTGGCGCGGTCGTATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACGCTGAGTATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGGTATGCAGGCTCTGc                                                                                                                                    <  2:116056/142‑1 (MQ=255)
            ggTAAACAGCTGATATGTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCCCCGTGCTTGCTATTCAGGCTCTGAGCCGCAtctt                                                                                                                          >  1:33037/1‑143 (MQ=255)
             gTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATCTCACCGTGATTGATATGCAGGCGCTGCGTCGCCtcttc                                                                                                                         >  2:93027/1‑143 (MQ=255)
                      tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTTGGCATTCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGGGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                                >  2:515264/1‑143 (MQ=255)
                             tgtggcaggCGCGGGCGAATTGGGCTTTGTGCTTGCGGGAAGAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                                         <  2:375639/136‑1 (MQ=255)
                               gtGCTGGCGGGGTCGAGTGTGTCTGCGGGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                                       <  1:420097/142‑1 (MQ=255)
                                    ggCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGGCGCCTCTTCTACTTTCCCCGGTTCGCCGGACt                                                                                                  >  1:400319/1‑143 (MQ=255)
                                     gggcgGGCGAAGTTGGCTGTGTGCGTGGGGTAGAAGTGGACGCGCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATGTGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                 <  2:430117/141‑1 (MQ=255)
                                     gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTGTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                 <  1:347283/143‑1 (MQ=255)
                                     gcgcGGGCGGGTGTGGCTTGGGGCGTGCGGTAGAAGTGCCAGTTCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                 <  2:82947/143‑1 (MQ=255)
                                      cgcgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTGCTcc                                                                                                >  2:667716/1‑143 (MQ=255)
                                         ggtccgaTTTGCCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCGCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccacc                                                                                             <  2:210316/137‑1 (MQ=255)
                                           tCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTGCGCCGTACTccaccac                                                                                           >  1:76965/1‑143 (MQ=255)
                                             gAATTTGTCTTTGTGCTGGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacca                                                                                         <  2:451439/143‑1 (MQ=255)
                                                           gCTTGCGGTAAGAGGGGCAGTGGCCGGAGAATTTGTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCGTCTACGTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATg                                                                           <  2:599117/143‑1 (MQ=255)
                                                             ctGCGGTTAGGGTGGGAGTGCGCGGAGAATTGTTTCGGTGACGGGGTGATCTGCATTTCGCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                         <  2:651190/142‑1 (MQ=255)
                                                                     aacgTGGCAGTGCCCGGAGAATGTTTTCGGTGACGGCGTGCTCTGCCTTTCACCGTGATTGATATGCAGGCCCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCt                                                                 <  2:42400/140‑1 (MQ=255)
                                                                                           ttcttCGGTGACGAGGCTACCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATc                                           <  1:664698/140‑1 (MQ=255)
                                                                                                   tGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTGCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATGCTGTGGCGCTTTCAGTCCGGCATCGagcaaac                                   >  1:377853/1‑143 (MQ=255)
                                                                                                                      tCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCATCTACTGTCGGCGATTAGCGGTGCTCCACGCCCATCCGAAGGAAGAAGACATACTGCTGTGCCGCATTCAGTACACTACCGAACATACACCCATGGCAACGCCaaac                >  2:325686/1‑140 (MQ=255)
                                                                                                                              attgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGACCAACGACTCGTTGAACCGCCATGCCCAGCTaa        >  2:620944/1‑143 (MQ=255)
                                                                                                                                  atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGACCAAACACTCTTTGAAACGCCATGCCCAGGTAACGAc    >  1:266650/1‑143 (MQ=255)
                                                                                                                                    atGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGACGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAAAACTCTTCGACACGCCATGCCCAGGCACAGAcca  >  2:529047/1‑142 (MQ=255)
                                                                                                                                           |                                                                                                                                       
GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCT  >  NC_000913/2132648‑2132922

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: