Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,682,573 G→A 100% L211L (CTG→TTG)  ilvD ← dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,682,5730GA88.0% 56.7 / 2.1 25L211L (CTG→TTG) ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (12/10);  total (15/10)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGCTTACGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCG  >  W3110S.gb/3682515‑3682635
                                                          |                                                              
tgctTACGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTc                                                      >  1:3023593/1‑69 (MQ=255)
  ctTACGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCa                                                     >  1:535816/1‑68 (MQ=255)
  ctTACGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCa                                                     >  1:1568625/1‑68 (MQ=255)
  ctTACGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCa                                                     >  1:524792/1‑68 (MQ=255)
  ctTACGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCa                                                     >  1:2228288/1‑68 (MQ=255)
        gTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAGCGCTTCGGTCAGGCAGtt                                              >  1:1105680/1‑69 (MQ=255)
        gTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGtt                                              <  1:2729505/69‑1 (MQ=255)
                    tGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAg                                   >  1:94372/1‑68 (MQ=255)
                    tGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGc                                  >  1:899083/1‑69 (MQ=255)
                     gCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGc                                  >  1:167004/1‑68 (MQ=255)
                        agcagcGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTa                              <  1:685867/69‑1 (MQ=255)
                             cGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGa                                        <  1:3620812/54‑1 (MQ=255)
                             cGAGCCGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACAt                         <  1:1946661/69‑1 (MQ=255)
                                 ccGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCg                     <  1:2276469/69‑1 (MQ=255)
                                 ccGTTGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCg                     <  1:1846688/69‑1 (MQ=255)
                                    ttGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAg                  <  1:2099150/69‑1 (MQ=255)
                                    ttGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAg                  <  1:3448021/69‑1 (MQ=255)
                                    ttGCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAg                  <  1:1522732/69‑1 (MQ=255)
                                      gCCCGGCTGCGACAGGCCCAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGc                 >  1:348216/1‑68 (MQ=255)
                                      gCCCGGCTGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGc                 >  1:47036/1‑68 (MQ=255)
                                             tGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGa                                        >  1:1777986/1‑38 (MQ=255)
                                             tGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAAc           >  1:2326714/1‑67 (MQ=255)
                                             tGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAAc           >  1:745376/1‑67 (MQ=255)
                                             tGCGACAGGCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAAc           >  1:1773852/1‑67 (MQ=255)
                                                    ggCCCAACGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCg  <  1:1459605/69‑1 (MQ=255)
                                                          |                                                              
TGCTTACGGTCGGCGTGGGTTGCCAGCAGCGAGCCGTTGCCCGGCTGCGACAGGCCCAGCGCTTCGGTCAGGCAGTTCATTGAGTTAGCGGTAAACATCCCGGAGCAGGAACCGCAGGTCG  >  W3110S.gb/3682515‑3682635

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