Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 528,351 T→C 18.6% S201S (TCA→TCG ybhN ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE528,3510TC18.6% 119.5 / 17.3 80S201S (TCA→TCGybhNconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base C (7/8);  total (40/40)
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

ATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAAC  >  minE/528298‑528414
                                                     |                                                               
aTTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGc                                                  <  1:2592601/69‑1 (MQ=255)
   aCCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCttt                                               >  1:74162/1‑69 (MQ=255)
   aCCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCttt                                               >  1:1298297/1‑69 (MQ=255)
   aCCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCttt                                               >  1:133982/1‑69 (MQ=255)
   aCCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCttt                                               >  1:2932365/1‑69 (MQ=255)
   aCCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCttt                                               >  1:2345295/1‑69 (MQ=255)
   aCCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCttt                                               >  1:200096/1‑69 (MQ=255)
   aCCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCtt                                                >  1:2998600/1‑68 (MQ=255)
       tCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGcc                                           <  1:852288/69‑1 (MQ=255)
       tCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGcc                                           <  1:3264466/69‑1 (MQ=255)
       tCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGcc                                           <  1:3024811/69‑1 (MQ=255)
       tCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGcc                                           <  1:3011210/69‑1 (MQ=255)
        ccAGTTAACGCTGGAAATCAGCATCTGAGCTAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGccc                                          >  1:21142/1‑69 (MQ=255)
        ccAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGGATTTCCATGAAGGCAGCACCAGCTTTTGccc                                          >  1:2500410/1‑69 (MQ=255)
        ccAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGgcag                                                         >  1:3051920/1‑54 (MQ=255)
        ccAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCCCCAGCTTTTGccc                                          >  1:839411/1‑69 (MQ=255)
        ccAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGccc                                          >  1:471772/1‑69 (MQ=255)
         cAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCt                                         >  1:1225795/1‑69 (MQ=255)
         cAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCt                                         >  1:2504188/1‑69 (MQ=255)
         cAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCt                                         >  1:672330/1‑69 (MQ=255)
            ttAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGGCCTTTg                                      >  1:710281/1‑69 (MQ=255)
            ttAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTg                                      >  1:658619/1‑69 (MQ=255)
            ttAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTg                                      >  1:1487979/1‑69 (MQ=255)
            ttAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTg                                      >  1:1486136/1‑69 (MQ=255)
            ttAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTg                                      >  1:1894406/1‑69 (MQ=255)
            ttAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTg                                      >  1:1284477/1‑69 (MQ=255)
            ttAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTg                                      >  1:1630122/1‑69 (MQ=255)
            ttAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTg                                      >  1:303552/1‑69 (MQ=255)
              aaCGCTGGAAATCAGCCTCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCTGCTTTTGCCCTTTGa                                     >  1:61322/1‑68 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCCTGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:256673/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGTTGGTCa                               <  1:438661/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2849918/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2781327/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2402503/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2269022/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:1547189/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:3069460/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2224890/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2186541/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:682021/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2076031/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:1956723/69‑1 (MQ=255)
                   tGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:172897/69‑1 (MQ=255)
                     gAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:3203217/67‑1 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:1344906/62‑1 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:3353835/62‑1 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2933874/62‑1 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2621032/62‑1 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:453125/62‑1 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:243025/62‑1 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:3341414/62‑1 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCg                        >  1:3234268/1‑69 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCg                        >  1:1146922/1‑69 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCg                        >  1:13041/1‑69 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCg                        >  1:2926078/1‑69 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCg                        >  1:2313683/1‑69 (MQ=255)
                          cAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCg                        >  1:1826368/1‑69 (MQ=255)
                                 tGAGCGATGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:516255/55‑1 (MQ=255)
                                 tGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:524301/55‑1 (MQ=255)
                                 tGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:698626/55‑1 (MQ=255)
                                 tGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:2843473/55‑1 (MQ=255)
                                  gAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCa                               <  1:1063043/54‑1 (MQ=255)
                                     tgAGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTc                 <  1:2515654/64‑1 (MQ=255)
                                     tgAGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTc                 <  1:2498643/64‑1 (MQ=255)
                                     tgAGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTc                 <  1:498412/64‑1 (MQ=255)
                                     tgAGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTc                 <  1:2463527/64‑1 (MQ=255)
                                     tgAGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTc                 <  1:2459631/64‑1 (MQ=255)
                                     tgAGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTc                 <  1:1889508/64‑1 (MQ=255)
                                     tgAGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTc                 <  1:1919953/64‑1 (MQ=255)
                                     tgAGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTc                 <  1:2246798/64‑1 (MQ=255)
                                       aGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGaaa           >  1:1636798/1‑69 (MQ=255)
                                       aGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGaaa           >  1:796825/1‑69 (MQ=255)
                                       aGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGaaa           >  1:1820253/1‑69 (MQ=255)
                                       aGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGaaa           >  1:1575856/1‑69 (MQ=255)
                                       aGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGaaa           >  1:2326568/1‑69 (MQ=255)
                                       aGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGaaa           >  1:2374790/1‑69 (MQ=255)
                                       aGGGCGAATTTCCACGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGaa            >  1:1227580/1‑68 (MQ=255)
                                                 tCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAAc  >  1:521836/1‑68 (MQ=255)
                                                 tCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAAc  >  1:1763535/1‑68 (MQ=255)
                                                 tCCATGAAGGAAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAAc  >  1:1646224/1‑68 (MQ=255)
                                                     |                                                               
ATTACCATCCAGTTAACGCTGGAAATCAGCATCTGAGCGAGGGCGAATTTCCATGAAGGCAGCACCAGCTTTTGCCCTTTGATGGTCATATGGCGGTGCTTCGCGAAAGCGCAAAAC  >  minE/528298‑528414

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: