Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 323,771 C→T 50.0% G61G (GGC→GGT htpG → molecular chaperone HSP90 family

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE323,7710CT50.0% ‑2.0 / 53.0 60G61G (GGC→GGThtpGmolecular chaperone HSP90 family
Reads supporting (aligned to +/- strand):  ref base C (20/10);  new base T (19/11);  total (39/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

CCTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCGCTCTCTAACCCGGACCTGTACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCGATAACG  >  minE/323704‑323832
                                                                   |                                                             
ccTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCGCTCTCTAACCCGGACCTGTACGAAGGTGATGGCg                                                              <  1:3039988/69‑1 (MQ=255)
ccTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCGCTCTCTAACCCGGACCTGTACGAAGGTGATGGCg                                                              <  1:2368393/69‑1 (MQ=255)
ccTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCGCTCTCTAACCCGGACCTGTACGAAGGTGATGGCg                                                              <  1:2343785/69‑1 (MQ=255)
ccTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCGCTCTCTAACCCGGACCTGTACGAAGGTGATGGCg                                                              <  1:2169860/69‑1 (MQ=255)
ccTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCGCTCTCTAACCCGGACCTGTACGAAGGTGATGGCg                                                              <  1:1634815/69‑1 (MQ=255)
ccTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCGCTCTCTAACCCGGACCTGTACGAAGGTGATGGCg                                                              <  1:511975/69‑1 (MQ=255)
 cTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCGCTCTCTAACCCGGACCTGTACGAAGGTGATGGCg                                                              <  1:2105823/68‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:3162715/69‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:2905506/69‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:2794486/69‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:220360/69‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:486067/69‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:1974789/69‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:1890602/69‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:16873/69‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:773479/69‑1 (MQ=255)
                           tCCGCGCGCTTTCTAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGCCTCTTTCGATa                                   <  1:2590564/69‑1 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGAt                                    >  1:1213780/3‑59 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:540227/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:3150608/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:3014836/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:2920996/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:3286877/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:2691383/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:2663176/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:2658311/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:1227794/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:2318604/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:551045/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:2190357/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:581491/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:210505/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:2032354/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:567301/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAg                           >  1:1725471/3‑68 (MQ=255)
                                    tttctAACCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTAGTGTCTCTTTCGATAAAGACAAg                           >  1:3429597/3‑68 (MQ=255)
                                     ctctAACCCGGACCTGTACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCg                         >  1:643241/1‑69 (MQ=255)
                                     ctctAACCCGGACCTGTACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCg                         <  1:2468973/69‑1 (MQ=255)
                                     ctctAACCCGGACCTGTACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCg                         >  1:2115077/1‑69 (MQ=255)
                                     ctctAACCCGGACCTGTACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCg                         >  1:1100042/1‑69 (MQ=255)
                                     ctctAACCCGGACCTGTACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCg                         <  1:1948892/69‑1 (MQ=255)
                                      tctAACCCGGACCTGTACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCg                         <  1:223146/68‑1 (MQ=255)
                                         aaCCCGGACCTGTACGAAGGTGATGGTGAACTACGCGTTCGTGTCTCTTTCGATa                                   <  1:1174087/55‑1 (MQ=255)
                                                     tACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGCCAAGCGTACGCTGACCATCTcc         >  1:2728936/1‑69 (MQ=255)
                                                     tACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTcc         >  1:3265303/1‑69 (MQ=255)
                                                     tACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTcc         >  1:399364/1‑69 (MQ=255)
                                                     tACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTcc         >  1:1113515/1‑69 (MQ=255)
                                                     tACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTcc         >  1:2856317/1‑69 (MQ=255)
                                                     tACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTcc         >  1:2799817/1‑69 (MQ=255)
                                                     tACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTcc         >  1:1906541/1‑69 (MQ=255)
                                                     tACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTcc         >  1:1738754/1‑69 (MQ=255)
                                                       cGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCg        >  1:1948099/1‑68 (MQ=255)
                                                       cGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCg        >  1:1690935/1‑68 (MQ=255)
                                                       cGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCg        >  1:1593239/1‑68 (MQ=255)
                                                       cGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCg        >  1:1251114/1‑68 (MQ=255)
                                                       cGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCg        >  1:942203/1‑68 (MQ=255)
                                                            gTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCGATAACg  >  1:3091112/1‑69 (MQ=255)
                                                            gTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCGATAACg  >  1:2241755/1‑69 (MQ=255)
                                                            gTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCGATAACg  >  1:1979062/1‑69 (MQ=255)
                                                            gTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCGATAACg  >  1:1600924/1‑69 (MQ=255)
                                                                   |                                                             
CCTCCGATGCGGCGGACAAGCTGCGTTTCCGTGCGCTCTCTAACCCGGACCTGTACGAAGGTGATGGCGAACTACGCGTTCGTGTCTCTTTCGATAAAGACAAGCGTACGCTGACCATCTCCGATAACG  >  minE/323704‑323832

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: