Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,224,954 C→T 49.4% A201A (GCC→GCT prmA → methylase for 50S ribosomal subunit protein L11

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,224,9540CT49.4% ‑2.4 / 85.9 93A201A (GCC→GCTprmAmethylase for 50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base C (28/19);  new base T (16/30);  total (44/49)
Fisher's exact test for biased strand distribution p-value = 2.25e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAATG  >  minE/2224887‑2225017
                                                                   |                                                               
tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                                >  1:3157534/1‑69 (MQ=255)
tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                                >  1:2322007/1‑69 (MQ=255)
tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                                >  1:2772431/1‑69 (MQ=255)
tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                                >  1:705883/1‑69 (MQ=255)
tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                                >  1:1203342/1‑69 (MQ=255)
    tGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGaa                                                              <  1:592433/67‑1 (MQ=255)
    tGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGaa                                                              <  1:2572042/67‑1 (MQ=255)
     gcagcaAAAGCCATTGTTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                           <  1:3259328/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                           <  1:1006006/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                           <  1:2139701/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                           <  1:3404743/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                           <  1:2718619/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                           <  1:3096298/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                           <  1:1217643/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                           <  1:2515695/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                           <  1:2654362/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                           <  1:887969/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                           <  1:3312036/69‑1 (MQ=255)
     gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                           <  1:564507/69‑1 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:2812545/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:2986724/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:1029492/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:2668988/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:2628509/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:2594231/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:419927/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:565986/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:737843/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:1611452/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:1452275/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:2141375/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:1183518/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:2117780/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:1360518/1‑69 (MQ=255)
       agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTAAGGCCAGCCGCGACAACGCTGAACGTaa                                                         >  1:1428918/1‑69 (MQ=255)
           aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGc                                                     >  1:3360304/1‑69 (MQ=255)
           aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGc                                                     >  1:1981995/1‑69 (MQ=255)
           aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGc                                                     >  1:1982804/1‑69 (MQ=255)
           aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGc                                                     >  1:3101952/1‑69 (MQ=255)
                    ggTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGAc                                            >  1:760232/1‑69 (MQ=255)
                             atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATgcg                                                     >  1:2865622/1‑49 (MQ=255)
                             atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                                   >  1:657575/1‑69 (MQ=255)
                             atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                                   >  1:3282971/1‑69 (MQ=255)
                             atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                                   >  1:2601749/1‑69 (MQ=255)
                             atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                                   >  1:740937/1‑69 (MQ=255)
                             atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                                   >  1:787211/1‑69 (MQ=255)
                             atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                                   >  1:1978515/1‑69 (MQ=255)
                             atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGa                                    >  1:2271410/1‑68 (MQ=255)
                                gATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAActc                                <  1:3073289/69‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:3310504/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:1476745/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:1428525/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:573600/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:1787586/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:3150546/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:1335138/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:756469/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:826054/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:1178271/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:3084247/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2506879/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2960758/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2946853/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2930166/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2778577/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:1985899/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2641055/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2050408/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2711643/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2688792/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2094606/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:2307810/54‑1 (MQ=255)
                                 aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                              <  1:222063/54‑1 (MQ=255)
                                                        cgcgATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCaga         >  1:2716268/1‑68 (MQ=255)
                                                          cgATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCagaaga      <  1:2998952/69‑1 (MQ=255)
                                                          cgATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCagaaga      <  1:2108712/69‑1 (MQ=255)
                                                          cgATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCagaaga      <  1:450196/69‑1 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAAt   >  1:477203/1‑68 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAAt   >  1:559409/1‑68 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAAt   >  1:563359/1‑68 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAAt   >  1:3432458/1‑68 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAAt   >  1:330662/1‑68 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAAt   >  1:2095384/1‑68 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAAt   >  1:3266613/1‑68 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAAt   >  1:2357999/1‑68 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAATg  <  1:232271/69‑1 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAATg  <  1:228358/69‑1 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAATg  <  1:2617482/69‑1 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAATg  <  1:556518/69‑1 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAATg  <  1:2995609/69‑1 (MQ=255)
                                                              aaCGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACAGAAAGATCAGCCAGAAGAAAt   >  1:2520114/1‑68 (MQ=255)
                                                               aCGCCGAACGTAATGGCGTTTCTTACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAATg  <  1:1254990/68‑1 (MQ=255)
                                                                 gCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCagaaga      <  1:2092557/62‑1 (MQ=255)
                                                                   |                                                               
TGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGAAATG  >  minE/2224887‑2225017

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: