Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,553,151 G→A 100% A33A (GCC→GCT)  dsdC ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,553,1510GA84.5% 48.7 / 1.8 32A33A (GCC→GCT) dsdCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base A (11/16);  total (13/19)
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.01e-01
Rejected as polymorphism: E-value score below prediction cutoff.

TCCAGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCCAGGGCGAAGGACTGATGCCTGGCAGCCACTTCAAAAGTATGCATTTTTGATAATTG  >  minE/1553089‑1553204
                                                              |                                                     
tcCAGCTGATTGATACGGTGACTTACCGCACT‑GGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:2099443/68‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:909825/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:3338504/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:3295857/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:3267985/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:3147615/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:525318/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:2952699/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:2678128/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:71628/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:2429192/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:2024564/66‑1 (MQ=255)
   aGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 <  1:2008296/66‑1 (MQ=255)
    gCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCCAGGGCGAAGGACt                                             <  1:794936/69‑1 (MQ=255)
       gattgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCCAGGGCGAAGGACTGAt                                          >  1:2192992/1‑69 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGGCAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:451935/1‑69 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAg                                                 >  1:3384092/1‑60 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:692230/1‑69 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:422114/1‑69 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:729552/1‑69 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:1686778/1‑69 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:3262771/1‑69 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:3037143/1‑69 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:1717023/1‑69 (MQ=255)
         ttgatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:1348835/1‑69 (MQ=255)
         tagatACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCGAGAGCGAAGGACTGATGc                                        >  1:2003452/3‑69 (MQ=255)
                 gTGACTTACCGCACTTGGGCTCAGCGACAATTCCTCTGCCGCCAGGGCGAAGTACTGATGCCTGGCAGc                                >  1:2818128/1‑69 (MQ=255)
                 gTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCCAGGGCGAAGGACTGATGCCTGGCAGc                                <  1:939444/69‑1 (MQ=255)
                                     tCAGCGACAATTCCTCTGCCGCCAGGGCGAAGGACTGATGCCTGGCAGCCACTTCAAAAGTATGCAttt            <  1:2964426/69‑1 (MQ=255)
                                               ttCCTCTGCCGCGAGAGCGAAGGACTGATGCCTGGCAGCCACTTCAAAAGTATGCAGTTTTGATAGTTg  <  1:343325/69‑1 (MQ=255)
                                               ttCCTCTGCCGCGAGAGCGAAGGACTGATGCCTGGCAGCCACTTCAAAAGTATGCAGTTTTGATAGTTg  <  1:590783/69‑1 (MQ=255)
                                               ttCCTCTGCCGCGAGAGCGAAGGACTGATGCCTGGCAGCCACTTCAAAAGTATGCAGTTTTGATAGTTg  <  1:2519232/69‑1 (MQ=255)
                                                              |                                                     
TCCAGCTGATTGATACGGTGACTTACCGCACTGGGGCTCAGCGACAATTCCTCTGCCGCCAGGGCGAAGGACTGATGCCTGGCAGCCACTTCAAAAGTATGCATTTTTGATAATTG  >  minE/1553089‑1553204

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: ‑