Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,111,828:1 +C 100% coding (917/1536 nt) mdoG → glucan biosynthesis protein, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,111,8281.C75.0% 14.9 / 1.6 8coding (917/1536 nt)mdoGglucan biosynthesis protein, periplasmic
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (6/0);  total (6/2)
Fisher's exact test for biased strand distribution p-value = 3.57e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic indel expands or contracts a homopolymer stretch.

AATGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCGCAAGGCTTTGGTCTGTTGCAGCGCGGTCGTGATTTCTCCCGCTTTGAAGATCTCGATG  >  W3110S.gb/1111758‑1111890
                                                                       |                                                              
aaTGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGaaaa                                                                 >  1:1156066/1‑71 (MQ=255)
aaTGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGaaaa                                                                 >  1:677358/1‑71 (MQ=255)
aaTGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGaaaa                                                                 >  1:660971/1‑71 (MQ=255)
aaTGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGaaaa                                                                 >  1:1614999/1‑71 (MQ=255)
 aTGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGaaaa                                                                 <  1:118011/70‑1 (MQ=255)
 aTGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGaaaa                                                                 <  1:1587185/70‑1 (MQ=255)
   ggCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAAccc                                                             >  1:1512320/1‑71 (MQ=255)
   ggCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAAccc                                                             >  1:393033/1‑71 (MQ=255)
   ggCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAAcc                                                              >  1:974721/1‑70 (MQ=255)
                            aTAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCCGCAAGGCTTTGGTCTGTTGCAgcg                                     >  1:1478/1‑71 (MQ=255)
                            aTAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCCGCAAGGCTTTGGTCTGTTGCAgcg                                     >  1:1454809/1‑71 (MQ=255)
                            aTAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCCGCAAGGCTTTGGTCTGTTGCAgcg                                     >  1:1379599/1‑71 (MQ=255)
                            aTAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCCGCAAGGCTTTGGTCTGTTGCAgcg                                     >  1:1251727/1‑71 (MQ=255)
                            aTAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCCGCAAGGCTTTGGTCTGTTGCAgcg                                     >  1:734700/1‑71 (MQ=255)
                            aTAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCCGCAAGGCTTTGGTCTGTTGCAgc                                      >  1:748839/1‑70 (MQ=255)
                                                              cATGGAAAACCCGCAAGGCTTTGGTCTGTTGCAGCGCGGTCGTGATTTCTCCCGCTTTGAAGATCTCgatg  <  1:438061/71‑1 (MQ=255)
                                                              cATGGAAAACCCGCAAGGCTTTGGTCTGTTGCAGCGCGGTCGTGATTTCTCCCGCTTTGAAGATCTCgatg  <  1:1246765/71‑1 (MQ=255)
                                                                       |                                                              
AATGGCGAATGGATCTGGCGTCCGTTGAATAACCCGAAACATTTAGCGGTCAGCAGCTTCTCCATGGAAAACCCGCAAGGCTTTGGTCTGTTGCAGCGCGGTCGTGATTTCTCCCGCTTTGAAGATCTCGATG  >  W3110S.gb/1111758‑1111890

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