Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,665,980 T→C 42.4% Q259Q (CAA→CAG guaB ← IMP dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,665,9800TC42.4% 13.2 / 43.0 59Q259Q (CAA→CAGguaBIMP dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (7/27);  new base C (8/17);  total (15/44)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.19e-01

GCTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAC  >  minE/1665915‑1666045
                                                                 |                                                                 
gCTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGc                                                                <  1:161073/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1222750/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:2497408/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:628635/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:3305362/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:2588796/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1404789/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:2343865/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1970578/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:2163986/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1762655/69‑1 (MQ=255)
 cTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:3373365/69‑1 (MQ=255)
 cTGTTGCCACGGTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCa                                                               <  1:1888529/69‑1 (MQ=255)
  tGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:19303/69‑1 (MQ=255)
  tGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:3310222/69‑1 (MQ=255)
  tGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:1354960/69‑1 (MQ=255)
  tGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:2011805/69‑1 (MQ=255)
   gTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAg                                                              <  1:355172/68‑1 (MQ=255)
      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:2929705/67‑1 (MQ=255)
      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:2545189/67‑1 (MQ=255)
      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:2510776/67‑1 (MQ=255)
      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:1732689/67‑1 (MQ=255)
      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:620743/67‑1 (MQ=255)
      gCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTa                                                            <  1:840657/67‑1 (MQ=255)
                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:3211438/2‑66 (MQ=255)
                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:2878472/2‑66 (MQ=255)
                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:485990/2‑66 (MQ=255)
                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:858483/2‑66 (MQ=255)
                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:997326/2‑66 (MQ=255)
                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:1356041/2‑66 (MQ=255)
                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:1704295/2‑66 (MQ=255)
                       gaTTTGCAGGTCAGGATATTTAGCGCGAGTTTCACGGATACGCTGCAGAACACCTTCAGAGTGAcc                                            >  1:1834286/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:2591836/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)
                                               aCGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTc                 <  1:3001855/69‑1 (MQ=255)
                                                cGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTc                 <  1:1158920/68‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTTCGGCAAc  <  1:2369249/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3092124/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1897878/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3071324/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1816757/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1813324/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3068670/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3014893/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:3001844/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1648556/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:2058640/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:2635992/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1350574/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:1248622/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:2606107/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:941699/69‑1 (MQ=255)
                                                              aCGCTGCAGAACACCTTCAGAGTGACCGTGAGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAc  <  1:2529081/69‑1 (MQ=255)
                                                                 |                                                                 
GCTGTTGCCACGTTGCCGCCGATAATTTGCAGATCCGGATATTTAGCACGGGTTTCACGGATACGTTGCAGTACACCTTCTGAGTGACCGTGGGAGGAGTCGATCAGCAGAACGTCAACGCCTGCGGCAAC  >  minE/1665915‑1666045

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: