Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,974,485 A→G 30.2% A262A (GCA→GCG metK → methionine adenosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,974,4850AG30.2% 67.3 / 35.2 76A262A (GCA→GCGmetKmethionine adenosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base A (36/17);  new base G (14/9);  total (50/26)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.02e-01

TGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAA  >  minE/1974417‑1974551
                                                                    |                                                                  
tGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCa                                                                    <  1:935139/69‑1 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTtctccg                                                             >  1:2060898/1‑61 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:3215994/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:2960879/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:3418355/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:2548129/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:628568/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:715876/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:2060899/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:848540/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:1804815/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:1427821/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:923293/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:896607/1‑64 (MQ=255)
             gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                            >  1:1354889/1‑64 (MQ=255)
               cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                                       <  1:320879/67‑1 (MQ=255)
               cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                                       <  1:1653393/67‑1 (MQ=255)
               cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                                       <  1:2528568/67‑1 (MQ=255)
                gTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATc                                                    >  1:1986681/1‑69 (MQ=255)
                       tATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                             <  1:3083886/69‑1 (MQ=255)
                       tATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                             >  1:1428852/1‑69 (MQ=255)
                             tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                             <  1:854794/63‑1 (MQ=255)
                             tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                             <  1:3046012/63‑1 (MQ=255)
                             tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                             <  1:3283788/63‑1 (MQ=255)
                             tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                             <  1:1418921/63‑1 (MQ=255)
                             tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                             <  1:1431015/63‑1 (MQ=255)
                             tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                             <  1:2132095/63‑1 (MQ=255)
                             tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                             <  1:718575/63‑1 (MQ=255)
                                  ccTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGtt                                  <  1:2567405/69‑1 (MQ=255)
                                  ccTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGtt                                  <  1:1824285/69‑1 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:539605/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:3139853/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:313293/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:1276951/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:42752/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:1664341/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:2630749/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:3190460/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:3229188/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:2817366/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:2945154/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                 >  1:2145291/1‑68 (MQ=255)
                                    tACGGCGGCATGGCGCGTCACGGAGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGAc                                      >  1:2172804/1‑63 (MQ=255)
                                        gcggcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAg                            >  1:1214761/1‑69 (MQ=255)
                                        gcggcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAg                            >  1:734721/1‑69 (MQ=255)
                                         cggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGc                           <  1:53537/69‑1 (MQ=255)
                                         cggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGc                           <  1:654477/69‑1 (MQ=255)
                                         cggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGc                           <  1:222714/69‑1 (MQ=255)
                                          ggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                                       <  1:2657417/40‑1 (MQ=25)
                                           gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt                         >  1:8868/1‑69 (MQ=255)
                                           gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt                         >  1:3334707/1‑69 (MQ=255)
                                           gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt                         >  1:40918/1‑69 (MQ=255)
                                           gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt                         >  1:110733/1‑69 (MQ=255)
                                           gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt                         >  1:1378215/1‑69 (MQ=255)
                                               ggCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgc                     <  1:2430509/69‑1 (MQ=255)
                                               ggCGCGTCACGGTGGCGGTGCACTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACgc                     <  1:3292727/69‑1 (MQ=255)
                                                      cACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt              <  1:2756557/69‑1 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                                        >  1:3069050/1‑42 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                                        >  1:2103211/1‑42 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                                        >  1:2233488/1‑42 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                                        >  1:242123/1‑42 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                                        >  1:2582660/1‑42 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                                        >  1:121259/1‑42 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                                        >  1:2644193/1‑42 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGtt             >  1:2103713/1‑69 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt              <  1:2078837/68‑1 (MQ=255)
                                                       aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt              <  1:885863/68‑1 (MQ=255)
                                                             gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc        >  1:1296277/1‑68 (MQ=255)
                                                             gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc        >  1:1088949/1‑68 (MQ=255)
                                                             gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc        >  1:152024/1‑68 (MQ=255)
                                                             gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc        >  1:81371/1‑68 (MQ=255)
                                                             gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc        >  1:2502417/1‑68 (MQ=255)
                                                             gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc        >  1:2611672/1‑68 (MQ=255)
                                                                  gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa  <  1:1890090/69‑1 (MQ=255)
                                                                  gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa  <  1:791583/69‑1 (MQ=255)
                                                                  gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa  <  1:662400/69‑1 (MQ=255)
                                                                  gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa  <  1:2982599/69‑1 (MQ=255)
                                                                    |                                                                  
TGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAA  >  minE/1974417‑1974551

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: