Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 384,346 C→T 32.1% F50F (TTC→TTT ybdB → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE384,3460CT32.1% 50.2 / 27.3 59F50F (TTC→TTTybdBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/20);  new base T (13/6);  total (33/26)
Fisher's exact test for biased strand distribution p-value = 2.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGTCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGAT  >  minE/384282‑384414
                                                                |                                                                    
cccGTCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGtt                                                                  >  1:3013303/1‑69 (MQ=255)
cccGTCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGtt                                                                  >  1:2742379/1‑69 (MQ=255)
cccGTCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGtt                                                                  >  1:2325104/1‑69 (MQ=255)
cccGTCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGtt                                                                  >  1:134329/1‑69 (MQ=255)
cccGTCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGtt                                                                  >  1:2213422/1‑69 (MQ=255)
cccGTCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGtt                                                                  >  1:1678400/1‑69 (MQ=255)
cccGTCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGt                                                                   >  1:3089879/1‑68 (MQ=255)
    tCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtact                                                              <  1:1697030/69‑1 (MQ=255)
     cttggCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:1060457/66‑1 (MQ=255)
     cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtactg                                                             <  1:2664233/69‑2 (MQ=255)
     cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:110161/69‑1 (MQ=255)
     cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:2161729/69‑1 (MQ=255)
     cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:1357877/69‑1 (MQ=255)
     cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:1410400/69‑1 (MQ=255)
     cTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTtacta                                                             <  1:1305483/69‑1 (MQ=255)
            atgatgTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATggcg                                                      >  1:621782/1‑69 (MQ=255)
            atgatgTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATggcg                                                      >  1:3185648/1‑69 (MQ=255)
                                  gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:1624058/69‑1 (MQ=255)
                                  gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:2352750/69‑1 (MQ=255)
                                  gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:25466/69‑1 (MQ=255)
                                  gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:1294107/69‑1 (MQ=255)
                                  gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:2796177/69‑1 (MQ=255)
                                  gCCGGTTGATACCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                <  1:2895583/69‑1 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGTCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:1311923/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGa                                >  1:2864896/1‑59 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2932245/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:3020643/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2566593/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:698418/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:820897/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:164285/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:1770295/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2179374/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2048183/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2108540/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTTGGCCTGCTGCATGGCGGCGCGTAGGCGGCGCTGGCGGAAACGCTGGGa                      >  1:2052158/1‑69 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:1911974/68‑1 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:996928/68‑1 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:3398627/68‑1 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:1144426/68‑1 (MQ=255)
                                            aCCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTggg                       <  1:2465067/68‑1 (MQ=255)
                                                   cTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATgg               >  1:2070544/1‑69 (MQ=255)
                                                   cTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATgg               >  1:1710759/1‑69 (MQ=255)
                                                   cTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATgg               >  1:57925/1‑69 (MQ=255)
                                                     catcaGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGc              >  1:1276006/1‑68 (MQ=255)
                                                     catcaGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGc              >  1:2233600/1‑68 (MQ=255)
                                                     catcaGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGc              >  1:914172/1‑68 (MQ=255)
                                                     catcaGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGc              >  1:902590/1‑68 (MQ=255)
                                                        caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:1905624/38‑1 (MQ=255)
                                                        caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:1364017/38‑1 (MQ=255)
                                                        caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:2125408/38‑1 (MQ=255)
                                                        caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:65806/38‑1 (MQ=255)
                                                        caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:2799516/38‑1 (MQ=255)
                                                        caGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:898145/38‑1 (MQ=255)
                                                         aGCCGTTCGGTTTACTACATGGCGGCGCGTCggcggc                                         <  1:2978114/37‑1 (MQ=255)
                                                            cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGAttt       >  1:3247668/1‑68 (MQ=255)
                                                            cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGAttt       >  1:610729/1‑68 (MQ=255)
                                                            cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGAttt       >  1:1425223/1‑68 (MQ=255)
                                                            cGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGAttt       >  1:918300/1‑68 (MQ=255)
                                                                cGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTatgat  >  1:2420593/1‑69 (MQ=255)
                                                                cGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTatgat  >  1:2284065/1‑69 (MQ=255)
                                                                cGGTTTACTACATGGCGGCGCGTAGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTatgat  >  1:2196675/1‑69 (MQ=255)
                                                                |                                                                    
CCCGTCTGGGCGATGATGTGCTGGAAGCCGAAATGCCGGTTGATACCCGTACTCATCAGCCGTTCGGTTTACTACATGGCGGCGCGTCGGCGGCGCTGGCGGAAACGCTGGGATCGATGGCCGGATTTATGAT  >  minE/384282‑384414

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: