Predicted mutation
evidence seq id position mutation annotation gene description
RA W3110S.gb 3,681,551 Δ1 bp coding (1653/1851 nt) ilvD ← dihydroxyacid dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,681,5490G.75.0% 47.1 / 28.0 28coding (1655/1851 nt)ilvDdihydroxyacid dehydratase
Reads supporting (aligned to +/- strand):  ref base G (1/6);  new base . (0/21);  total (1/27)
Fisher's exact test for biased strand distribution p-value = 2.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.09e-01

CGCTTCACGACGCGCCGCCAGTTCGGCATCGCTTACCTGTAACTGAATGCCACGGTTCGGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGC  >  W3110S.gb/3681491‑3681588
                                                          |                                       
cgcTTCACGACGCGCCGCCAGTTCGGCATCGCTTACCTGTAACTGAATGCCACGGTTCGGGa                                      >  1:1366918/1‑62 (MQ=255)
                                ttACCTGTAACTGAATGCCACGGTTCGGGATGTCGATAGCGATCAGGTCACCATCTTCAATc      <  1:811433/62‑1 (MQ=255)
                                ttACCTGTAACTGAATGCCACGGTTCGGGATGTCGATAGCGATCAGGTCACCATCTTCAATc      <  1:683972/62‑1 (MQ=255)
                                ttACCTGTAACTGAATGCCACGGTTCGGGATGTCGATAGCGATCAGGTCACCATCTTCAATc      <  1:1907522/62‑1 (MQ=255)
                                ttACCTGTAACTGAATGCCACGGTTCGGGATGTCGATAGCGATCAGGTCACCATCTTCAATc      <  1:2846189/62‑1 (MQ=255)
                                ttACCTGTAACTGAATGCCACGGTTCGGGATGTCGATAGCGATCAGGTCACCATCTTCAATc      <  1:2436555/62‑1 (MQ=255)
                                 tACCTGTAACTGAATGCCACGGTTCGGGATGTCGATAGCGATCAGGTCACCATCTTCAATc      <  1:2235149/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:2246182/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:982212/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:899624/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:857193/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:462154/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:3065912/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:2989954/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:2972341/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:2844567/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:281022/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:1004862/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:2107710/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:2099083/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:2063316/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:1947881/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:1724492/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:1694448/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:1676636/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:1412850/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:123443/61‑1 (MQ=255)
                                    cTGTAACTGAATGCCACGGTTCGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGc  <  1:1097523/61‑1 (MQ=255)
                                                          |                                       
CGCTTCACGACGCGCCGCCAGTTCGGCATCGCTTACCTGTAACTGAATGCCACGGTTCGGGATGTCGATAGCGATCAGGTCACCATCTTCAATCAGGC  >  W3110S.gb/3681491‑3681588

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: