Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 253,013 A→G 54.2% N125N (AAT→AAC yajB ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE253,0130AG54.2% 3.6 / 58.1 72N125N (AAT→AACyajBconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (16/17);  new base G (24/15);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 3.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGG  >  minE/252945‑253074
                                                                    |                                                             
ccGTTTAACACGTTCTGGATGAAATCCATATCTCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:1133767/1‑69 (MQ=255)
ccGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:2021513/1‑69 (MQ=255)
ccGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:1613007/1‑69 (MQ=255)
ccGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:157739/1‑69 (MQ=255)
ccGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTaa                                                               >  1:1415454/1‑69 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:379933/68‑1 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:2672577/68‑1 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:1818970/68‑1 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:3197189/68‑1 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:1437563/68‑1 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:2968124/68‑1 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:2817579/68‑1 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:615763/68‑1 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:2654531/68‑1 (MQ=255)
        cacGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        <  1:1195761/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:305653/1‑55 (MQ=255)
                     aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:1595548/1‑55 (MQ=255)
                     aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:2498018/1‑55 (MQ=255)
                     aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:1874934/1‑55 (MQ=255)
                     aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:2950221/1‑55 (MQ=255)
                     aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:2297069/1‑55 (MQ=255)
                     aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:539004/1‑55 (MQ=255)
                     aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:1935865/1‑55 (MQ=255)
                     aaaTCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCa                                                        >  1:1100485/1‑55 (MQ=255)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2975837/69‑2 (MQ=255)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2768636/69‑2 (MQ=255)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2166070/69‑2 (MQ=255)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2589113/69‑2 (MQ=255)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2003731/69‑2 (MQ=255)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2853650/69‑2 (MQ=255)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2855718/69‑2 (MQ=255)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:2917401/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)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:1254750/69‑2 (MQ=255)
                                          accaaccaCTGCTCTGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAAta                     <  1:1002401/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:924416/1‑63 (MQ=19)
                                                       ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:508192/1‑63 (MQ=19)
                                                       ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:136049/1‑63 (MQ=19)
                                                       ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:2588595/1‑63 (MQ=19)
                                                       ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:423/1‑63 (MQ=19)
                                                       ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:395572/1‑63 (MQ=19)
                                                       ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:3295866/1‑63 (MQ=19)
                                                       ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatca              >  1:1615491/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:2727952/1‑63 (MQ=19)
                                                       ctGACCATAAGTAGTTGTTCAGATTGATAAAGCGTTGCGGTGAGTCTGGAAGAATAGtcatc               >  1:920606/1‑62 (MQ=19)
                                                             acaAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCggg  <  1:3397077/69‑1 (MQ=255)
                                                               aaGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCggg  <  1:1818901/67‑1 (MQ=255)
                                                                    |                                                             
CCGTTTAACACGTTCTGGATGAAATCCATATCGCGATAGCGCACCAGCCACTGCTCTGACCACAAGTAATTGTTCAGATTGATAAAACGTGGCGGTGAGTCCGGCAAAATCGTCATCACTTGCTCGCGGG  >  minE/252945‑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: