Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 289,285 G→A 26.1% T643T (ACG→ACA lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE289,2850GA26.1% 70.2 / 20.4 61T643T (ACG→ACAlonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base G (25/20);  new base A (8/8);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 7.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

AAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCA  >  minE/289218‑289350
                                                                   |                                                                 
aaGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACgg                                                                  >  1:836847/1‑69 (MQ=255)
aaGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACgg                                                                  >  1:399938/1‑69 (MQ=255)
           aCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcg                                                       >  1:2428942/1‑69 (MQ=255)
           aCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcg                                                       >  1:2834026/1‑69 (MQ=255)
           aCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcg                                                       >  1:821351/1‑69 (MQ=255)
           aCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgc                                                        >  1:905907/1‑68 (MQ=255)
           aCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgc                                                        >  1:2478921/1‑68 (MQ=255)
           aCCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgc                                                        >  1:2845227/1‑68 (MQ=255)
             cTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcg                                                       <  1:2735318/67‑1 (MQ=255)
             cTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcg                                                       <  1:1243243/67‑1 (MQ=255)
             cTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcg                                                       <  1:2322496/67‑1 (MQ=255)
             cTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcg                                                       <  1:2200297/67‑1 (MQ=255)
             cTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcg                                                       <  1:310870/67‑1 (MQ=255)
                taCCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcgcgtg                                                   >  1:2225484/1‑68 (MQ=255)
                  ccGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcgcgtgcg                                                 >  1:166754/1‑68 (MQ=255)
                  ccGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTcgtgcgcgtgc                                                  >  1:2027336/1‑67 (MQ=255)
                                gAAGTGATGCAGGAGTCTATTCAGGCGGCGTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTggg                                      >  1:2077347/1‑65 (MQ=255)
                                gAAGTGATGCAGGAGTCTATTCAGGCGGCGTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTggg                                      >  1:3329285/1‑65 (MQ=255)
                                gAAGTGATGCAGGAGTCTATTCAGGCGGCGTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTggg                                      >  1:2939539/1‑65 (MQ=255)
                                gAAGTGATGCAGGAGTCTATTCAGGCGGCGTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTggg                                      >  1:764517/1‑65 (MQ=255)
                                gAAGTGATGCAGGAGTCTATTCAGGCGGCGTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTggg                                      >  1:546090/1‑65 (MQ=255)
                                gAAGTGATGCAGGAGTCTATTCAGGCGGCGTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTggg                                      >  1:880730/1‑65 (MQ=255)
                                gAAGTGATGCAGGAGTCTATTCAGGCGGCGTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTggg                                      >  1:1670946/1‑65 (MQ=255)
                                gAAGTGATGCAGGAGTCTATTCAGGCGGCGTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTggg                                      >  1:1148358/1‑65 (MQ=255)
                                     gATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAAcc                              >  1:732517/1‑68 (MQ=255)
                                      aTGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCt                            >  1:3241194/1‑69 (MQ=255)
                                      aTGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCt                            >  1:3382738/1‑69 (MQ=255)
                                      aTGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCt                            >  1:27390/1‑69 (MQ=255)
                                      aTGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCt                            >  1:1551901/1‑69 (MQ=255)
                                              gTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGAttttt                     >  1:1534370/1‑68 (MQ=255)
                                              gTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGAttttt                     >  1:2201664/1‑68 (MQ=255)
                                                        gcggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:1141961/69‑1 (MQ=255)
                                                        gcggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:2770340/69‑1 (MQ=255)
                                                        gcggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:782128/69‑1 (MQ=255)
                                                        gcggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:458002/69‑1 (MQ=255)
                                                        gcggcgTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:183888/69‑1 (MQ=255)
                                                        gcggcgTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:1364912/69‑1 (MQ=255)
                                                        gcggcgTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:1525682/69‑1 (MQ=255)
                                                        gcggcgTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:3420420/69‑1 (MQ=255)
                                                        gcggcgTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:2108357/69‑1 (MQ=255)
                                                        gcggcgTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:1042109/69‑1 (MQ=255)
                                                        gcggcgTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:233016/69‑1 (MQ=255)
                                                        gcggcgTTAACAGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:3097220/69‑1 (MQ=255)
                                                         cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGt          <  1:3182216/68‑1 (MQ=255)
                                                         cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTg         <  1:670724/69‑1 (MQ=255)
                                                         cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTg         <  1:1465356/69‑1 (MQ=255)
                                                         cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTg         <  1:757650/69‑1 (MQ=255)
                                                         cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTg         <  1:744014/69‑1 (MQ=255)
                                                         cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTg         <  1:719500/69‑1 (MQ=255)
                                                         cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTg         <  1:3155109/69‑1 (MQ=255)
                                                         cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTg         <  1:572252/69‑1 (MQ=255)
                                                         cggcgTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTg         <  1:1783958/69‑1 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGCCAtcca  >  1:2198268/1‑69 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  >  1:1733049/1‑69 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  >  1:730344/1‑69 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  >  1:2675253/1‑69 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  >  1:3123679/1‑69 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  <  1:2378701/69‑1 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  >  1:280676/1‑69 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  >  1:2195499/1‑69 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  >  1:3189261/1‑69 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  >  1:244045/1‑69 (MQ=255)
                                                                aaCGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACAtcca  <  1:2694300/69‑1 (MQ=255)
                                                                   |                                                                 
AAGGCAAACTGACCTATACCGGTTCGCTCGGCGAAGTGATGCAGGAGTCCATTCAGGCGGCGTTAACGGTGGTTCGTGCGCGTGCGGAAAAACTGGGGATCAACCCTGATTTTTACGAAAAACGTGACATCCA  >  minE/289218‑289350

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: