Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,021,226 G→A 40.2% V273I (GTT→ATT)  yqhG → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,021,2260GA40.2% 30.4 / 43.8 67V273I (GTT→ATT) yqhGconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (17/23);  new base A (21/6);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 5.72e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

CCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGC  >  minE/2021161‑2021288
                                                                 |                                                              
ccgATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGttt                                                             <  1:2282533/69‑1 (MQ=255)
 cgATTCCGCCCGCATTCATGTTCGCTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGAtttt                                                            <  1:870518/69‑1 (MQ=255)
 cgATTCCGCCCGCATTCATGTTCGCTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGAtttt                                                            <  1:662483/69‑1 (MQ=255)
 cgATTCCGCCCGCATTCATGTTCGCTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGAtttt                                                            <  1:2322578/69‑1 (MQ=255)
 cgATTCCGCCCGCATTCATGTTCGCTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGAtttt                                                            <  1:2389359/69‑1 (MQ=255)
 cgATTCCGCCCGCATTCATGTTCGCTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGAtttt                                                            <  1:1559069/69‑1 (MQ=255)
 cgATTCCGCCCGCATTCATGTTCGCTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGAtttt                                                            <  1:3006872/69‑1 (MQ=255)
                 cATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAg                                              >  1:385122/1‑67 (MQ=255)
                   tGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCa                                           <  1:2747974/68‑1 (MQ=255)
                   tGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCa                                           <  1:3194509/68‑1 (MQ=255)
                   tGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCa                                           <  1:1147782/68‑1 (MQ=255)
                   tGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCa                                           <  1:2218333/68‑1 (MQ=255)
                   tGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCa                                           <  1:2066427/68‑1 (MQ=255)
                   tGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAACAACTGGTTTTACAGCGGCGCGAAGGCa                                           <  1:2853518/68‑1 (MQ=255)
                     ttCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGt                                        >  1:798981/1‑69 (MQ=255)
                     ttCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGt                                        >  1:452349/1‑69 (MQ=255)
                     ttCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGt                                        >  1:3278398/1‑69 (MQ=255)
                     ttCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGt                                        >  1:1590493/1‑69 (MQ=255)
                     ttCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGt                                        >  1:2042425/1‑69 (MQ=255)
                                               aacGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:2144270/69‑1 (MQ=255)
                                               aacGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:1903251/69‑1 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGTCGCGAAGGTAAg                                         >  1:627198/1‑42 (MQ=37)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGc                              >  1:2143316/1‑53 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGc                              >  1:3134760/1‑53 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:686886/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:374932/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1188883/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1203318/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1254003/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:3077291/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1308938/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1008023/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:2820249/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:2757964/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:2678802/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:423368/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1770414/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:529859/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1710996/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1639990/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1587161/1‑54 (MQ=255)
                                               aacGAAGAGCAAGAACTGATTTTACAGCGGCGCGAAGGTAAGTGGGAAATCGCt                             >  1:1582095/1‑54 (MQ=255)
                                                acGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              >  1:3129057/1‑68 (MQ=255)
                                                acGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              >  1:3186705/1‑68 (MQ=255)
                                                acGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              >  1:319870/1‑68 (MQ=255)
                                                acGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              >  1:2564911/1‑68 (MQ=255)
                                                acGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              >  1:356520/1‑68 (MQ=255)
                                                acGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATACGCCCg              >  1:447652/1‑68 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGAttt                        >  1:2997223/1‑56 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGAttt                        >  1:1536322/1‑56 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGAttt                        >  1:1794792/1‑56 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGAtt                         >  1:2428068/1‑55 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:3395656/66‑1 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:332785/66‑1 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:310535/66‑1 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:2419787/66‑1 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:2113277/66‑1 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:200589/66‑1 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:1862315/66‑1 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:1517575/66‑1 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:1260445/66‑1 (MQ=255)
                                                  gAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCg              <  1:1133637/66‑1 (MQ=255)
                                                       gCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGc        >  1:918591/1‑67 (MQ=255)
                                                           gAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGc  <  1:989201/69‑1 (MQ=255)
                                                           gAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGc  <  1:874137/69‑1 (MQ=255)
                                                           gAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGc  <  1:1371912/69‑1 (MQ=255)
                                                             aCTTGTTTTACAGCGGCGCGAAGGCAAGTTGTAAATCGCTGATTTTATCCGCCCg              <  1:2500956/55‑1 (MQ=255)
                                                                 |                                                              
CCGATTCCGCCCGCATTCATGTTCGTTTTACGCTGACCAATGGCAACAACGAAGAGCAAGAACTGGTTTTACAGCGGCGCGAAGGCAAGTGGGAAATCGCTGATTTTATCCGCCCGAACAGCGGCAGC  >  minE/2021161‑2021288

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: