Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,561,939 G→A 66.3% I161I (ATC→ATT lldR ← DNA‑binding transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,561,9390GA66.3% 57.5 / 42.1 77I161I (ATC→ATTlldRDNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (9/17);  new base A (18/33);  total (27/50)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.25e-01

GACGGCTATGCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGATA  >  minE/2561875‑2561999
                                                                |                                                            
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:1615381/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:2248837/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:2040530/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:1991730/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:834881/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:1903091/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:1824898/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:3397640/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:3087645/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:455750/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:148076/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:2554013/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:1427998/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:512866/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:1353516/67‑1 (MQ=255)
ggcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:1318795/67‑1 (MQ=255)
gACGGCTATGCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGt                                                          <  1:1540256/69‑1 (MQ=255)
gACGGCTATGCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGt                                                          <  1:3245771/69‑1 (MQ=255)
gACGGCTATGCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGt                                                          <  1:1899147/69‑1 (MQ=255)
gACGGCTATGCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGt                                                          <  1:3087753/69‑1 (MQ=255)
gACGGCTATGCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGt                                                          <  1:3115338/69‑1 (MQ=255)
 gcGGCTATGCTTCACAGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAACACAATGt                                                          <  1:1511998/67‑1 (MQ=255)
 aCGGCTATGCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGt                                                          <  1:1082850/68‑1 (MQ=255)
     cTATGCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTAtgt                                                     <  1:3154053/69‑1 (MQ=255)
         gCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAgg                                                 >  1:3004532/1‑69 (MQ=255)
                  aGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCaa                                        >  1:1438114/1‑69 (MQ=255)
                  aGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCaa                                        >  1:864236/1‑69 (MQ=255)
                  aGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCaa                                        >  1:2861526/1‑69 (MQ=255)
                  aGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCaa                                        >  1:1139805/1‑69 (MQ=255)
                  aGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCaa                                        >  1:159272/1‑69 (MQ=255)
                  aGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCaa                                        >  1:460194/1‑69 (MQ=255)
                  aGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCa                                         >  1:1411732/1‑68 (MQ=255)
                              cATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGaa                            <  1:3264253/69‑1 (MQ=255)
                              cATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGaa                            <  1:1261891/69‑1 (MQ=255)
                              cATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGaa                            <  1:2239404/69‑1 (MQ=255)
                                  gaagaaACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAaacg                        <  1:579311/69‑1 (MQ=255)
                                  gaagaaACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAaacg                        <  1:2795586/69‑1 (MQ=255)
                                  gaagaaACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAaacg                        <  1:228173/69‑1 (MQ=255)
                                  gaagaaACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAaacg                        <  1:1678752/69‑1 (MQ=255)
                                  gaagaaACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAaacg                        <  1:1320019/69‑1 (MQ=255)
                                  gaagaaACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAaacg                        <  1:1238603/69‑1 (MQ=255)
                                     gaaACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAaacgaac                     <  1:1685031/69‑1 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:3175380/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:382436/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:481361/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:766880/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:2137769/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:2689719/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:2648586/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:2625245/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:1299573/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:2403425/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:1478204/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:2349771/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:811426/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:164735/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:2013543/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:1907620/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:1780525/1‑64 (MQ=255)
                                          cgcgCATGGTTTGCAGCAACACAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAaacgaac                     >  1:178405/1‑64 (MQ=255)
                                                  gTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAg        >  1:2665685/1‑69 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:695380/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:2331824/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:877092/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:797324/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:790659/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:950059/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:581141/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:3341/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:2808090/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:2597336/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:2355790/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:2210136/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:1670561/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:1484263/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:1111163/68‑3 (MQ=255)
                                                         gcaacaCAATGTTATGTGAAGCTTCGGCAATCGACAGATGAAAACGAACGTCCGCTTGTGAGGCGaga  <  1:197938/68‑3 (MQ=255)
                                                                |                                                            
GACGGCTATGCTTCACTGAGGATTGCAGGACATCGAAGAAACCGCGCATGGTTTGCAGCAGCACGATGTTATGTGAGGCTTCGGCAATCGCCAGATGAAAACGAACGTCCGCTTGTGAGGCGATA  >  minE/2561875‑2561999

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: