Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 423,948 C→T 44.4% E316E (GAG→GAA nagC ← DNA‑binding transcriptional dual regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE423,9480CT44.4% 10.3 / 40.7 54E316E (GAG→GAAnagCDNA‑binding transcriptional dual regulator
Reads supporting (aligned to +/- strand):  ref base C (23/7);  new base T (13/11);  total (36/18)
Fisher's exact test for biased strand distribution p-value = 9.29e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

GCGGATTAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTC  >  minE/423890‑424016
                                                          |                                                                    
gCGGATTAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACt                                                             <  1:2546939/68‑1 (MQ=255)
gCGGATTAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACt                                                             <  1:1589776/68‑1 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:1224392/60‑2 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:127276/60‑2 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:2140068/60‑2 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:359384/60‑2 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:2841774/60‑2 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:2722717/60‑2 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:3150095/60‑2 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:2992403/60‑2 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:1873873/60‑2 (MQ=255)
     ttAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:3090353/60‑2 (MQ=255)
      tAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATATTCAAtc                                                                <  1:3092996/59‑2 (MQ=255)
                    gCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGcc                                         >  1:2257439/1‑68 (MQ=255)
                    gCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCt                                        >  1:2375872/1‑69 (MQ=255)
                    gCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCt                                        >  1:2831701/1‑69 (MQ=255)
                    gCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCt                                        >  1:3310372/1‑69 (MQ=255)
                    gCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCt                                        >  1:2639576/1‑69 (MQ=255)
                       aTGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTg                                     >  1:2040861/1‑69 (MQ=255)
                       aTGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTg                                     >  1:611437/1‑69 (MQ=255)
                            gagggTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcg                                 >  1:1776778/4‑68 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:89124/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:2747669/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:2632527/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:2748706/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:279974/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:2831648/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:3256724/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:3260802/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:444162/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:713633/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcgc                                >  1:2583271/1‑69 (MQ=255)
                            gATGGTTTTACCCAGATGACGACCGACATATTCAATCACTTCCGAGGCCAGGCTATCGCCTTTGTTcg                                 >  1:120051/1‑68 (MQ=255)
                              tGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCgg                              >  1:3352662/1‑69 (MQ=255)
                              tGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCgg                              >  1:250918/1‑69 (MQ=255)
                              tGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCgg                              >  1:2264083/1‑69 (MQ=255)
                                   ttACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCgg                              >  1:768015/1‑64 (MQ=255)
                                            tGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtt                >  1:2987232/1‑69 (MQ=255)
                                             gacgacCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtt                >  1:14394/1‑68 (MQ=255)
                                             gacgacCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtt                >  1:1354006/1‑68 (MQ=255)
                                             gacgacCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtt                >  1:1248505/1‑68 (MQ=255)
                                              acgacCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt              <  1:2835399/69‑1 (MQ=255)
                                              acgacCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt              <  1:2321829/69‑1 (MQ=255)
                                              acgacCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt              <  1:2324239/69‑1 (MQ=255)
                                              acgacCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGtttt              <  1:2513664/69‑1 (MQ=255)
                                                 accgacATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGAt           >  1:852223/1‑69 (MQ=255)
                                                 accgacATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGAt           >  1:804161/1‑69 (MQ=255)
                                                 accgacATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGAt           >  1:2615373/1‑69 (MQ=255)
                                                 accgacATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGAt           >  1:2764638/1‑69 (MQ=255)
                                                 accgacATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGAt           >  1:470054/1‑69 (MQ=255)
                                                 accgacATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGAt           >  1:3432243/1‑69 (MQ=255)
                                                 accgacATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGAt           >  1:2662999/1‑69 (MQ=255)
                                                      cATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCgg                              <  1:228784/45‑1 (MQ=255)
                                                       aTACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTAGATGGTGCa     >  1:1135707/1‑69 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGtt                                   >  1:33088/1‑36 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGc                             >  1:1449464/1‑42 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGc                             >  1:345287/1‑42 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGc                             <  1:1496774/42‑1 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGc                             <  1:1745862/42‑1 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGc                             >  1:628983/1‑42 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGc                             >  1:1929542/1‑42 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGc                             >  1:2845229/1‑42 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGc                             >  1:2129673/1‑42 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGc                             >  1:2814588/1‑42 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAgtc  >  1:3312162/1‑69 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAgtc  >  1:166855/1‑69 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAgtc  >  1:1917738/1‑69 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAgt   >  1:1360210/1‑68 (MQ=255)
                                                          cTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAgt   >  1:576326/1‑68 (MQ=255)
                                                          |                                                                    
GCGGATTAAATAAGTTGATAGCAATGGCGATGGTTTTACCCAGATGACGACCGACATACTCAATTACTTCCGACGCCAGACTATCGCCTTTGTTCGCGGCTTTGCAGATAGTTTTGATGGTGCAGTC  >  minE/423890‑424016

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