Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 253,007 C→T 50.6% L127L (TTG→TTA yajB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE253,0070CT50.6% ‑0.1 / 71.0 77L127L (TTG→TTAyajBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (22/16);  new base T (24/15);  total (46/31)
Fisher's exact test for biased strand distribution p-value = 8.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGG  >  minE/252941‑253074
                                                                  |                                                                   
cATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACcac                                                                     <  1:1871907/67‑1 (MQ=255)
  tGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGt                                                                 <  1:3337756/69‑1 (MQ=255)
    ccGTTTAACACGTTCTGGATGAAATCCATATCTCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:1133767/1‑69 (MQ=255)
    ccGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:2021513/1‑69 (MQ=255)
    ccGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:1415454/1‑69 (MQ=255)
    ccGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:1613007/1‑69 (MQ=255)
    ccGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:157739/1‑69 (MQ=255)
    ccGTATAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTa                                                                >  1:1344380/1‑68 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:1818970/68‑1 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:2672577/68‑1 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:2654531/68‑1 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:2817579/68‑1 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:2968124/68‑1 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:3197189/68‑1 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:379933/68‑1 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:1437563/68‑1 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:1195761/68‑1 (MQ=255)
            cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:615763/68‑1 (MQ=255)
                  cTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGa                                                  >  1:1597218/1‑68 (MQ=255)
                  cTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGa                                                  >  1:1664991/1‑68 (MQ=255)
                  cTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGa                                                  >  1:2818090/1‑68 (MQ=255)
                    cgATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATaa                                               >  1:2751746/2‑69 (MQ=255)
                        gAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATaaaatg                                           >  1:2472950/1‑67 (MQ=255)
                        gAAATCCATATCGCGATAGCGCACAAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACg                                           >  1:1799779/1‑69 (MQ=255)
                         aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:539004/1‑55 (MQ=255)
                         aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:2297069/1‑55 (MQ=255)
                         aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:1935865/1‑55 (MQ=255)
                         aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:2950221/1‑55 (MQ=255)
                         aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:2498018/1‑55 (MQ=255)
                         aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:1874934/1‑55 (MQ=255)
                         aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:305653/1‑55 (MQ=255)
                         aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:1100485/1‑55 (MQ=255)
                         aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:1595548/1‑55 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:1002401/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2768636/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2589113/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2853650/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2166070/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2855718/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:1254750/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2003731/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2917401/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2975837/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:1809219/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:3194728/69‑2 (MQ=255)
                                              accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:490115/69‑2 (MQ=255)
                                               ccAACCACTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2429672/68‑2 (MQ=255)
                                               ccAACCACTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:1763570/68‑2 (MQ=255)
                                                   ccACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGtcat                <  1:1264403/69‑1 (MQ=255)
                                                   ccACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGtcat                <  1:1528940/69‑1 (MQ=255)
                                                        gCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACtt           >  1:1833913/1‑69 (MQ=255)
                                                         ctctGACCACAAGTAATTGTTCAGATTGATAAAACGTGGc                                       <  1:2318949/40‑1 (MQ=255)
                                                         ctctGACCACAAGTAATTGTTCAGATTGATAAAACGTGGc                                       <  1:2079556/40‑1 (MQ=255)
                                                         ctctGACCACAAGTAATTGTTCAGATTGATAAAACGTGGc                                       <  1:1432341/40‑1 (MQ=255)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCg                                      >  1:2431560/1‑39 (MQ=18)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:88859/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:666062/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:1258107/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:545916/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:1615491/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:508192/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:136049/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:423/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:395572/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:924416/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:3295866/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:2727952/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:3189399/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:3180355/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:3172654/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:1901592/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:1908354/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:1916792/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:1931129/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:2203408/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:1000458/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:2588595/1‑63 (MQ=19)
                                                           ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatc               >  1:920606/1‑62 (MQ=19)
                                                                 acaAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCggg  <  1:3397077/69‑1 (MQ=255)
                                                                  |                                                                   
CATGCCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGG  >  minE/252941‑253074

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: