Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,378,032 C→T 100% G95G (GGC→GGT ycjR → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,378,0320CT75.0% 30.2 / 8.9 20G95G (GGC→GGTycjRhypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (5/0);  new base T (0/15);  total (5/15)
Fisher's exact test for biased strand distribution p-value = 6.45e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCT  >  W3110S.gb/1377977‑1378100
                                                       |                                                                    
gCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCt                                                       >  1:2807455/1‑71 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:816721/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:1250525/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:726125/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:706197/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:506946/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:3359981/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:2796529/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:2406385/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:2115492/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:1720714/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:1640797/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:1634456/68‑1 (MQ=255)
              gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:1517937/68‑1 (MQ=255)
               aacgcATTCTTGAAGCGCTGGCAGAAGTTGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:2854038/68‑1 (MQ=255)
                acgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:2919041/68‑1 (MQ=255)
                                                     ggCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCt  >  1:2073260/1‑71 (MQ=255)
                                                     ggCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCt  >  1:3140332/1‑71 (MQ=255)
                                                     ggCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCt  >  1:1775320/1‑71 (MQ=255)
                                                     ggCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCt  >  1:762885/1‑71 (MQ=255)
                                                       |                                                                    
GCTTAAAGCAGATCGAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCT  >  W3110S.gb/1377977‑1378100

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