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% 70.2 / NA 25I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base C (14/11);  total (14/11)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132652‑2132923
                                                                                                                                       |                                                                                                                                        
cggTACAGGGGAAACAGCTGCTATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                   <  1:181062/142‑1 (MQ=255)
  gTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCg                                                                                                                                 >  2:52988/1‑143 (MQ=255)
        ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtctt                                                                                                                           >  1:420911/1‑143 (MQ=255)
                              cTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                     <  2:349780/143‑1 (MQ=255)
                                    cgGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGTAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                               <  1:546290/143‑1 (MQ=255)
                                      gTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACTGGGAGAATTTTGTCGGAAACGAGGTGGTGAAATGGACACCACGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTTCTccacca                                                                                             >  1:417549/1‑143 (MQ=255)
                                                   ttGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGaa                                                                                <  1:14169/143‑1 (MQ=255)
                                                        cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTACCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGt                                                                           <  1:320/143‑1 (MQ=255)
                                                         gggCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                                          <  1:541853/141‑1 (MQ=255)
                                                             ggTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGACCTCGtca                                                                      >  2:377925/1‑141 (MQ=255)
                                                                   aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCtgt                                                                >  1:413070/1‑143 (MQ=255)
                                                                             ctcGGACATTTTTCTCTGTGACGTGGTGATCTTCATTTCACCGTGATTTATATTCAGTCTCTCCGTCGCCTCTTCCACTTTCCACGGTTCGCCTTCCTCCACCACCATCCGCAGGACGATTTCGTAATCCTGTGCCGCTTTCa                                                      <  2:177531/141‑1 (MQ=255)
                                                                                  agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                 >  1:1987/1‑143 (MQ=255)
                                                                                      ttttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGAGGTCGTAATCCTGTGGCGCTTTCAGGTCGGTAt                                             >  1:7694/1‑143 (MQ=255)
                                                                                              gTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                     >  1:402770/1‑143 (MQ=255)
                                                                                                 aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                <  2:115467/49‑1 (MQ=255)
                                                                                                 aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                                >  1:115467/1‑49 (MQ=255)
                                                                                                        gATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGATTTCAGTTCGGTATCGAACAAACACACtttac                           >  2:428215/1‑141 (MQ=255)
                                                                                                         aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACCAACACTCTTTGaa                          >  1:226836/1‑143 (MQ=255)
                                                                                                            tGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCCTTCAGTTCGGTAACGAACAAACACTCTTTCAAACg                       >  2:434251/1‑143 (MQ=255)
                                                                                                               aTTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGGCGTAATCCTGTGCCGCTTTCAGTCCGGGGTCGAACCAACACACTTTGAAGAGCCa                    >  2:422544/1‑143 (MQ=255)
                                                                                                               aTTTCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCa                    <  1:114087/143‑1 (MQ=255)
                                                                                                                      cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTACTCCTGTGACGGTTTCAGTTCGGTAGCGAACAAACACTCATTGAAACGACATGCCAAg             >  2:16538/1‑143 (MQ=255)
                                                                                                                              atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc     <  1:102213/143‑1 (MQ=255)
                                                                                                                                 tGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTg  <  1:52988/143‑1 (MQ=255)
                                                                                                                                       |                                                                                                                                        
TGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTG  >  NC_000913/2132652‑2132923

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