Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,377,802 A→T 20.0% Y196N (TAC→AAC)  ubiB ← 2‑octaprenylphenol hydroxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,377,8020AT20.0% 49.9 / 9.3 25Y196N (TAC→AAC) ubiB2‑octaprenylphenol hydroxylase
Reads supporting (aligned to +/- strand):  ref base A (9/11);  new base T (2/3);  total (11/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

TTGCGCCGAAGCTGAATGGCGTTGGCAGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTGGGCGCAGACGGCGACCATCCGGCAGCAAACGCGGCACCCAGCGAG  >  minE/2377732‑2377867
                                                                      |                                                                 
ttGCGCCGAAGCTGAATGGCGTTGGCAGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTa                                                                   >  1:2098347/1‑71 (MQ=255)
  gcgcCGAAGCTGAATGGCGTTGGCAGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTACt                                                                 >  1:2041645/1‑71 (MQ=255)
   cgcCGAAGCTGAATGGCGTTGGCAGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTACTc                                                                <  1:2081741/71‑1 (MQ=255)
               aTGGCGTTGGCAGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTc                                                    <  1:1863551/71‑1 (MQ=255)
                       ggCAGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTggg                                             <  1:379516/70‑1 (MQ=255)
                         cAGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTTCTCGCGCACCACTTCGGTTGGgcgc                                          <  1:2184799/71‑1 (MQ=255)
                         cAGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTTCTCGCGCACCACTTCGGTTGGgcgc                                          <  1:2115537/71‑1 (MQ=255)
                          aGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTTCTCGCGCACCACTTCGGTTGGgcgc                                          <  1:898756/70‑1 (MQ=255)
                          aGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTGGGCGCa                                         <  1:1328019/71‑1 (MQ=255)
                          aGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTGGGCGCa                                         <  1:1133777/71‑1 (MQ=255)
                              tCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTTCTCGCGCACCACTTCGGTTGGGCGCAGACg                                     >  1:826350/1‑71 (MQ=255)
                              tCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTGGGCGCAGa                                       >  1:1454965/1‑69 (MQ=255)
                                   cagcaAATTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTGGGCGCAGACGGCGAc                                <  1:106209/71‑1 (MQ=255)
                                     gcaAATTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTGGGCGCAGACGGCGACCa                              >  1:67880/1‑71 (MQ=255)
                                       aaaTTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTg                                               >  1:1770164/1‑52 (MQ=255)
                                         aTTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTc                                                    <  1:1953329/45‑1 (MQ=255)
                                         aTTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTc                                                    <  1:76808/45‑1 (MQ=255)
                                           tCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCgg                                                  >  1:373538/1‑45 (MQ=255)
                                           tCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCgg                                                  >  1:1180219/1‑45 (MQ=255)
                                               ttCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTgg                                              <  1:146359/45‑1 (MQ=255)
                                                     aatcaatGTCTTTTCGTACTCGCGCACCACTTCGGTTGGGCGCAGACGGCGACCATCCGGCAGCAAAcgcg              >  1:162744/1‑71 (MQ=255)
                                                       tcaatGTCTTTTCGTACTCGCGCACCACTTCGGtt                                                >  1:759873/1‑35 (MQ=255)
                                                             tCTTTTCGTACTCGCTCACCACTTCGGTTGGGCGCAGACGGCGACCATCCGGCAGCAAAcg                <  1:1682868/61‑1 (MQ=255)
                                                               ttttCGTACTCGCGCACCACTTCGGTTGGGCGCAGACGGCGACCATCCGGCAGCAAAcgcg              >  1:692672/1‑61 (MQ=255)
                                                                 ttCGTTCTCGCGCACCACTTCGGTTGGGCGCAGACGGCGACCATCCGGCAGCAAACGCGGCACCCAGCGAg  >  1:2117246/1‑71 (MQ=255)
                                                                  tCGTACTCGCGCACCACTTCGGGTGGGCGCAGACGGCGACCATCCGGCAGCAAACGCGGCACCCAGCg    <  1:786772/68‑1 (MQ=255)
                                                                      |                                                                 
TTGCGCCGAAGCTGAATGGCGTTGGCAGATTCCCGCAGCAAATTCAGTTCATCAATCAATGTCTTTTCGTACTCGCGCACCACTTCGGTTGGGCGCAGACGGCGACCATCCGGCAGCAAACGCGGCACCCAGCGAG  >  minE/2377732‑2377867

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: