Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 287,566 C→T 40.5% T70T (ACC→ACT lon → DNA‑binding ATP‑dependent protease La

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE287,5660CT40.5% 29.8 / 46.1 69T70T (ACC→ACTlonDNA‑binding ATP‑dependent protease La
Reads supporting (aligned to +/- strand):  ref base C (25/16);  new base T (8/20);  total (33/36)
Fisher's exact test for biased strand distribution p-value = 1.36e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.82e-01

GGTCGCGCAGAAAGAAGCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACGGCACCGT  >  minE/287512‑287619
                                                      |                                                     
ggTCGCGCAGAAAGAAGCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGc                                         <  1:668775/69‑1 (MQ=255)
      gcAGAAAGAAGCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCta                                    <  1:3159249/68‑1 (MQ=255)
            agaagCTTCAACGGATGAGCCGGGTTTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:2305033/69‑1 (MQ=255)
            agaagCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:106479/69‑1 (MQ=255)
            agaagCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:1202033/69‑1 (MQ=255)
            agaagCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:292008/69‑1 (MQ=255)
            agaagCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:2455574/69‑1 (MQ=255)
            agaagCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:2344283/69‑1 (MQ=255)
            agaagCTTCAACGGATGAGCCGGGTGTAAACGATCTTATCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:1291892/69‑1 (MQ=255)
            agaagCTTCAACGGATGAGCCGGGAGTAAACGATCTTATCACCGTCGGGACCGTGGCCTCTATATAGCa                             <  1:1959057/69‑1 (MQ=255)
             gaagCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:1406592/68‑1 (MQ=255)
             gaagCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:1643107/68‑1 (MQ=255)
                  ttCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCa                             <  1:1790506/63‑1 (MQ=255)
                    cAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCta                                    <  1:1642262/54‑1 (MQ=255)
                      aCGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     <  1:472338/67‑1 (MQ=255)
                      aCGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     <  1:2697751/67‑1 (MQ=255)
                         gATGAGCCGGGTGTAAACGATCTTTTCCCCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGCAACTg                >  1:1841466/1‑69 (MQ=255)
                         gATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTg                >  1:855466/1‑69 (MQ=255)
                         gATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTg                >  1:188765/1‑69 (MQ=255)
                         gATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTg                >  1:1147803/1‑69 (MQ=255)
                         gATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTg                >  1:3146991/1‑69 (MQ=255)
                         gATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTg                >  1:2562886/1‑69 (MQ=255)
                         gATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTg                >  1:1492528/1‑69 (MQ=255)
                         gATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTg                >  1:2327403/1‑69 (MQ=255)
                         gATGAGCCGGGTGTAAACGATCTTTTCACCGACGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTg                >  1:2484347/1‑69 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:1406343/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:1243300/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:2941702/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:2974740/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:2601473/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:374603/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:3141711/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:1500027/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:3190697/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:1846163/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:3233498/69‑1 (MQ=255)
                            gAGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:3404516/69‑1 (MQ=255)
                             aGCCGGGTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCt             <  1:875079/68‑1 (MQ=255)
                                 gggTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACgg        <  1:309169/69‑1 (MQ=255)
                                 gggTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACgg        <  1:2641817/69‑1 (MQ=255)
                                 gggTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACgg        <  1:3063057/69‑1 (MQ=255)
                                 gggTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACgg        <  1:3060446/69‑1 (MQ=255)
                                 gggTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACgg        <  1:2222268/69‑1 (MQ=255)
                                 gggTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACgg        <  1:1977697/69‑1 (MQ=255)
                                 gggTGTAAACGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACgg        <  1:1855848/69‑1 (MQ=255)
                                       aaaCGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTg                      >  1:3390828/1‑49 (MQ=255)
                                       aaaCGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:1750543/1‑50 (MQ=255)
                                       aaaCGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:754892/1‑50 (MQ=255)
                                       aaaCGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:2967874/1‑50 (MQ=255)
                                       aaaCGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:3050556/1‑50 (MQ=255)
                                       aaaCGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:2627744/1‑50 (MQ=255)
                                       aaaCGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:2714490/1‑50 (MQ=255)
                                       aaaCGATCTTTTCACTGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:2045295/1‑50 (MQ=255)
                                        aaCGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGaaa                   >  1:1116416/1‑51 (MQ=255)
                                        aaCGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACGGCACCGt  >  1:1259294/1‑68 (MQ=255)
                                        aaCGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACGGCACCGt  >  1:1688633/1‑68 (MQ=255)
                                        aaCGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACGGCACCGt  >  1:2538374/1‑68 (MQ=255)
                                        aaCGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACGGCACCGt  >  1:3040662/1‑68 (MQ=255)
                                        aaCGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACGGCACCGt  >  1:984634/1‑68 (MQ=255)
                                        aaCGATCATTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACGGCACCGt  >  1:1255165/1‑68 (MQ=255)
                                                 ttCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:984111/1‑40 (MQ=255)
                                                 ttCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:993208/1‑40 (MQ=255)
                                                 ttCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:328806/1‑40 (MQ=255)
                                                 ttCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:2952139/1‑40 (MQ=255)
                                                 ttCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:2778375/1‑40 (MQ=255)
                                                 ttCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:2600531/1‑40 (MQ=255)
                                                 ttCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:2518687/1‑40 (MQ=255)
                                                 ttCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:1726042/1‑40 (MQ=255)
                                                 ttCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGa                     >  1:1106609/1‑40 (MQ=255)
                                                      |                                                     
GGTCGCGCAGAAAGAAGCTTCAACGGATGAGCCGGGTGTAAACGATCTTTTCACCGTCGGGACCGTGGCCTCTATATTGCAGATGCTGAAACTGCCTGACGGCACCGT  >  minE/287512‑287619

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: