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,7870AC98.1% 25.7 / ‑6.8 14I204S (ATC→AGC) wcaAputative colanic acid biosynthesis glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base C (6/7);  total (7/7)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAC  >  NC_000913/2132658‑2132920
                                                                                                                                  |                                                                                                                                     
agGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCctc                                                                                                                          <  1:268752/143‑1 (MQ=255)
   gTAATCAGCTGATTTTTTTTGTTGGCGTGGGCGAAGAGGTCGGAGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCtcttc                                                                                                                       <  1:15139/143‑1 (MQ=255)
         aGCTGACATTTTTTGCTGGCGCGGCCGAATATGTCTTTCTGCTTGCGCTATTCTTGGACGCACCCCGCTTACTTTTTCGGTGTCGAGGTGATCTGACCTTCTATCATTTTGTTATGCAGGCTCTGCGTCGCCTCTTCTACttt                                                                                                                 <  2:461175/143‑1 (MQ=255)
           cTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCGCTGCGTCGCACCTTCTACTTTcc                                                                                                               >  2:253868/1‑143 (MQ=255)
                               ggTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCCCTGCTACTCTGCACGGTTCGCCCTACttgcac                                                                                           >  2:250056/1‑138 (MQ=255)
                                             ttGGGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTACATTTGACGATAGTTGTTATGCAGGCTCTGCGTCGCCTGTTCTACTTTCCACGGTTCGGCGTACTGCACCACCATCCGCAGGaa                                                                             >  1:66286/1‑143 (MQ=255)
                                                  cTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGCAGATGt                                                                        >  1:188375/1‑143 (MQ=255)
                                                         tAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTATACTTTCCACAGTTCGCCGTACTCAACCACCATCCGCAGGAAGATGTCGTAATc                                                                 >  1:289807/1‑143 (MQ=255)
                                                                            gggcttttttCGGTGACGTGGTGATCTGCATTTCCCCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                              <  2:156912/139‑1 (MQ=255)
                                                                            agaATTTTTTCGGTGACGATGTGATCTGTATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTGCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCg                                              <  1:60759/143‑1 (MQ=255)
                                                                                          gACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACGATCCGGAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGaacaaaca                                <  1:89289/143‑1 (MQ=255)
                                                                                                  gATCTGCATTTCACCGTGATTAATATGCAGGCTCTGCGTCGCATCTTCTACTTTCCACGGTTCGCCGTCCTTCACCACCATCCGCAGGAAGAACTACTAATCATCTACCGCTCGCAGTCCGGTATCGAACAAacaccccttaa                        >  2:52779/1‑136 (MQ=255)
                                                                                                                 cGTGATTGTTATGCAGGCTCTGCGTCGCCTCTTCTACTGTCCACGGTTCGCCGTTCTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAg          <  2:262311/143‑1 (MQ=255)
                                                                                                                         atatGCAGGCTCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGAACGATGACGTAATCCTGTGCCGCTTTCCGTTCGGTATCGAACAAACCCTCTCTGACACCCCATACCCAGCTaaccac  >  2:341389/1‑139 (MQ=255)
                                                                                                                                  |                                                                                                                                     
AGGGTAAACAGCTGATATTTTTTGCTGGCGCGGTCGAATTTGTCTTTGTGCTTGCGGTAAAAGTGGAAGTACCCGGAGAATTTTTTCGGTGACGAGGTGATCTGCATTTCACCGTGATTGATATGCAGGATCTGCGTCGCCTCTTCTACTTTCCACGGTTCGCCGTACTCCACCACCATCCGCAGGAAGATGTCGTAATCCTGTGCCGCTTTCAGTTCGGTATCGAACAAACACTCTTTGAAACGCCATGCCCAGGTAAAGAC  >  NC_000913/2132658‑2132920

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: