Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,665,995 T→A 44.1% S254S (TCA→TCT)  guaB ← IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,665,9950TA44.1% 14.3 / 62.6 59S254S (TCA→TCT) guaBIMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (26/7);  new base A (8/18);  total (34/25)
Fisher's exact test for biased strand distribution p-value = 4.39e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

AATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCAA  >  minE/1665938‑1666056
                                                         |                                                             
gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                       >  1:997326/2‑66 (MQ=255)
gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                       >  1:858483/2‑66 (MQ=255)
gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                       >  1:485990/2‑66 (MQ=255)
gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                       >  1:3211438/2‑66 (MQ=255)
gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                       >  1:1356041/2‑66 (MQ=255)
gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                       >  1:2878472/2‑66 (MQ=255)
gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                       >  1:1834286/2‑66 (MQ=255)
gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                                       >  1:1704295/2‑66 (MQ=255)
   ttGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTggg                                                  <  1:2114374/68‑1 (MQ=255)
           ccGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCg                                          >  1:2786754/1‑68 (MQ=255)
                  tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcaggag                                  >  1:2963779/1‑69 (MQ=255)
                  tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                  >  1:1066890/1‑69 (MQ=255)
                  tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                  >  1:612765/1‑69 (MQ=255)
                  tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                  >  1:378671/1‑69 (MQ=255)
                  tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                  >  1:367638/1‑69 (MQ=255)
                  tttAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATcagcag                                  >  1:2591836/1‑69 (MQ=255)
                        aCGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTc                            <  1:3001855/69‑1 (MQ=255)
                         cGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTc                            <  1:1158920/68‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTTCGGCAAc             <  1:2369249/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:2529081/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:3001844/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:1350574/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:3014893/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:3068670/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:3071324/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:3092124/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:2635992/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:2606107/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:1248622/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:1648556/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:2058640/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:1897878/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:941699/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:1816757/69‑1 (MQ=255)
                                       aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:1813324/69‑1 (MQ=255)
                                             cAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             <  1:3076647/63‑1 (MQ=255)
                                                tACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTc    <  1:1276177/69‑1 (MQ=255)
                                                tACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTc    <  1:2299277/69‑1 (MQ=255)
                                                tACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTc    <  1:1871025/69‑1 (MQ=255)
                                                 acacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTc    <  1:282002/68‑1 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGt                             >  1:1081655/1‑42 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc             >  1:1018539/1‑58 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:1220522/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:1654803/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:984878/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:945103/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:1650849/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:19404/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:1115660/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:2001039/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:1427856/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:2469107/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:1541330/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:2668009/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:2700329/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:2734812/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCaa  >  1:1441178/1‑69 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCa   >  1:728669/1‑68 (MQ=255)
                                                  cacCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCa   >  1:2819680/1‑68 (MQ=255)
                                                         |                                                             
AATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAACCAGCGCGTCAA  >  minE/1665938‑1666056

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: ‑