Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,004,409 T→C 57.2% T396A (ACC→GCC)  pntA ← pyridine nucleotide transhydrogenase, alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,004,4090TC57.2% 4.8 / 17.2 28T396A (ACC→GCC) pntApyridine nucleotide transhydrogenase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base T (7/5);  new base C (10/6);  total (17/11)
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

CAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAA  >  minE/1004351‑1004477
                                                          |                                                                    
cAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTc                                                             >  1:2856016/1‑68 (MQ=255)
             tGATTGCCAACGCCATCAACGCGTATTTACGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACtt                                               <  1:1184351/69‑1 (MQ=255)
             tGATTGCCAACGCCATCAACGCGTATTTACGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACtt                                               <  1:1289870/69‑1 (MQ=255)
             tGATTGCCAACGCCATCAACGCGTATTTACGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACtt                                               <  1:1656535/69‑1 (MQ=255)
             tGATTGCCAACGCCATCAACGCGTATTTACGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACtt                                               <  1:3136136/69‑1 (MQ=255)
             tGATTGCCAACGCCATCAACGCGTATTTACGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACtt                                               <  1:2966570/69‑1 (MQ=255)
              gATTGCCAACGCCATCAACGCGTATTTACGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACtt                                               <  1:2786118/68‑1 (MQ=255)
                              aaCGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGcc                                       >  1:625316/1‑60 (MQ=255)
                              aaCGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTtgt                              >  1:1677835/1‑69 (MQ=255)
                              aaCGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTtgt                              >  1:1044220/1‑69 (MQ=255)
                              aaCGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTtgt                              >  1:906233/1‑69 (MQ=255)
                              aaCGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTtgt                              >  1:758719/1‑69 (MQ=255)
                              aaCGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTtgt                              >  1:248474/1‑69 (MQ=255)
                                        tACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCg                    <  1:91998/69‑1 (MQ=255)
                                        tACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCg                    <  1:3259595/69‑1 (MQ=255)
                                        tACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCg                    <  1:2237296/69‑1 (MQ=255)
                                        tACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCg                    <  1:1880823/69‑1 (MQ=255)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTGAGTTTTCACTTCCGGTGCCGCtt                                   >  1:2334961/1‑53 (MQ=39)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCttt                                  >  1:2654772/1‑54 (MQ=255)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCttt                                  >  1:2357483/1‑54 (MQ=255)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCttt                                  >  1:2924256/1‑54 (MQ=255)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCttt                                  >  1:3105130/1‑54 (MQ=255)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCttt                                  >  1:3250350/1‑54 (MQ=255)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCttt                                  >  1:1986469/1‑54 (MQ=255)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCttt                                  >  1:702193/1‑54 (MQ=255)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCttt                                  >  1:1957711/1‑54 (MQ=255)
                                         aCGCCACGGTGAGCAGGCACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCttt                                  >  1:1774074/1‑54 (MQ=255)
                                              aCGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCg                    <  1:2035578/63‑1 (MQ=255)
                                                          tACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGaa  >  1:2587009/1‑69 (MQ=255)
                                                          tACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGaa  >  1:2519945/1‑69 (MQ=255)
                                                          tACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGa   >  1:568316/1‑68 (MQ=255)
                                                          |                                                                    
CAGCCAAAAAGAATGATTGCCAGCGCCATCAACGCGTATTTACGCCACGGTGAGCAGGTACATTTTTCCTCAGTTTTCACTTCCGGTGCCGCTTTTTGTGCCGCCTGCGGCTGAGCTGATACCTGAA  >  minE/1004351‑1004477

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: