Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 1,448,281 A→T 17.8% A229A (GCA→GCT atoS → sensory histidine kinase in two‑component regulatory system with AtoC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE1,448,2810AT17.8% 107.9 / 16.1 56A229A (GCA→GCTatoSsensory histidine kinase in two‑component regulatory system with AtoC
Reads supporting (aligned to +/- strand):  ref base A (17/29);  new base T (8/2);  total (25/31)
Fisher's exact test for biased strand distribution p-value = 1.72e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTCTCCCGTCGCCTGAGCGCCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAAT  >  minE/1448221‑1448347
                                                            |                                                                  
tttCTCCCGTCGCCTGAGCGCCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACTAAatat                                                            <  1:3573/69‑1 (MQ=255)
tttCTCCCGTCGCCTGAGCGCCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAatat                                                            <  1:994182/69‑1 (MQ=255)
tttCTCCCGTCGCCTGAGCGCCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAatat                                                            <  1:1421341/69‑1 (MQ=255)
tttCTCCCGTCGCCTGAGCGCCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAatat                                                            <  1:2037492/69‑1 (MQ=255)
tttCTCCCGTCGCCTGAGCGCCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAatat                                                            <  1:1306529/69‑1 (MQ=255)
tttCTCCCGTCGCCTGAGCGCCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAatat                                                            <  1:2203153/69‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:279997/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:2704929/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:339953/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:2571729/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:520123/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:2283136/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:2196333/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:2011218/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:1235836/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:1824965/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:934812/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:971081/68‑1 (MQ=255)
                  cgcCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:1714057/68‑1 (MQ=255)
                  cgcCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAATATTCCCACTCGATTACcaca                                          >  1:730101/1‑69 (MQ=255)
                  cgcCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAATATTCCCACTCGATTACcaca                                          >  1:974897/1‑69 (MQ=255)
                  cgcCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAATATTCCCACTCGATTACcaca                                          >  1:445183/1‑69 (MQ=255)
                  cgcCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAATATTCCCACTCGATTACcaca                                          >  1:696725/1‑69 (MQ=255)
                  cgcCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAATATTCCCACTCGATTACcaca                                          >  1:2684811/1‑69 (MQ=255)
                  cgcCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAATATTCCCACTCGATTACcaca                                          >  1:2237367/1‑69 (MQ=255)
                  cgcCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAATATTCCCACTCGATTACcaca                                          >  1:767831/1‑69 (MQ=255)
                  cgcCAATATAGATATCATCACCGATGGCCTCTCGACTCTGGCTCAAAATATTCCCACTCGATTACcac                                           >  1:763882/1‑68 (MQ=255)
                    ccAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:1545178/66‑1 (MQ=255)
                          ttGATATCATCACCGATGGCATCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:1993679/60‑1 (MQ=255)
                                 catcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAAt                           >  1:2412882/1‑69 (MQ=255)
                                 catcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAAt                           >  1:1021587/1‑69 (MQ=255)
                                 catcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAAt                           >  1:2780503/1‑69 (MQ=255)
                                 catcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAAt                           >  1:3108870/1‑69 (MQ=255)
                                 catcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAAt                           >  1:2601468/1‑69 (MQ=255)
                                 catcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAAt                           >  1:2333370/1‑69 (MQ=255)
                                 catcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAAt                           >  1:2282686/1‑69 (MQ=255)
                                 catcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAAt                           >  1:1933271/1‑69 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGAt                                                 >  1:1116861/1‑46 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:2293486/52‑1 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:11771/52‑1 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:1563963/52‑1 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:1668017/52‑1 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:218371/52‑1 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:3023794/52‑1 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:2562798/52‑1 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATg                          >  1:2428181/1‑69 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATg                          >  1:2608756/1‑69 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATg                          >  1:812446/1‑69 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATg                          >  1:1734525/1‑69 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATg                          >  1:965302/1‑69 (MQ=255)
                                  atcaCCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATg                          >  1:3302153/1‑69 (MQ=255)
                                  atcaCCGATGGCCTATCGACTCTGGCACAAAATATTCCCACTCGATTACcac                                           <  1:748186/52‑1 (MQ=255)
                                              ctctCGCCTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAAtca                >  1:237137/1‑67 (MQ=255)
                                              ctctCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAAtca                >  1:2013014/1‑67 (MQ=255)
                                                     ctctGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCagag         <  1:2423672/67‑1 (MQ=255)
                                                     ctctGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCagag         <  1:3224788/67‑1 (MQ=255)
                                                            acaAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTtaat  <  1:2644382/67‑1 (MQ=255)
                                                            |                                                                  
TTTCTCCCGTCGCCTGAGCGCCAATATTGATATCATCACCGATGGCCTCTCGACTCTGGCACAAAATATTCCCACTCGATTACCACAATTGCCCGGTGAAATGGGGCAAATCAGTCAGAGTGTTAAT  >  minE/1448221‑1448347

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: