Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,019,133 C→T 35.6% P340L (CCG→CTG)  yqhD → alcohol dehydrogenase, NAD(P)‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,019,1330CT35.6% 25.7 / 20.9 42P340L (CCG→CTG) yqhDalcohol dehydrogenase, NAD(P)‑dependent
Reads supporting (aligned to +/- strand):  ref base C (16/11);  new base T (2/13);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 8.07e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.62e-01

TCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAG  >  minE/2019066‑2019200
                                                                   |                                                                   
tCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCg                                                                    >  1:3263226/1‑69 (MQ=255)
 ccGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGa                                                                   >  1:2951231/1‑69 (MQ=255)
 ccGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGa                                                                   >  1:466863/1‑69 (MQ=255)
 ccGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGa                                                                   >  1:446297/1‑69 (MQ=255)
 ccGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGa                                                                   >  1:1957186/1‑69 (MQ=255)
 ccGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGa                                                                   >  1:3421982/1‑69 (MQ=255)
 ccGATGATGAGCGTATTGACACCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGa                                                                   >  1:2619560/1‑69 (MQ=255)
            cGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTcc                                                        <  1:1530647/69‑1 (MQ=255)
            cGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTcc                                                        <  1:2273515/69‑1 (MQ=255)
             gcaTTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:1584967/67‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTTCTGACCCACCTCTCCg                                                       <  1:2118806/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:403755/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:2921929/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:3118372/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:383833/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:2822356/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:548003/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:655346/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:968352/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:2433890/69‑1 (MQ=255)
             gTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:2116149/69‑1 (MQ=255)
               aTTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCg                                                       <  1:1224791/67‑1 (MQ=255)
               aTTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGAc                                                     <  1:3112798/69‑1 (MQ=255)
               aTTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGAc                                                     <  1:2061088/69‑1 (MQ=255)
                  gACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACct                                                            >  1:50907/1‑59 (MQ=255)
                  gACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGa                                                      >  1:2617603/1‑65 (MQ=255)
                  gACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTAc                                                  >  1:2720402/1‑69 (MQ=255)
                        gcgATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTg                                            <  1:693346/69‑1 (MQ=255)
                                   aaCCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTc                                 <  1:1256899/69‑1 (MQ=255)
                                             ttCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCt                       >  1:2312092/1‑69 (MQ=255)
                                             ttCTTTGAGCAATTAGGCGTGCTGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCt                       >  1:3324303/1‑69 (MQ=255)
                                              tCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCtt                      >  1:1110961/1‑69 (MQ=255)
                                              tCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCtt                      >  1:2935174/1‑69 (MQ=255)
                                                           aGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAAc          >  1:2138891/1‑68 (MQ=255)
                                                           aGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAAc          >  1:2148521/1‑68 (MQ=255)
                                                           aGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAAc          >  1:2260715/1‑68 (MQ=255)
                                                           aGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAAc          >  1:1065368/1‑68 (MQ=255)
                                                              cGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTgg       <  1:371482/68‑1 (MQ=255)
                                                              cGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTgg       <  1:1135759/68‑1 (MQ=255)
                                                                tGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAag    <  1:322146/69‑1 (MQ=255)
                                                                tGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAag    <  1:76343/69‑1 (MQ=255)
                                                                 gCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAag    <  1:2471146/68‑1 (MQ=255)
                                                                  ccGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAagag  <  1:1430244/69‑1 (MQ=255)
                                                                  ccGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAagag  <  1:671737/69‑1 (MQ=255)
                                                                   |                                                                   
TCCGATGATGAGCGTATTGACGCCGCGATTGCCGCAACCCGCAATTTCTTTGAGCAATTAGGCGTGCCGACCCACCTCTCCGACTACGGTCTGGACGGCAGCTCCATCCCGGCTTTGCTGAAAAAACTGGAAGAG  >  minE/2019066‑2019200

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: