Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,464,528 C→T 37.4% S62S (AGC→AGT atpB → F0 sector of membrane‑bound ATP synthase, subunit a

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,464,5280CT37.4% 25.8 / 27.0 48S62S (AGC→AGTatpBF0 sector of membrane‑bound ATP synthase, subunit a
Reads supporting (aligned to +/- strand):  ref base C (17/13);  new base T (6/12);  total (23/25)
Fisher's exact test for biased strand distribution p-value = 1.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.91e-01

TCAATATTGACTCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGCGGTGTGCCAGGTAAGTTTCAGACCGCGATTGAGCTGGTGATCGGC  >  minE/2464464‑2464594
                                                                |                                                                  
tcaatATTGACTCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCcgtagcgta                                                                 >  1:814629/1‑68 (MQ=255)
        gACTCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCaaaaa                                                         <  1:2917140/68‑1 (MQ=255)
        gACTCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCaaaaa                                                         <  1:2549328/68‑1 (MQ=255)
        gACTCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCaaaaa                                                         <  1:2360918/68‑1 (MQ=255)
        gACTCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCaaaaa                                                         <  1:757247/68‑1 (MQ=255)
        gACTCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAg                                                        >  1:1864896/1‑69 (MQ=255)
           tCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCg                                                     <  1:2863718/69‑1 (MQ=255)
           tCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCg                                                     <  1:495661/69‑1 (MQ=255)
           tCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCg                                                     <  1:170026/69‑1 (MQ=255)
           tCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCg                                                     <  1:749206/69‑1 (MQ=255)
            ccATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCg                                                     >  1:10237/1‑68 (MQ=255)
            ccATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCg                                                     <  1:745138/68‑1 (MQ=255)
             cATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCg                                                     >  1:1888053/1‑67 (MQ=255)
             cATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCg                                                     >  1:1684890/1‑67 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAg                                                        >  1:1568094/1‑60 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:1776264/1‑69 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:1931486/1‑69 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:559901/1‑69 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:1609370/1‑69 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:490844/1‑69 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:22680/1‑69 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:1581254/1‑69 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:3363957/1‑69 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:1535536/1‑69 (MQ=255)
                 ttcttcTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:2798242/1‑69 (MQ=255)
                 ttattctCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGc                                               >  1:466801/4‑69 (MQ=255)
                         ggtggtGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGcc                                       >  1:1051871/1‑69 (MQ=255)
                         ggtggtGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGcc                                       >  1:2913635/1‑69 (MQ=255)
                         ggtggtGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGcc                                       >  1:912934/1‑69 (MQ=255)
                         ggtggtGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGcc                                       >  1:1236865/1‑69 (MQ=255)
                         ggtggtGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGcc                                       >  1:1951186/1‑69 (MQ=255)
                         ggtggtGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGCCGTGcc                                       >  1:2631172/1‑69 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:2307074/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:19419/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:689547/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:2039268/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:2251824/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:2821417/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:465947/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:417408/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:2380945/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:3308366/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:3137957/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGTGTAGCCAAAAAGGCAACCAGCGGCGTGCCAGGTAAGttt                             <  1:2586080/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGCGGTGTGCCAGGTAAGttt                             <  1:71188/67‑1 (MQ=255)
                                     tCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGCGGTGTGCCAGGTAAGttt                             <  1:1142123/67‑1 (MQ=255)
                                             tCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGCGGTGTGCCAGGTAAGTTTCAGACcgcg                    <  1:3401413/68‑1 (MQ=255)
                                             tCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGCGGTGTGCCAGGTAAGTTTCAGACcgcg                    <  1:1633615/68‑1 (MQ=255)
                                                     ttATTCCGTAGCGTAGCCAAAAAGGCGACCAGCGGTGTGCCAGGTAAGTTTCAGACCGCGATTGAGCTg           >  1:1884471/1‑69 (MQ=255)
                                                              agcgtagcCAAAAAGGCGACCAGCGGTGTGCCAGGTAAGTTTCAGACCGCGATTGAGCTGGTGATCGGc  <  1:402823/69‑1 (MQ=255)
                                                              agcgtagcCAAAAAGGCGACCAGCGGTGTGCCAGGTAAGTTTCAGACCGCGATTGAGCTGGTGATCGGc  <  1:2626123/69‑1 (MQ=255)
                                                                |                                                                  
TCAATATTGACTCCATGTTCTTCTCGGTGGTGCTGGGTCTGTTGTTCCTGGTTTTATTCCGTAGCGTAGCCAAAAAGGCGACCAGCGGTGTGCCAGGTAAGTTTCAGACCGCGATTGAGCTGGTGATCGGC  >  minE/2464464‑2464594

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: