Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,974,491 T→C 29.2% S264S (TCT→TCC metK → methionine adenosyltransferase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,974,4910TC29.2% 71.1 / 32.9 75S264S (TCT→TCCmetKmethionine adenosyltransferase 1
Reads supporting (aligned to +/- strand):  ref base T (36/17);  new base C (13/9);  total (49/26)
Fisher's exact test for biased strand distribution p-value = 5.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATC  >  minE/1974430‑1974557
                                                             |                                                                  
gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                                  >  1:1804815/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:3215994/1‑64 (MQ=255)
gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                                  >  1:2060899/1‑64 (MQ=255)
gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                                  >  1:715876/1‑64 (MQ=255)
gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                                  >  1:2960879/1‑64 (MQ=255)
gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                                  >  1:848540/1‑64 (MQ=255)
gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                                  >  1:1427821/1‑64 (MQ=255)
gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                                  >  1:1354889/1‑64 (MQ=255)
gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                                  >  1:923293/1‑64 (MQ=255)
gTCGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCgg                                                                  >  1:896607/1‑64 (MQ=255)
  cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                                             <  1:320879/67‑1 (MQ=255)
  cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                                             <  1:2528568/67‑1 (MQ=255)
  cGTAAAATCATCGTTGATACCTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAg                                                             <  1:1653393/67‑1 (MQ=255)
   gTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATc                                                          >  1:1986681/1‑69 (MQ=255)
          tATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                                   >  1:1428852/1‑69 (MQ=255)
          tATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                                   <  1:3083886/69‑1 (MQ=255)
                tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                                   <  1:3046012/63‑1 (MQ=255)
                tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                                   <  1:3283788/63‑1 (MQ=255)
                tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                                   <  1:2132095/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)
                tGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCaaa                                                   <  1:1431015/63‑1 (MQ=255)
                     ccTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGtt                                        <  1:1824285/69‑1 (MQ=255)
                     ccTACGGCGGCATGGCGCGTCATGGTGGCGGTGCGTTCTCCGGTAAAGATCCATCAAAAGTGGACCGtt                                        <  1:2567405/69‑1 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:3190460/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:2945154/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:2817366/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:3229188/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:3139853/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:2630749/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:313293/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:42752/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:1276951/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:2145291/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:539605/1‑68 (MQ=255)
                       tACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTc                                       >  1:1664341/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:110733/1‑69 (MQ=255)
                              gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt                               >  1:8868/1‑69 (MQ=255)
                              gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt                               >  1:1378215/1‑69 (MQ=255)
                              gcATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCt                               >  1:3334707/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:2644193/1‑42 (MQ=255)
                                          aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTgg                                              >  1:121259/1‑42 (MQ=255)
                                          aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGtt                   >  1:2103713/1‑69 (MQ=255)
                                          aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt                    <  1:885863/68‑1 (MQ=255)
                                          aCGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGt                    <  1:2078837/68‑1 (MQ=255)
                                                gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc              >  1:2611672/1‑68 (MQ=255)
                                                gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc              >  1:1088949/1‑68 (MQ=255)
                                                gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc              >  1:2502417/1‑68 (MQ=255)
                                                gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc              >  1:1296277/1‑68 (MQ=255)
                                                gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc              >  1:152024/1‑68 (MQ=255)
                                                gCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTc              >  1:81371/1‑68 (MQ=255)
                                                     gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa        <  1:2982599/69‑1 (MQ=255)
                                                     gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa        <  1:791583/69‑1 (MQ=255)
                                                     gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa        <  1:1890090/69‑1 (MQ=255)
                                                     gCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGaaa        <  1:662400/69‑1 (MQ=255)
                                                          ctctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACAt   >  1:2039310/1‑69 (MQ=255)
                                                          ctctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACAt   >  1:879424/1‑69 (MQ=255)
                                                          ctctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACAt   >  1:884643/1‑69 (MQ=255)
                                                          ctctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACAt   >  1:414008/1‑69 (MQ=255)
                                                          ctctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACAt   >  1:3215614/1‑69 (MQ=255)
                                                            ctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATc  <  1:1867269/68‑1 (MQ=255)
                                                            ctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATc  <  1:834509/68‑1 (MQ=255)
                                                            ctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATc  <  1:2395163/68‑1 (MQ=255)
                                                            ctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATc  <  1:131346/68‑1 (MQ=255)
                                                            ctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATc  <  1:2927626/68‑1 (MQ=255)
                                                            ctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATc  <  1:1270500/68‑1 (MQ=255)
                                                            ctGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATc  <  1:964211/68‑1 (MQ=255)
                                                             tGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATc  <  1:716367/67‑1 (MQ=255)
                                                             |                                                                  
GTCGTAAAATTATCGTTGATACCTACGGCGGCATGGCGCGTCACGGTGGCGGTGCATTCTCTGGTAAAGATCCATCAAAAGTGGACCGTTCCGCAGCCTACGCAGCACGTTATGTCGCGAAAAACATC  >  minE/1974430‑1974557

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: