Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 429,936 A→G 38.0% S181S (TCA→TCG glnS → glutamyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE429,9360AG38.0% 33.2 / 45.3 71S181S (TCA→TCGglnSglutamyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (22/22);  new base G (25/2);  total (47/24)
Fisher's exact test for biased strand distribution p-value = 2.18e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

AAAATGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGCA  >  minE/429871‑430003
                                                                 |                                                                   
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2994747/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2266410/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2261406/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2299093/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2403011/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2518431/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2544054/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:258885/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2701220/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2775794/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2783343/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:294460/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:2961367/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:3108354/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:1734548/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:1623218/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:33587/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:1562050/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:3397771/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:3405021/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:389397/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:773522/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:795863/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:866878/1‑69 (MQ=255)
aaaaTGCGTACCGGCGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCg                                                                  >  1:1087211/1‑69 (MQ=255)
    tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:1109641/1‑69 (MQ=255)
    tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:3143217/1‑69 (MQ=255)
    tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:2612619/1‑69 (MQ=255)
    tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:702240/1‑69 (MQ=255)
    tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:2438195/1‑69 (MQ=255)
    tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:1981752/1‑69 (MQ=255)
    tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:1565534/1‑69 (MQ=255)
    tGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTa                                                              >  1:1786528/1‑69 (MQ=255)
        tGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGt                                                          <  1:1959853/69‑1 (MQ=255)
               ggtTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATgcgc                                                   >  1:642476/1‑69 (MQ=255)
                           ggTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCGTTTATCGTGATGCGCGATCCGGTGCTg                                       <  1:768106/69‑1 (MQ=255)
                           ggTAAAGCCTGCCTGCGTGCGAAGATCGACATGGCTTCGCCGTTTATCGTGATGCGCGATCCGGTGCTg                                       <  1:3114090/69‑1 (MQ=255)
                             tAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTa                                     >  1:2369693/1‑69 (MQ=255)
                             tAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTa                                     >  1:383954/1‑69 (MQ=255)
                             tAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTa                                     >  1:165257/1‑69 (MQ=255)
                             tAAAGCATGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTa                                     >  1:52923/1‑69 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:2259646/41‑1 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:872574/41‑1 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:1431230/41‑1 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:760422/41‑1 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:1528355/41‑1 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:2389804/41‑1 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:2053594/41‑1 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:2732282/1‑69 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:1480220/1‑69 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:386252/1‑69 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:2792997/1‑69 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:380953/1‑69 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:2206456/1‑69 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:1248436/1‑69 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:946653/1‑69 (MQ=255)
                              aaaGCCTGCCTGCGTGCGAAAAGCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTAc                                    >  1:3142454/1‑69 (MQ=255)
                               aaGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGtt                                                                <  1:579818/40‑1 (MQ=255)
                                            tgcgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTg                       <  1:1467800/68‑1 (MQ=255)
                                            tgcgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTg                       <  1:1763067/68‑1 (MQ=255)
                                            tgcgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTg                       <  1:195312/68‑1 (MQ=255)
                                            tgcgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGAGATCCGGTGCTGTACCGTATTAAATTTg                       <  1:2134135/68‑1 (MQ=255)
                                              cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:261337/69‑1 (MQ=255)
                                              cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1332111/69‑1 (MQ=255)
                                              cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1553901/69‑1 (MQ=255)
                                              cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1883385/69‑1 (MQ=255)
                                              cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1870957/69‑1 (MQ=255)
                                              cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:3093585/69‑1 (MQ=255)
                                              cgAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:1825980/69‑1 (MQ=255)
                                               gAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTg                    <  1:634950/68‑1 (MQ=255)
                                                        cATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACAccacca          <  1:570933/69‑1 (MQ=255)
                                                                 aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:2672017/1‑68 (MQ=255)
                                                                 aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:3193641/1‑68 (MQ=255)
                                                                 aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:1786949/1‑68 (MQ=255)
                                                                 aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:1806644/1‑68 (MQ=255)
                                                                 aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:227283/1‑68 (MQ=255)
                                                                 aCCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGca  >  1:2917844/1‑68 (MQ=255)
                                                                 |                                                                   
AAAATGCGTGCCGGTGGTTTTGAAGAAGGTAAAGCCTGCCTGCGTGCGAAAATCGACATGGCTTCACCGTTTATCGTGATGCGCGATCCGGTGCTGTACCGTATTAAATTTGCTGAACACCACCAGACTGGCA  >  minE/429871‑430003

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: