Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,654,260 T→C 51.6% I126I (ATT→ATC purM → phosphoribosylaminoimidazole synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,654,2600TC51.6% 0.4 / 53.0 62I126I (ATT→ATCpurMphosphoribosylaminoimidazole synthetase
Reads supporting (aligned to +/- strand):  ref base T (7/23);  new base C (14/18);  total (21/41)
Fisher's exact test for biased strand distribution p-value = 1.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

CTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCC  >  minE/1654203‑1654319
                                                         |                                                           
cTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGtt                                                   <  1:1961034/68‑1 (MQ=255)
cTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGtt                                                   <  1:2171318/68‑1 (MQ=255)
cTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGtt                                                   <  1:3100389/68‑1 (MQ=255)
     aCGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGc                                             >  1:1098217/1‑69 (MQ=255)
     aCGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGc                                             >  1:2053821/1‑69 (MQ=255)
     aCGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGc                                             >  1:1368415/1‑69 (MQ=255)
     aCGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGc                                             >  1:1338881/1‑69 (MQ=255)
     aCGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGc                                             >  1:1222408/1‑69 (MQ=255)
     aCGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGc                                             >  1:537446/1‑69 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                        <  1:533837/69‑1 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                        <  1:3065052/69‑1 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                        <  1:2720461/69‑1 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                        <  1:758955/69‑1 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                        <  1:956039/69‑1 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                        <  1:1009093/69‑1 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                        <  1:963121/69‑1 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                        <  1:1862861/69‑1 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:1832311/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:1148331/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:1190546/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:642729/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:622617/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:451261/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:3374382/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:3087621/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:1301167/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:2998079/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:2686703/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:2403670/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:1563850/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:2307981/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:184041/69‑2 (MQ=255)
          aCCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:1413728/69‑2 (MQ=255)
           ccGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCa                                        <  1:2551581/68‑1 (MQ=255)
              gAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGtt                                                   <  1:2771220/54‑1 (MQ=255)
                 aaCTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGtt                                 <  1:2719768/69‑1 (MQ=255)
                 aaCTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGtt                                 <  1:120791/69‑1 (MQ=255)
                 aaCTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGtt                                 <  1:2995661/69‑1 (MQ=255)
                    tGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:787441/59‑2 (MQ=255)
                              tACCGCTTCAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcg                                        <  1:1930603/49‑2 (MQ=39)
                               aCCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGgtgggtgga                   <  1:1254876/69‑2 (MQ=255)
                                  gCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGa                 <  1:475694/68‑1 (MQ=255)
                                  gCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGa                 <  1:376433/68‑1 (MQ=255)
                                  gCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGa                 <  1:1448157/68‑1 (MQ=255)
                                  gCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGa                 <  1:116505/68‑1 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGGAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:2927943/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATcgg                                       >  1:2364642/1‑40 (MQ=38)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:1647953/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:782039/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:1246116/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:1824735/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:3324030/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:409887/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:3392501/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:2232599/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:2291997/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:1301203/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:2632510/1‑67 (MQ=255)
                                      cAGCGGTGATCAGCGGCATCGCGGAAGGTTGTCTGCAATCGGGCTGTTCACTGGTGGGTGGCGAAAc              >  1:2332358/1‑67 (MQ=255)
                                          ggTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGa        <  1:3013107/69‑1 (MQ=255)
                                          ggTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGa        <  1:2420300/69‑1 (MQ=255)
                                                cAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGcc  >  1:837375/1‑69 (MQ=255)
                                                         |                                                           
CTATTACGCAACCGGAAAACTGGATGTTGATACCGCTTCAGCGGTGATCAGCGGCATTGCGGAAGGTTGTCTGCAATCAGGCTGTTCACTGGTGGGTGGCGAAACGGCAGAAATGCC  >  minE/1654203‑1654319

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: