Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,683,429 A→G 50.0% G515G (GGT→GGC yfhM ← conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,683,4290AG50.0% ‑5.1 / 32.2 38G515G (GGT→GGCyfhMconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (2/17);  new base G (2/17);  total (4/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGACACAGCTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTT  >  minE/1683362‑1683480
                                                                   |                                                   
cTGACACAGCTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCAc                                                    <  1:1942713/69‑1 (MQ=255)
cTGACACAGCTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCAc                                                    <  1:1398979/69‑1 (MQ=255)
           tCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGgtag                                         <  1:1242996/69‑1 (MQ=255)
           tCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGgtag                                         <  1:2013370/69‑1 (MQ=255)
           tCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGgtag                                         <  1:1729372/69‑1 (MQ=255)
           tCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGgtag                                         <  1:1408902/69‑1 (MQ=255)
              cGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGgtagt                                        <  1:2949719/67‑1 (MQ=255)
                   tGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGGTAGTACCccac                                 <  1:789096/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:673280/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:414778/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:791682/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:2556078/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:403102/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:3122028/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:2901494/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:2880732/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:2745676/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:1065061/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:1323820/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:1437284/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:2267665/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:1514926/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:2069768/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:2152414/69‑1 (MQ=255)
                    gggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCGCCATACAGGTAGTACCccacc                                <  1:2070866/69‑1 (MQ=255)
                     ggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGgtag                                         >  1:215683/1‑59 (MQ=255)
                     ggCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGGTAGTACCccacca                               >  1:290144/1‑69 (MQ=255)
                                                 aTTACCATTAGCAGGTGCGCCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGc     >  1:2669488/1‑67 (MQ=255)
                                                 aTTACCATTAGCAGGTGCGCCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGc     >  1:3016613/1‑67 (MQ=255)
                                                  ttACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGtt  <  1:89274/69‑1 (MQ=255)
                                                  ttACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGtt  <  1:2420504/69‑1 (MQ=255)
                                                  ttACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGtt  <  1:1545185/69‑1 (MQ=255)
                                                  ttACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGtt  <  1:2217765/69‑1 (MQ=255)
                                                  ttACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGtt  <  1:2851173/69‑1 (MQ=255)
                                                  ttACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGtt  <  1:2797041/69‑1 (MQ=255)
                                                  ttACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGtt  <  1:2745461/69‑1 (MQ=255)
                                                  ttACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGtt  <  1:2621814/69‑1 (MQ=255)
                                                  ttACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGtt  <  1:2401454/69‑1 (MQ=255)
                                                                   |                                                   
CTGACACAGCTTCACGCAGTGGGCGCAGGAAAAGTTGCCCTTGCAAAGTATTACCATTAGCAGGTGCACCATACAGGTAGTACCCCACCACGGAGAATTTCACTTCATCTTTCGGCGTT  >  minE/1683362‑1683480

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: