Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 450,085 G→A 63.3% L85L (CTG→TTG)  kdpB ← potassium translocating ATPase, subunit B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE450,0850GA63.3% 36.9 / 41.0 68L85L (CTG→TTG) kdpBpotassium translocating ATPase, subunit B
Reads supporting (aligned to +/- strand):  ref base G (19/6);  new base A (32/11);  total (51/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

TTTTTTCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCCACAGCCAACCGCTAATGGCCGCGCTAAAC  >  minE/450032‑450149
                                                     |                                                                
ttttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAAtt                                                   >  1:3273872/1‑69 (MQ=255)
ttttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAAtt                                                   >  1:286500/1‑69 (MQ=255)
ttttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAAtt                                                   >  1:1813802/1‑69 (MQ=255)
ttttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAAtt                                                   >  1:475937/1‑69 (MQ=255)
ttttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAAtt                                                   >  1:2127498/1‑69 (MQ=255)
ttttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAAtt                                                   >  1:2044645/1‑69 (MQ=255)
ttttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAAt                                                    >  1:1611762/1‑68 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTa                                                  >  1:336759/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTa                                                  >  1:1741609/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTa                                                  >  1:3131022/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTa                                                  >  1:1589071/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTa                                                  >  1:3014997/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTa                                                  >  1:2162203/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTa                                                  >  1:1315117/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTa                                                  >  1:2339796/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAACTCCTCGGCGAAATTa                                                  >  1:182785/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAAtt                                                   >  1:1498558/1‑68 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAAtt                                                   >  1:2877169/1‑68 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTa                                                  >  1:2881082/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTa                                                  >  1:336315/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTa                                                  >  1:2241515/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTa                                                  >  1:1466450/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTa                                                  >  1:290223/1‑69 (MQ=255)
 tttttCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTa                                                  >  1:833313/1‑69 (MQ=255)
                   cTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATc                                <  1:2832965/69‑1 (MQ=255)
                   cTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATc                                <  1:3330041/69‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCcaca                            >  1:1757980/1‑68 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCcaca                            >  1:2608918/1‑68 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCcaca                            >  1:1605902/1‑68 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCcaca                            >  1:755042/1‑68 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:3266560/67‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:2536967/67‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:2989297/67‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:958777/67‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:1688204/67‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:1681538/67‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:1668986/67‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:1467563/67‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:1189977/67‑1 (MQ=255)
                        ggCCTGCGCTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:895720/67‑1 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:392618/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:44110/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:50417/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:61459/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:721457/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:2437675/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:2793809/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:1091420/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:1452982/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:1725050/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:1872685/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:2441991/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTa                                        >  1:1001291/1‑48 (MQ=255)
                                cTTTACTGCGGCCTTCTGCCAACGCCTAGGCGAAATTAGCGAACAGTa                                        >  1:557193/1‑48 (MQ=39)
                                 tttACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCCACAGCCAACCGCt                  <  1:256855/69‑1 (MQ=255)
                                   tACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGa                                  >  1:154130/1‑51 (MQ=255)
                                   tACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGa                                  >  1:1080369/1‑51 (MQ=255)
                                   tACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGa                                  >  1:2808720/1‑51 (MQ=255)
                                   tACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGa                                  >  1:1141914/1‑51 (MQ=255)
                                   tACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGa                                  >  1:1255917/1‑51 (MQ=255)
                                   tACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGa                                  >  1:2633415/1‑51 (MQ=255)
                                   tACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGa                                  >  1:1463041/1‑51 (MQ=255)
                                   tACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGa                                  >  1:1627929/1‑51 (MQ=255)
                                   tACTGCGGCCTTCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGa                                  >  1:1785705/1‑51 (MQ=255)
                                             ttCTGCCAACGCCTCGGCGAAATTAGCGAACAGTACGGTGACCcac                             <  1:2222698/46‑1 (MQ=255)
                                                    aGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCCACAGCCAACCGCTAATGGCCGCGCTAAAc  <  1:3417381/66‑1 (MQ=255)
                                                    aGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCCACAGCCAACCGCTAATGGCCGCGCTAAAc  <  1:2295532/66‑1 (MQ=255)
                                                    aGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCCACAGCCAACCGCTAATGGACGCGCTAAAc  <  1:757039/66‑1 (MQ=255)
                                                     |                                                                
TTTTTTCACCCCTTTCAGACTGTTGGCCTGCGCTTTACTGCGGCCTTCTGCCAGCGCCTCGGCGAAATTAGCGAACAGTACGGTGATCCACAGCCAACCGCTAATGGCCGCGCTAAAC  >  minE/450032‑450149

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: