Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,378,050 G→A 100% A101A (GCG→GCA ycjR → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,378,0500GA79.0% 33.5 / 5.3 19A101A (GCG→GCAycjRhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (4/0);  new base A (0/15);  total (4/15)
Fisher's exact test for biased strand distribution p-value = 2.58e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCT  >  W3110S.gb/1377991‑1378100
                                                           |                                                  
gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:2796529/68‑1 (MQ=255)
gaacgcATTCTTGAAGCGCTGGCAGAAGTGGGCGGTAAAGGTATCGTCGTTCCGGCTGCATGGGGCATGtt                                         <  1:816721/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:1250525/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)
                                                           |                                                  
GAGCGCATTCTCGAAGCGCTGGCAGAAGTGGGCGGTAAAGGCATCGTCGTTCCGGCTGCGTGGGGCATGTTTACCTTCCGCTTACCGCCGATGACCTCGCCGCGTAGCCT  >  W3110S.gb/1377991‑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: