Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,475,894 T→C 38.6% G342G (GGT→GGC pstS → phosphate transporter subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,475,8940TC38.6% 17.5 / 26.6 44G342G (GGT→GGCpstSphosphate transporter subunit
Reads supporting (aligned to +/- strand):  ref base T (14/13);  new base C (13/4);  total (27/17)
Fisher's exact test for biased strand distribution p-value = 1.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

TGTAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCATCCGG  >  minE/2475840‑2475950
                                                      |                                                        
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTaat                                             >  1:1847705/2‑66 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:1070836/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:1048603/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:903891/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:769023/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:497106/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:2846159/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:2657719/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:1645933/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:2640731/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:2576746/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:2443059/2‑68 (MQ=255)
cgtAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACt                                             >  1:2327961/2‑68 (MQ=255)
          cAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACTAATAAAACTcc                                  <  1:982252/69‑1 (MQ=255)
          cAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACTAATAAAACTcc                                  <  1:2411608/69‑1 (MQ=255)
          cAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACTAATAAAACTcc                                  <  1:2091154/69‑1 (MQ=255)
           aGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGCAAGCCGCTGTACTAATAAAACTcc                                  <  1:570448/68‑1 (MQ=255)
                cgcgCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATaaaa                                      <  1:2098308/59‑1 (MQ=255)
                cgcgCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATaaaa                                      <  1:869780/59‑1 (MQ=255)
                cgcgCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATaaaa                                      <  1:2997036/59‑1 (MQ=255)
                cgcgCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATaaaa                                      <  1:297575/59‑1 (MQ=255)
                cgcgCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATaaaa                                      <  1:2078442/59‑1 (MQ=255)
                cgcgCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATaaaa                                      <  1:2346959/59‑1 (MQ=255)
                           gAAGACCAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGgtgt                  >  1:946749/1‑68 (MQ=255)
                                ccAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACt                                    >  1:2039805/1‑45 (MQ=255)
                                ccAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTAcg            >  1:185087/1‑69 (MQ=255)
                                ccAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTAcg            >  1:1473846/1‑69 (MQ=255)
                                ccAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTAcg            >  1:1353648/1‑69 (MQ=255)
                                ccAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTAcg            >  1:727033/1‑69 (MQ=255)
                                ccAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTAc             >  1:2127817/1‑68 (MQ=255)
                                     atTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCa       <  1:1575876/69‑1 (MQ=255)
                                     atTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCa       <  1:520589/69‑1 (MQ=255)
                                      ttAAAGACAGTAGCGGTAAGCCGCTGTACtaataa                                        >  1:1276337/1‑35 (MQ=255)
                                      ttAAAGACAGTAGCGGTAAGCCGCTGTACtaataa                                        >  1:1429555/1‑35 (MQ=255)
                                      ttAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTATGCCGCAt      >  1:859472/1‑69 (MQ=255)
                                      ttAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCAt      >  1:1267950/1‑69 (MQ=255)
                                      ttAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCAt      >  1:795569/1‑69 (MQ=255)
                                      ttAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCAt      >  1:2634133/1‑69 (MQ=255)
                                           gACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCATCCgg  <  1:2379422/68‑1 (MQ=255)
                                           gACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCATCCgg  <  1:846885/68‑1 (MQ=255)
                                           gACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCATCCgg  <  1:24330/68‑1 (MQ=255)
                                           gACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCATCCgg  <  1:2701050/68‑1 (MQ=255)
                                           gACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCATCCgg  <  1:2546514/68‑1 (MQ=255)
                                                 aGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCATCCgg  >  1:521334/1‑62 (MQ=255)
                                                      tAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCATCCgg  <  1:2508910/57‑1 (MQ=255)
                                                      |                                                        
TGTAGTTGAACAGGTTCGCGCTGCGTGGAAGACCAATATTAAAGACAGTAGCGGTAAGCCGCTGTACTAATAAAACTCCAGGCCGGGTACGGTGTTTTACGCCGCATCCGG  >  minE/2475840‑2475950

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: