Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,530,282 A→G 100% F177F (TTT→TTC yfcA ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,530,2820AG82.0% 67.3 / 6.1 44F177F (TTT→TTCyfcAconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (7/1);  new base G (21/15);  total (28/16)
Fisher's exact test for biased strand distribution p-value = 2.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

AACCGCCGATGTTTGACGTTGCGTTGAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCAT  >  minE/1530221‑1530349
                                                             |                                                                   
aaCCGCCGATGTTTGACGTTGCGTTGAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:1326486/1‑64 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:2858379/1‑39 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:923426/1‑39 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:1899623/1‑39 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:640305/1‑39 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:1194066/1‑39 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:1210375/1‑39 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:1233595/1‑39 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:2211399/1‑39 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTaaa                                                                   >  1:604639/1‑39 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGt                                      >  1:1323100/1‑68 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGt                                      >  1:1906548/1‑68 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGt                                      >  1:2039808/1‑68 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGt                                      >  1:524043/1‑68 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGt                                      >  1:235905/1‑68 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGt                                      >  1:646284/1‑68 (MQ=255)
                         gAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGt                                      >  1:779702/1‑68 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAg                                                >  1:1282052/1‑50 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2974349/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2206810/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2667643/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2278977/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2323040/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2716927/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2667803/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2541690/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:3023648/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2087820/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:458461/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:2006830/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:501132/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:1888568/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:1553391/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:1431160/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:1354266/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:1075562/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGat                              >  1:109437/1‑67 (MQ=255)
                                 tGGCGTGAGCAGTGGCTTTGGCGAGGTTGAATCAGCACAGCGTAACGAAGGCCAGCGCGTAAAACGatc                             >  1:734135/1‑67 (MQ=255)
                                        aGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGc                      <  1:2426915/69‑1 (MQ=255)
                                        aGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGc                      <  1:1113643/69‑1 (MQ=255)
                                        aGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGc                      <  1:1124288/69‑1 (MQ=255)
                                        aGCAGTGGCTTTGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGc                      <  1:1947570/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:2403595/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:1146687/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:549662/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:1740656/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:2044235/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:3191951/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:3098131/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:2240040/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:28432/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:277619/69‑1 (MQ=255)
                                                 tttGGCGAGGTTGAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGATCCGGCAGCCGGGCCaaa             <  1:2756437/69‑1 (MQ=255)
                                                     gCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCaaata           <  1:2149793/67‑3 (MQ=255)
                                                             aaaTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGa           >  1:1412153/1‑59 (MQ=255)
                                                             aaaTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCcatcat  >  1:3273385/1‑68 (MQ=255)
                                                             aaaTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCcatcat  >  1:3351158/1‑68 (MQ=255)
                                                             aaaTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCcatcat  >  1:2259480/1‑68 (MQ=255)
                                                             aaaTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCcatcat  >  1:702794/1‑68 (MQ=255)
                                                             aaaTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCcatcat  >  1:1106164/1‑68 (MQ=255)
                                                             aaaTCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCcatcat  >  1:960841/1‑68 (MQ=255)
                                                             |                                                                   
AACCGCCGATGTTTGACGTTGCGTTGAGTAATTTGGCGTGAGCCGTGGCTTTGGCGAGGTTAAATCCGCACAGCGTAACGAAGGCCAGCGCGTAAAACGACCCGGCAGCCGGGCCAAAGAATCCATCAT  >  minE/1530221‑1530349

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