Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,720,012 A→G 100% D523D (GAT→GAC yhfK ← conserved inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,720,0120AG83.8% 173.8 / 23.8 74D523D (GAT→GACyhfKconserved inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (0/12);  new base G (0/62);  total (0/74)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGCTCT  >  minE/2719956‑2720043
                                                        |                               
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:566748/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1816831/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1912279/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1958333/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1976395/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:2009279/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:2115201/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:2139770/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:2154356/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:2164180/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:2252245/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:267543/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:360122/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:394319/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:473349/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1040532/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:567967/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:578726/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:621371/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:651722/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:652535/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:670734/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:757653/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:766855/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:789885/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:844290/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:863968/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:926236/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:961977/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1549707/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1168975/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1235397/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1246196/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1255263/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1301331/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:136442/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1385500/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:139428/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1796293/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1138822/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1477116/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1495950/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1554530/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1011260/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1579717/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1053246/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:174765/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1624217/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1712337/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1624952/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:109169/67‑1 (MQ=255)
ggTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGGCATGGGCGt                       <  1:162189/67‑1 (MQ=255)
 ggtGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:697021/64‑1 (MQ=255)
 gTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:2018769/66‑1 (MQ=255)
 gTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:648183/66‑1 (MQ=255)
 gTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:2084172/66‑1 (MQ=255)
 gTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1453261/66‑1 (MQ=255)
 gTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:2104380/66‑1 (MQ=255)
 gTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1602980/66‑1 (MQ=255)
 gTTGCGGGTCCTCGCTCAGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:296041/66‑1 (MQ=255)
                tcaGAATCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:804504/49‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:923624/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:539380/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:1160864/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:856981/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:769927/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:736893/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:669108/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:528430/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:254192/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:156829/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:1728348/67‑1 (MQ=255)
                     aTCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGctct  <  1:1801621/67‑1 (MQ=255)
                      tCAAGCGAATCGCGTCCTGATAGGCTTCTAAAGCGTCATGGGCGt                       <  1:1789597/45‑1 (MQ=255)
                                                        |                               
GGTTGCGGATCCTCGCTAAGAATCAAGCGAATCGCTTCCTGATAGGCTTCTAAAGCATCATGGGCGTTTTTACGCAATAACCCGCTCT  >  minE/2719956‑2720043

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