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

GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132647‑2132925
                                                                                                                                            |                                                                                                                                          
gggtttGGTCGAGGGTAAAGAGCGTCTATTTTTGGCGGGCGCGGTCGAATTTGTCTTTGGGTTGGCGGTAAAAGTGGAAGTGTCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTc                                                                                                                                          <  2:379927/138‑1 (MQ=255)
  aTCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATACAGGCTCTg                                                                                                                                        >  2:421444/1‑143 (MQ=255)
     tGGTAGAGGGGAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGt                                                                                                                                     <  2:11731/143‑1 (MQ=255)
          gagGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCct                                                                                                                                >  1:307452/1‑143 (MQ=255)
             ggTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCgctt                                                                                                                             >  1:116771/1‑143 (MQ=255)
              ttgaaCATCTGAGATTTTTGGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                            <  2:729366/140‑1 (MQ=255)
                    cgCTGATAGTGTTTGCTGGCGCGGTCGAATGTGTCTTTGTGCTTGCGGTAAAAGTGGAAGGACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                      <  2:36652/142‑1 (MQ=255)
                     tcgtgtatTTTTTTCTTGCGGTGTCGAATCTGTCTTTCTGCTTGCGGTAGAAGTGTAAGGACTCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGTCTCTGCGTCGCCTCTTCTACTTTc                                                                                                                     <  2:524980/138‑1 (MQ=255)
                                 tGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCg                                                                                                         >  1:357031/1‑143 (MQ=255)
                                   gtGGCGCGGTCGAATTTGTGTGTGTGCTTGCGGGAGAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGGGGGGATCTGCATTTCACCGTGATTGAGATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTa                                                                                                       <  2:276415/142‑1 (MQ=255)
                                      gcgcGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTc                                                                                                    >  2:905103/1‑143 (MQ=255)
                                             ggagtttGCCTTGGTGCGTGCGGTAAGGGTGGAAGTGCCCGGAGAATTTTTTCGGGGACGAGGTGATCTGCATTTCACCGTGATAGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTccaccacc                                                                                             <  2:452760/139‑1 (MQ=255)
                                             cGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTATGAGTAGCCACTTCTTCTTTCCACGGTTCGCCGTACTccaccacc                                                                                             <  1:870229/143‑1 (MQ=255)
                                                     tCTTTGTGCTTGCGGTAAAAGTGGAAGTATCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATGTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAg                                                                                     <  1:756582/143‑1 (MQ=255)
                                                                    tAAACGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATc                                                                      <  1:156446/143‑1 (MQ=255)
                                                                        tgtgtatgTACCCTGAGAATTTTTTTTGTTACGTTGTGTTTTTTATTTCACCTTGATTGATATGCAGGCTCTTCGTCGTCTCTTCTACTTTCCTTGGTTCGCCGTTTTTTATCTTCATCCGCAGGAAGATGTCGTAATCCtgt                                                                  <  2:497955/136‑1 (MQ=255)
                                                                        aGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACGGTGATTGATATGCAGGCTCTGGGTCGCCTCATCTACTTTCGAGGGTTCGCCGGACTCCACCACCATCCGCAGGGAGAGGGCGTGAtgatgt                                                                  >  2:302009/1‑138 (MQ=255)
                                                                                  cccGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCa                                                        >  1:487842/1‑143 (MQ=255)
                                                                                    cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCCGGCAGATATCGTAACCCTCTCCCGcctcccgg                                                      >  2:210550/1‑136 (MQ=255)
                                                                                    cGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATACAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGGACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAgg                                                      >  1:10210/1‑142 (MQ=255)
                                                                                      gagaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTc                                                    <  1:896062/143‑1 (MQ=255)
                                                                                       tgaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   <  1:48450/142‑1 (MQ=255)
                                                                                       agaATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATCGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCTGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                                   <  1:690114/143‑1 (MQ=255)
                                                                                             ttttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCCTAAATCAGTACCGCTAGCAGCTCGGTCTCg                                             >  1:160855/1‑143 (MQ=255)
                                                                                              tttCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGa                                            <  1:905103/143‑1 (MQ=255)
                                                                                                   gTGCCGCGGTGATCTGCTGTTCGCCGTGATTGATCTGCAGGCGCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaa                                       <  2:209529/143‑1 (MQ=255)
                                                                                                   gTGACGAGTTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGATACCGaactga                                       <  2:507560/143‑4 (MQ=255)
                                                                                                              aTCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGCACAACCACTCTTTGag                            >  1:301981/1‑142 (MQ=255)
                                                                                                                   cATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCCTCCGCAGGACGCTGTCGTCCTCCTCTGCCGCTTTCAGTTCGGTATAGAACAAACACTGTGTGACAAGcc                       >  1:656513/1‑143 (MQ=255)
                                                                                                                           cGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCCCCGTACTCCACCACCCTCCGCAGCAAGATGTCATAAAAATATACCACTTTCAATTCGGTATCCAACAAACACGCATTGAAACGCCCTGACCac               >  1:867232/1‑142 (MQ=255)
                                                                                                                             ttattgatATGCAGGCTCTGCTTCGTCTCTTCTACTTTCTACGGTTCGCCGTAGTCCACCACCATCTGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCTAACAAACACTCTTTGAAACGCCATGCCCAGGt             <  2:509119/141‑1 (MQ=255)
                                                                                                                                ttgatATGCAGGCTCTGCTTTTCCTCTTCTACTTTCGCCGGTTCGCCGTACTCCACGACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTaaa          <  2:13904/143‑1 (MQ=255)
                                                                                                                                 tgatATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGTCGTAATCCTGTGCCGCTTTCAGTCCGGTATCGAACCAACACTCTTTGAAGCGCCATGCCCAGGTCAAg         >  1:709799/1‑143 (MQ=255)
                                                                                                                                   atagGCAGGCTCTGCGTCGCCTCCTCTACTTGCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAc       <  1:125228/143‑1 (MQ=255)
                                                                                                                                     atGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCt     <  2:116771/143‑1 (MQ=255)
                                                                                                                                     atGCAGGCTCTGCGTCGCCTCTTCTACTTGCCACGGTTCGCCGTACTCCACCACCGTCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACGAACGCACTTTGAAACGCCATGCCCAGGTAGAGAccc     >  1:38735/1‑142 (MQ=255)
                                                                                                                                        cAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCCGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTAGCGCACAAACACTCTTTGAACCGCCATGCCCAGGTAAAGACCTGGt  >  2:643587/1‑143 (MQ=255)
                                                                                                                                            |                                                                                                                                          
GGATCTGGTAGAGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGACCTGGT  >  NC_000913/2132647‑2132925

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: