Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,650,803 A→G 100% R235R (CGT→CGC hda ← ATPase regulatory factor involved in DnaA inactivation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,650,8030AG78.4% 109.1 / 20.5 78R235R (CGT→CGChdaATPase regulatory factor involved in DnaA inactivation
Reads supporting (aligned to +/- strand):  ref base A (16/1);  new base G (61/0);  total (77/1)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.61e-01
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATT  >  minE/1650735‑1650863
                                                                    |                                                            
gccgcATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTa                                                              <  1:648661/69‑1 (MQ=255)
gccgcATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTa                                                              <  1:2879680/69‑1 (MQ=255)
gccgcATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTa                                                              <  1:2700897/69‑1 (MQ=255)
gccgcATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTa                                                              <  1:2090460/69‑1 (MQ=255)
 ccgcATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCCTTCACAAACGGAATGGTCAGCTTa                                                              <  1:1651548/68‑1 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTgc                                                        >  1:3176559/1‑59 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGt                                                   >  1:1382905/1‑64 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATc                                               >  1:1823325/1‑68 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATc                                               >  1:1173262/1‑68 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATc                                               >  1:2625321/1‑68 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATc                                               >  1:1920685/1‑68 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:2327444/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:1783275/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:184771/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:1910099/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:1922601/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:2015073/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:2072844/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:252150/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:722005/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:728666/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:750395/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:939961/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:961580/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:1025188/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:1039699/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:1610786/1‑69 (MQ=255)
                aaCACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCGTTGCGCGGTAATCg                                              >  1:1229254/1‑69 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2640712/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2733622/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2730052/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1770714/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2740155/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2744221/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:276654/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2781321/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2819659/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2909963/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:3155812/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1289635/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:327460/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:412978/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:469020/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:594597/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:693416/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:709794/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1167604/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:879922/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2244055/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:170759/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1556143/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1951150/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1506491/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2076383/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1406858/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1365871/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2246543/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1379202/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2557611/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2551090/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:1003092/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2541422/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTGCg                                                            >  1:2525728/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCa                                         >  1:649687/1‑68 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCAc                                        >  1:385194/1‑69 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCAc                                        >  1:687369/1‑69 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCAc                                        >  1:2162005/1‑69 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCAc                                        >  1:1536533/1‑69 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCACGGTAATCGACGCAc                                        >  1:1332726/1‑69 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCACCTTGCg                                                            >  1:2382570/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTAAGCTTGCg                                                            >  1:1198970/1‑49 (MQ=255)
                      ttATCTACAACTTCAGAATTTCTTTAACAAACGGAATGGTCAGCTTGCg                                                            >  1:2521500/1‑49 (MQ=255)
                                   cAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGAt                            >  1:1493101/1‑68 (MQ=255)
                                   cAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGAt                            >  1:1618636/1‑68 (MQ=255)
                                         ttctttcACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGt                     <  1:1223176/69‑1 (MQ=255)
                                          tctttcACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCa                                  >  1:2680515/1‑55 (MQ=255)
                                          tctttcACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCa                                  >  1:2707017/1‑55 (MQ=255)
                                          tctttcACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCa                                  >  1:660184/1‑55 (MQ=255)
                                          tctttcACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCa                                  >  1:3403964/1‑55 (MQ=255)
                                          tctttcACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCa                                  >  1:238862/1‑55 (MQ=255)
                                                            gTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCAtt  >  1:3053901/1‑69 (MQ=255)
                                                            gTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCAtt  >  1:1291075/1‑69 (MQ=255)
                                                            gTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCAtt  >  1:2786487/1‑69 (MQ=255)
                                                                    |                                                            
GCCGCATCCGACAATAAACACCTTATCTACAACTTCAGAATTTCTTTCACAAACGGAATGGTCAGCTTACGTTGCGCGGTAATCGACGCACGATCCAACTGATCCAACGTCATAAATAGCGTGCGCATT  >  minE/1650735‑1650863

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