Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 882,776 T→C 31.1% G175G (GGA→GGG yciT ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE882,7760TC31.1% 64.4 / 37.5 77G175G (GGA→GGGyciTpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (16/37);  new base C (18/6);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 4.08e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTTCGCACTCTTTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGC  >  minE/882711‑882840
                                                                 |                                                                
gCTTCGCACTCTTTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATcc                                                                <  1:2412897/68‑1 (MQ=255)
gCTTCGCACTCTTTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATcc                                                                <  1:3340295/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:1010099/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:784796/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:69574/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:137960/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:2495966/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:459376/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:1514158/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:2326010/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:2923349/68‑1 (MQ=255)
                 cAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:390475/68‑1 (MQ=255)
                   gCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:185564/66‑1 (MQ=255)
                    cACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATc                                           <  1:3188122/69‑1 (MQ=255)
                    cACGGCATTGACCACAGCGGTACGCATAATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATc                                           <  1:582535/69‑1 (MQ=255)
                             gACCACATCGGTACGCATCATATCGCGACCGGTAAAt                                                                  >  1:1726983/1‑37 (MQ=255)
                             gACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCaa                                   >  1:1065515/1‑68 (MQ=255)
                             gACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCaa                                   >  1:1697952/1‑68 (MQ=255)
                             gACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCaa                                   >  1:31994/1‑68 (MQ=255)
                             gACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCaa                                   >  1:1693976/1‑68 (MQ=255)
                             gACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCaa                                   >  1:1428836/1‑68 (MQ=255)
                             gACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCaa                                   >  1:3180819/1‑68 (MQ=255)
                             gACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCaa                                   >  1:1562939/1‑68 (MQ=255)
                             gACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCaa                                   >  1:1575209/1‑68 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:2968298/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:2503880/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:2537289/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:2579015/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:2834711/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:974034/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:3090639/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:3410868/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:3401787/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:238569/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:2390264/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:1757080/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:134211/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:135727/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:2034124/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:1736481/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTCAGGTTGCCAGCCATCaa                                         >  1:1583296/1‑56 (MQ=255)
                                   aTCGGTACGCATCATATCGCGACCGGTAAACCCAGTTTAAGGTTGCCAGCCATCaa                                         >  1:3032456/1‑56 (MQ=255)
                                    tCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTgccagc                                               <  1:2761269/49‑1 (MQ=255)
                                      ggTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:2780816/61‑1 (MQ=255)
                                      ggTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:3428937/61‑1 (MQ=255)
                                      ggTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:2844548/61‑1 (MQ=255)
                                      ggTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:2879901/61‑1 (MQ=255)
                                      ggTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:1692314/61‑1 (MQ=255)
                                      ggTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:2289034/61‑1 (MQ=255)
                                      ggTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:859656/61‑1 (MQ=255)
                                      ggTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:1984243/61‑1 (MQ=255)
                                            catcatATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGa                    <  1:1441254/68‑1 (MQ=255)
                                            catcatATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGa                    <  1:537119/68‑1 (MQ=255)
                                            catcatATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGa                    <  1:1090902/68‑1 (MQ=255)
                                            catcatATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGa                    <  1:2417662/68‑1 (MQ=255)
                                            catcatATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGa                    <  1:2517098/68‑1 (MQ=255)
                                            catcatATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGa                    <  1:2662308/68‑1 (MQ=255)
                                                atatCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:421878/51‑1 (MQ=255)
                                                atatCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGc               >  1:2308712/1‑69 (MQ=255)
                                                atatCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGc               >  1:601976/1‑69 (MQ=255)
                                                atatCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGc               >  1:512142/1‑69 (MQ=255)
                                                ataaCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:433353/51‑1 (MQ=255)
                                                    cgcgACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:2031669/47‑1 (MQ=255)
                                                     gcgACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTg          >  1:237953/1‑69 (MQ=255)
                                                     gcgACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTg          >  1:139695/1‑69 (MQ=255)
                                                     gcgACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTg          >  1:2971584/1‑69 (MQ=255)
                                                     gcgACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTg          >  1:2992283/1‑69 (MQ=255)
                                                     gcgACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTg          >  1:3095206/1‑69 (MQ=255)
                                                        aCCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATa                                 <  1:1935022/43‑1 (MQ=255)
                                                          cGGTAAACCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGa     <  1:1984387/69‑1 (MQ=255)
                                                          cGGTAAACCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGa     <  1:372278/69‑1 (MQ=255)
                                                          cGGTAAACCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGa     <  1:2213303/69‑1 (MQ=255)
                                                          cGGTAAACCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGa     <  1:2893994/69‑1 (MQ=255)
                                                           ggTAAACCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGa     <  1:2441510/68‑1 (MQ=255)
                                                           ggTAAACCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGa     <  1:1819615/68‑1 (MQ=255)
                                                            gTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATg   <  1:2247860/69‑1 (MQ=255)
                                                            gTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATg   <  1:423510/69‑1 (MQ=255)
                                                             tAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGc  <  1:2881698/69‑1 (MQ=255)
                                                                 |                                                                
GCTTCGCACTCTTTTTCCAGCACGGCATTGACCACATCGGTACGCATCATATCGCGACCGGTAAATCCAGTTTCAGGTTGCCAGCCATCAATACCAATAAATGCTTTGCTGAAATGCACCTGTTGGATGC  >  minE/882711‑882840

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: