Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,224,948 T→C 57.4% D199D (GAT→GAC 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,9480TC57.4% 25.7 / 62.0 82D199D (GAT→GACprmAmethylase for 50S ribosomal subunit protein L11
Reads supporting (aligned to +/- strand):  ref base T (19/16);  new base C (17/30);  total (36/46)
Fisher's exact test for biased strand distribution p-value = 1.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGA  >  minE/2224881‑2225013
                                                                   |                                                                 
tGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGCCCAGCCGCGAt                                                                   <  1:2940307/68‑1 (MQ=255)
 gAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   <  1:807084/67‑1 (MQ=255)
 gAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   <  1:3059381/67‑1 (MQ=255)
 gAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   <  1:1369559/67‑1 (MQ=255)
 gAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   <  1:3215457/67‑1 (MQ=255)
 gAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   <  1:1404458/67‑1 (MQ=255)
 gAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATaa                                                                 <  1:791071/69‑1 (MQ=255)
 gAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATaa                                                                 <  1:1284749/69‑1 (MQ=255)
 gAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATaa                                                                 <  1:1962306/69‑1 (MQ=255)
  aaaCTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATaa                                                                 <  1:2741754/68‑1 (MQ=255)
      tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                            >  1:2322007/1‑69 (MQ=255)
      tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                            >  1:2772431/1‑69 (MQ=255)
      tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                            >  1:3157534/1‑69 (MQ=255)
      tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                            >  1:705883/1‑69 (MQ=255)
      tGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCg                                                            >  1:1203342/1‑69 (MQ=255)
          tGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGaa                                                          <  1:2572042/67‑1 (MQ=255)
          tGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGaa                                                          <  1:592433/67‑1 (MQ=255)
           gcagcaAAAGCCATTGTTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                       <  1:3259328/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                       <  1:2718619/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                       <  1:2139701/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                       <  1:3096298/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                       <  1:3404743/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGt                                                       <  1:1006006/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                       <  1:3312036/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                       <  1:2654362/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                       <  1:564507/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                       <  1:2515695/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                       <  1:887969/69‑1 (MQ=255)
           gcagcaAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGt                                                       <  1:1217643/69‑1 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAAc                                                                >  1:360222/1‑58 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:1452275/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:737843/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:2668988/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:2117780/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:1611452/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:2628509/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:419927/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:2594231/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:2986724/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:2812545/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:1183518/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:1360518/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:565986/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:2141375/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:1029492/1‑69 (MQ=255)
             agcaAAAGCCATTGGTATTGATATTGATCCGCAGGCGATTAAGGCCAGCCGCGACAACGCTGAACGTaa                                                     >  1:1428918/1‑69 (MQ=255)
                 aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   >  1:1643738/1‑51 (MQ=255)
                 aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   >  1:3123034/1‑51 (MQ=255)
                 aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   >  1:1389754/1‑51 (MQ=255)
                 aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   >  1:785658/1‑51 (MQ=255)
                 aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   >  1:1443255/1‑51 (MQ=255)
                 aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   >  1:3234379/1‑51 (MQ=255)
                 aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGAt                                                                   >  1:3246242/1‑51 (MQ=255)
                 aaaGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGc                                                 >  1:3101952/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)
                          ggTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGAc                                        >  1:760232/1‑69 (MQ=255)
                                   atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATgcg                                                 >  1:2865622/1‑49 (MQ=255)
                                   atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                               >  1:657575/1‑69 (MQ=255)
                                   atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                               >  1:1978515/1‑69 (MQ=255)
                                   atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                               >  1:740937/1‑69 (MQ=255)
                                   atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                               >  1:3282971/1‑69 (MQ=255)
                                   atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                               >  1:2601749/1‑69 (MQ=255)
                                   atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGaa                               >  1:787211/1‑69 (MQ=255)
                                   atCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGa                                >  1:2271410/1‑68 (MQ=255)
                                      gATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAActc                            <  1:3073289/69‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:3084247/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:573600/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:826054/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:1178271/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:756469/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:1335138/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2778577/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2094606/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2506879/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2050408/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:1985899/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2641055/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2688792/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2711643/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:1428525/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:1787586/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2930166/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2946853/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2960758/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:2307810/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:222063/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:1476745/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:3150546/54‑1 (MQ=255)
                                       aTCCGCAGGCGATTCAGGCCAGCCGCGACAACGCTGAACGTAATGGCGTTTCTg                                          <  1:3310504/54‑1 (MQ=255)
                                                              cgcgATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCaga     >  1:2716268/1‑68 (MQ=255)
                                                                cgATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCagaaga  <  1:450196/69‑1 (MQ=255)
                                                                cgATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCagaaga  <  1:2998952/69‑1 (MQ=255)
                                                                cgATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCagaaga  <  1:2108712/69‑1 (MQ=255)
                                                                   |                                                                 
TGAAACTGGGTGCAGCAAAAGCCATTGGTATTGATATCGATCCGCAGGCGATTCAGGCCAGCCGCGATAACGCCGAACGTAATGGCGTTTCTGACCGTCTGGAACTCTACTTACCGAAAGATCAGCCAGAAGA  >  minE/2224881‑2225013

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: