Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,701,087 A→G 100% A103A (GCT→GCC rbsC ← D‑ribose transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,701,0870AG78.3% 42.1 / 8.0 23A103A (GCT→GCCrbsCD‑ribose transporter subunit
Reads supporting (aligned to +/- strand):  ref base A (5/0);  new base G (18/0);  total (23/0)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.79e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

CTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATGCAGCAACTGCGCCGGTCAGCGCC  >  W3110S.gb/3701038‑3701156
                                                 |                                                                     
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:2413786/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:61834/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:3314813/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:3218271/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:291153/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:2803528/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:2765680/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:256980/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:2407201/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:2373502/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:209019/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:2010968/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:1965692/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:1824329/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:1558327/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:1381651/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:1264256/1‑71 (MQ=255)
cTACAATCACCCCAGTTACCGCACCAATTGCAGCACCTAACGCGAGAGCGGCAGCGACAGCCACCAGCGCa                                                  >  1:1264507/1‑71 (MQ=255)
                tACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCaa                                        >  1:1101727/1‑65 (MQ=255)
                tACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATg                                      >  1:2728393/1‑67 (MQ=255)
                tACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATg                                      >  1:2825738/1‑67 (MQ=255)
                tACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATg                                      >  1:1310256/1‑67 (MQ=255)
                tACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATg                                      >  1:601516/1‑67 (MQ=255)
                                                cagcagcGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATGCAGCAACTGCGCCGGTCAGCGcc  <  1:2923525/71‑1 (MQ=255)
                                                 |                                                                     
CTACAATCACCCCGGTTACCGCACCAATTGCGGCACCTAACGCGAGAGCAGCAGCGACAGCCACCAGCGCATTGACTTCAATGCCGACGATAGATGCAGCAACTGCGCCGGTCAGCGCC  >  W3110S.gb/3701038‑3701156

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