Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,974,473 C→T 46.1% H258H (CAC→CAT metK → methionine adenosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,974,4730CT46.1% 14.9 / 74.6 91H258H (CAC→CATmetKmethionine adenosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base C (30/19);  new base T (33/9);  total (63/28)
Fisher's exact test for biased strand distribution p-value = 1.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.15e-02

TGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTT  >  minE/1974413‑1974540
                                                            |                                                                   
tGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCgg                                                             <  1:2286095/69‑1 (MQ=255)
tGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCgg                                                             <  1:689477/69‑1 (MQ=255)
tGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCgg                                                             <  1:676020/69‑1 (MQ=255)
tGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCgg                                                             <  1:1988634/69‑1 (MQ=255)
tGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCgg                                                             <  1:515046/69‑1 (MQ=255)
  aCTGCGTTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:3084040/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGt                                                            >  1:2602986/1‑68 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:1547078/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:2391299/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:2566857/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:3157835/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:3354494/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:172884/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:3427384/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:368462/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:1601107/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:402622/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:506708/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:54751/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:592234/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:723831/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:838176/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTg                                                           >  1:1077384/1‑69 (MQ=255)
  aCTGCGGTCTGACTGGTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTAATGGTGGCGGTg                                                           >  1:1061804/1‑69 (MQ=39)
    tGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCa                                                         <  1:935139/69‑1 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTtctccg                                                  >  1:2060898/1‑61 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:2548129/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:2960879/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:1427821/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:628568/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:3418355/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:2060899/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:1354889/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:1804815/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:715876/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:848540/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:3215994/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:896607/1‑64 (MQ=255)
                 gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                 >  1:923293/1‑64 (MQ=255)
                   cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                            <  1:2528568/67‑1 (MQ=255)
                   cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                            <  1:320879/67‑1 (MQ=255)
                   cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                            <  1:1653393/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:2132095/63‑1 (MQ=255)
                                 tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                  <  1:3046012/63‑1 (MQ=255)
                                 tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                  <  1:3283788/63‑1 (MQ=255)
                                 tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                  <  1:1431015/63‑1 (MQ=255)
                                 tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                  <  1:1418921/63‑1 (MQ=255)
                                 tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                  <  1:718575/63‑1 (MQ=255)
                                 tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                  <  1:854794/63‑1 (MQ=255)
                                      ccTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGtt                       <  1:2567405/69‑1 (MQ=255)
                                      ccTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGtt                       <  1:1824285/69‑1 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:42752/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:539605/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:1664341/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:1276951/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:3229188/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:2630749/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:3190460/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:3139853/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:313293/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:2145291/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:2945154/1‑68 (MQ=255)
                                        tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                      >  1:2817366/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:222714/69‑1 (MQ=255)
                                             cggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGc                <  1:654477/69‑1 (MQ=255)
                                              ggcATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                            <  1:2657417/40‑1 (MQ=25)
                                               gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt              >  1:40918/1‑69 (MQ=255)
                                               gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt              >  1:8868/1‑69 (MQ=255)
                                               gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt              >  1:3334707/1‑69 (MQ=255)
                                               gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt              >  1:1378215/1‑69 (MQ=255)
                                               gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt              >  1:110733/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:2103211/1‑42 (MQ=255)
                                                           aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                             >  1:242123/1‑42 (MQ=255)
                                                           aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                             >  1:121259/1‑42 (MQ=255)
                                                           aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                             >  1:2233488/1‑42 (MQ=255)
                                                           aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                             >  1:2582660/1‑42 (MQ=255)
                                                           aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                             >  1:2644193/1‑42 (MQ=255)
                                                           aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                             >  1:3069050/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)
                                                            |                                                                   
TGACTGCGGTCTGACTGGTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTT  >  minE/1974413‑1974540

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: