Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,628,879 C→T 33.2% T109T (ACG→ACA yhhK ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,628,8790CT33.2% 31.3 / 19.0 42T109T (ACG→ACAyhhKconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (20/8);  new base T (7/7);  total (27/15)
Fisher's exact test for biased strand distribution p-value = 1.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTGACGATTAACACTTCTCCCAGCCGCCCTGTTGTGCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGATT  >  minE/2628815‑2628941
                                                                |                                                              
ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc                                                             >  1:2475894/1‑68 (MQ=255)
ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc                                                             >  1:2377731/1‑68 (MQ=255)
ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc                                                             >  1:2279889/1‑68 (MQ=255)
ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc                                                             >  1:2864037/1‑68 (MQ=255)
ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc                                                             >  1:803898/1‑68 (MQ=255)
ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc                                                             >  1:958429/1‑68 (MQ=255)
ttGACGATTAATGCTTCTCCCAGCCACCCTGTTGTGCCGTAAATCCCAGTGCTCGCATAAACGCTGtc                                                             >  1:1161108/1‑68 (MQ=255)
                                 gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATcctc                             <  1:1523523/67‑3 (MQ=255)
                                 gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc                             <  1:1422515/67‑1 (MQ=255)
                                 gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc                             <  1:1685311/67‑1 (MQ=255)
                                 gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc                             <  1:990470/67‑1 (MQ=255)
                                 gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc                             <  1:2993171/67‑1 (MQ=255)
                                 gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc                             <  1:870670/67‑1 (MQ=255)
                                 gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc                             <  1:898170/67‑1 (MQ=255)
                                 gtgCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGc                             <  1:3056215/67‑1 (MQ=255)
                                        aaaTCCCAGTGCTCGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca                    <  1:2727212/69‑1 (MQ=255)
                                        aaaTCCCAGTGCTCGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca                    <  1:2690202/69‑1 (MQ=255)
                                        aaaTCCCAGTGCTCGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca                    <  1:1895666/69‑1 (MQ=255)
                                        aaaTCCCAGTGCTCGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca                    <  1:3254474/69‑1 (MQ=255)
                                        aaaTCCCAGTGCTAGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca                    <  1:305723/69‑1 (MQ=255)
                                        aaaTCCCAGTGCTAGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca                    <  1:1030591/69‑1 (MQ=255)
                                        aaaTCCCAGTGCTAGCATAAACGCTGTCATCACGCCGCGATCTTCCACGCCAGCATCCGCCATccacca                    <  1:1157856/69‑1 (MQ=255)
                                                gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa            >  1:817649/1‑69 (MQ=255)
                                                gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa            >  1:2834944/1‑69 (MQ=255)
                                                gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa            >  1:2152505/1‑69 (MQ=255)
                                                gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa            >  1:2916705/1‑69 (MQ=255)
                                                gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa            >  1:1765769/1‑69 (MQ=255)
                                                gcgcCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGaaa            >  1:3174182/1‑69 (MQ=255)
                                                       cATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAgg     >  1:2747724/1‑69 (MQ=255)
                                                       cATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAgg     >  1:1853340/1‑69 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCa                           >  1:1396426/1‑43 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc                         >  1:99189/1‑45 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc                         >  1:246366/1‑45 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc                         >  1:2457516/1‑45 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc                         >  1:2103674/1‑45 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATc                         >  1:1203590/1‑45 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt  >  1:305355/1‑68 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt  >  1:679245/1‑68 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt  >  1:3001623/1‑68 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt  >  1:908715/1‑68 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt  >  1:1352616/1‑68 (MQ=255)
                                                           aaCGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGAtt  >  1:1051401/1‑68 (MQ=255)
                                                                |                                                              
TTGACGATTAACACTTCTCCCAGCCGCCCTGTTGTGCCGTAAACCCCAGCGCCTGCATAAACGCCGTCATCACACCGCGATCTTCCACGCCTGCATCCGCCATCCACCAGCATGAAACGCCAGGATT  >  minE/2628815‑2628941

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: