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

GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCC  >  NC_000913/2132648‑2132910
                                                                                                                                           |                                                                                                                           
gATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGctct                                                                                                                          >  1:429715/1‑143 (MQ=255)
  tCTGGTAGAGGGTAAACAGCTGCTATCATTAGCGGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                        <  1:18607/143‑1 (MQ=255)
  tCTGGGAGAGCGTAAACAGCCTATGTTGTGTGCTGGCGCGGTCGATTTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGTAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGc                                                                                                                        <  2:321461/143‑1 (MQ=255)
       gagagGGGAAACAGCTGATCACCTACGCTGGCGCGGGCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGc                                                                                                                   <  2:91778/142‑1 (MQ=255)
                    gCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTc                                                                                                      >  1:90733/1‑143 (MQ=255)
                      tGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCa                                                                                                    >  1:413248/1‑143 (MQ=255)
                       gaaagtttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCAc                                                                                                   <  1:492792/138‑1 (MQ=255)
                             ttttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAACTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCg                                                                                             <  1:138252/143‑1 (MQ=255)
                               ttGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGcc                                                                                           >  2:24405/1‑143 (MQ=255)
                                        ccgTCGAATTTGTCTTTGGGCCTGCGGTAAAAGTGGAAGTACGCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccac                                                                                  <  2:21862/141‑1 (MQ=255)
                                              agtttGTCTGTGGGCTTGGGGTAAGAGTGGAAGTACCCTGAGAATTTTTTCGGTGACGAGGTGAGCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCAt                                                                            <  2:501772/141‑1 (MQ=255)
                                                   gTCTGTGTGCTGGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCa                                                                       <  1:80429/143‑1 (MQ=255)
                                                        tgtgCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAg                                                                  <  1:94519/143‑1 (MQ=255)
                                                             ttGCGTTAAAAGTGGTTGTACCCTTAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCGCCGTGATTGATATCTTGGCTTTGTGTCGCCTCTTCTACTTTCCTCGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTc                                                             <  2:281114/143‑1 (MQ=255)
                                                                cGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGCTGTCGTa                                                          >  1:91541/1‑143 (MQ=255)
                                                                cGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTa                                                          >  1:341751/1‑143 (MQ=255)
                                                                 ggTAAAAGTGGAAGTGCCCGGAGAATGTGGTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTaa                                                         <  1:578894/143‑1 (MQ=255)
                                                                    aaaaGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCCTACTCCACCCCCATCCGAAGGAACATGTCGTAATcc                                                      >  2:159557/1‑143 (MQ=255)
                                                                                aCCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCCCCCTCCGCAGGAAGATGTCGTACACCTGTGCCCCttcc                                          >  2:218521/1‑141 (MQ=255)
                                                                                  ccGGAGAATTTTTTCGGTGACGAGAGGACCTGCATTGCACCGTGATTGATATGCCGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTAATCCACCACCATCCGCAGGAAGCTCTCGCAATCCTGTGCCGCTTTCAg                                        >  2:352209/1‑143 (MQ=255)
                                                                                        aaTTTTTTCTGTGACCAGGTAATCTCCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGTAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGt                                  <  1:306347/143‑1 (MQ=255)
                                                                                                cGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaaca                          <  1:251349/143‑1 (MQ=255)
                                                                                                    ggcgtggtggtCTGCTTTTCGCCGGTATTGGTGCGCCGGCTGTGCGTCGCCTCTTCTAGTTTCCAGGGTTCGCCGTACTCCACCACCTTCCGCAGGACGAGGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                      <  2:463065/138‑1 (MQ=255)
                                                                                                     aCGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAacacacac                     >  2:375919/1‑143 (MQ=255)
                                                                                                      cGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaacaat                    >  1:416198/1‑141 (MQ=255)
                                                                                                                        aCCGTGATTGACATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGccc  <  2:329868/143‑1 (MQ=255)
                                                                                                                                           |                                                                                                                           
GATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCC  >  NC_000913/2132648‑2132910

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