Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,530,159 C→T 100% G218G (GGG→GGA yfcA ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,530,1590CT77.2% 58.4 / 12.8 48G218G (GGG→GGAyfcAconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (2/9);  new base T (37/0);  total (39/9)
Fisher's exact test for biased strand distribution p-value = 3.28e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

ATCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCACTTTGCCGCCGAGAATAAACAGCAGTA  >  minE/1530100‑1530220
                                                           |                                                             
aTCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGc                                                      <  1:3179328/69‑1 (MQ=255)
 tCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGc                                                      >  1:1837521/1‑68 (MQ=255)
 tCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGc                                                      >  1:2041217/1‑68 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCCTGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:1695944/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCctacca                                                        >  1:3209488/1‑51 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:110387/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:444492/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:2714545/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:2611841/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:2566363/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:2424449/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:2278113/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:2230380/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:2193400/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:1935997/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:1164705/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:1764979/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:1384973/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:12825/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCCATGCGCGCGCCGAGGAACTGTCCTACCAACATCACAAATcctg                                         >  1:132774/1‑64 (MQ=255)
                tGCTCAACACCAGTCTTGACCCATGCGCGCGCCGAGGAACTGTCCGACCAACATCACAAATcctg                                         >  1:3256294/1‑63 (MQ=255)
                         ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAATCCCGTCGCCCAAATCa                             <  1:1342324/69‑1 (MQ=255)
                         ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa                             <  1:552087/69‑1 (MQ=255)
                         ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa                             <  1:361465/69‑1 (MQ=255)
                         ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa                             <  1:3414275/69‑1 (MQ=255)
                         ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa                             <  1:2464927/69‑1 (MQ=255)
                         ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa                             <  1:702636/69‑1 (MQ=255)
                         ccAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCa                             <  1:1931418/69‑1 (MQ=255)
                             tCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGc                                                      <  1:3230814/40‑1 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCCCCgaga               >  1:2883605/1‑56 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:3342856/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:995878/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:348363/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:39789/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:1045303/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:697596/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:883012/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:3289835/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:3226122/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:1162714/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:3055395/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:2884694/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:2638091/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:2500953/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:1732022/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:1896842/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:1905094/1‑67 (MQ=255)
                                                    ggAACTGTCCTACCAACATCACAAATCCTGTCGCCCAAATCACTTTGCCACCGAGAATAAAcagcagca  >  1:1926696/1‑67 (MQ=255)
                                                           |                                                             
ATCAGCTTTTGTCCTTTGCTCAACACCAGTCGTGAACCCATGCGCGCGCCGAGGAACTGCCCCACCAGCATGACAAAGCCCGTCGCCCAAATCACTTTGCCGCCGAGAATAAACAGCAGTA  >  minE/1530100‑1530220

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: