Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 223,577 A→G 50.7% G151G (GGA→GGG prfH → predicted peptide chain release factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE223,5770AG50.7% ‑1.6 / 61.3 69G151G (GGA→GGGprfHpredicted peptide chain release factor
Reads supporting (aligned to +/- strand):  ref base A (10/24);  new base G (10/25);  total (20/49)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTA  >  minE/223522‑223639
                                                       |                                                              
aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGa                                                   <  1:564295/69‑1 (MQ=255)
aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGa                                                   <  1:242571/69‑1 (MQ=255)
aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGa                                                   <  1:3200656/69‑1 (MQ=255)
aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGa                                                   <  1:1970346/69‑1 (MQ=255)
aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGa                                                   <  1:1511878/69‑1 (MQ=255)
aGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGa                                                   <  1:1470666/69‑1 (MQ=255)
               aCCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCg                                    <  1:2592441/69‑1 (MQ=255)
               aCCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCg                                    <  1:3427487/69‑1 (MQ=255)
               aCCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCg                                    <  1:1264747/69‑1 (MQ=255)
               aCCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCg                                    <  1:798941/69‑1 (MQ=255)
                   gcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                <  1:803609/69‑1 (MQ=255)
                   gcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                <  1:872464/69‑1 (MQ=255)
                   gcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                <  1:1672686/69‑1 (MQ=255)
                   gcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                <  1:185654/69‑1 (MQ=255)
                   gcgACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                <  1:33629/69‑1 (MQ=255)
                      aCGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                <  1:670004/66‑1 (MQ=255)
                               tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATcc                      >  1:171233/1‑67 (MQ=255)
                               tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATcc                      >  1:2093337/1‑67 (MQ=255)
                               tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATcc                      >  1:2375438/1‑67 (MQ=255)
                               tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATcc                      >  1:413148/1‑67 (MQ=255)
                               tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATcc                      >  1:719730/1‑67 (MQ=255)
                               tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATcc                      >  1:2561248/1‑67 (MQ=255)
                               tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATcc                      >  1:3362934/1‑67 (MQ=255)
                               tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATcc                      >  1:1467085/1‑67 (MQ=255)
                               tACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATc                       >  1:2104993/1‑66 (MQ=255)
                                   gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcag                 <  1:3315435/68‑1 (MQ=255)
                                   gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcag                 <  1:2296215/68‑1 (MQ=255)
                                   gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcag                 <  1:2899799/68‑1 (MQ=255)
                                   gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcag                 <  1:2173538/68‑1 (MQ=255)
                                   gggATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcag                 <  1:2142724/68‑1 (MQ=255)
                                      aTGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCata             <  1:2698317/69‑1 (MQ=255)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcatcaaa             <  1:379965/69‑3 (MQ=38)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:554427/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:3256902/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:2276868/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:2681157/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:3339109/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:568183/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:1021494/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:2446904/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:222986/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:2068552/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:1219769/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:1330902/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:754354/69‑3 (MQ=39)
                                      aTGGCCTTTATCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGgcagcaaa             <  1:145389/69‑3 (MQ=39)
                                         gCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTa                                >  1:2950743/1‑47 (MQ=255)
                                            tttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGtcagca               <  1:3064752/61‑1 (MQ=255)
                                            tttATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGtcagca               <  1:3068645/61‑1 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat    >  1:2030019/1‑69 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat    >  1:1125592/1‑69 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat    >  1:1611122/1‑69 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat    >  1:2026130/1‑69 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat    >  1:2053042/1‑69 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataat    >  1:2656521/1‑69 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataa     >  1:1979153/1‑68 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataa     >  1:808275/1‑68 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataa     >  1:3263737/1‑68 (MQ=255)
                                               aTCGAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAaataa     >  1:1323228/1‑68 (MQ=255)
                                                  gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCCAATAATTa  <  1:1038117/68‑1 (MQ=39)
                                                  gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa  <  1:515245/68‑1 (MQ=255)
                                                  gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa  <  1:2636626/68‑1 (MQ=255)
                                                  gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa  <  1:2477807/68‑1 (MQ=255)
                                                  gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa  <  1:1430513/68‑1 (MQ=255)
                                                  gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa  <  1:1359147/68‑1 (MQ=255)
                                                  gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa  <  1:1293488/68‑1 (MQ=255)
                                                  gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa  <  1:878786/68‑1 (MQ=255)
                                                  gAAGGGTAATGAAGGATGAAACTGCCGGAAGGCGAATAAACGCCATCCGGCAGCTAATCAAATAATTa  <  1:3414537/68‑1 (MQ=255)
                                                   aaGGGTAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCAAATCAAATAATTa  <  1:1603131/67‑1 (MQ=255)
                                                       |                                                              
AGATTGAACGTGGCAACCCGCGACGGACATTTACAGGGATGGCTTTTATCGAAGGATAATGAAGGATGAAACTGCCGGAAGGCGATTAAACGCCATCCGGCAGCATATCAAATAATTA  >  minE/223522‑223639

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: