Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,342,948 T→C 33.3% D225D (GAT→GAC glnL → sensory kinase in two‑component regulatory system with GlnG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*minE2,342,9480TC33.3% 40.8 / 29.1 57D225D (GAT→GACglnLsensory kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base T (15/23);  new base C (7/12);  total (22/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGGTGTCGATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAATATTGTG  >  minE/2342895‑2343011
                                                     |                                                               
aCGCTGGTGTCGATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGaa                                                  >  1:2698530/1‑69 (MQ=255)
aCGCTGGTGTCGATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGaa                                                  >  1:2270199/1‑69 (MQ=255)
aCGCTGGTGTCGATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGaa                                                  >  1:1795378/1‑69 (MQ=255)
aCGCTGGTGTCGATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGaa                                                  >  1:2274089/1‑69 (MQ=255)
           gATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGa                                       <  1:291740/69‑1 (MQ=255)
           gATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGa                                       >  1:1284414/1‑69 (MQ=255)
           gATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGa                                       <  1:2782090/69‑1 (MQ=255)
           gATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGa                                       >  1:578797/1‑69 (MQ=255)
            aTGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGa                                       <  1:2379125/68‑1 (MQ=255)
               gAAGTGCCGGACAACGTGCGGTTGATTCGTGATTAAGATCCCAGCCTACCGGa                                                   >  1:587286/1‑53 (MQ=255)
               gAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGa                                                   >  1:2046016/1‑53 (MQ=255)
               gAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGa                                                   >  1:2924895/1‑53 (MQ=255)
               gAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGa                                                   >  1:725157/1‑53 (MQ=255)
                 aCTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGa                                                   <  1:400153/51‑1 (MQ=255)
                    gCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCa                              <  1:1679110/69‑1 (MQ=255)
                    gCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCa                              <  1:638745/69‑1 (MQ=255)
                    gCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCa                              <  1:150300/69‑1 (MQ=255)
                    gCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCa                              <  1:2862467/69‑1 (MQ=255)
                         aCAACGTGCGGTTGATTCGTGATTACGATCCCCGCCTACCGGAACTGGCGCACGACCCCGATCAAATTg                         <  1:1173817/69‑1 (MQ=255)
                         aCAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTg                         <  1:2877280/69‑1 (MQ=255)
                         aCAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTg                         <  1:3303258/69‑1 (MQ=255)
                         aCAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTg                         <  1:3234438/69‑1 (MQ=255)
                         aCAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTg                         <  1:540998/69‑1 (MQ=255)
                           aaCGTGCGGTTATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAAc                      <  1:1518654/69‑1 (MQ=255)
                            aCGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAAc                      <  1:1865905/69‑1 (MQ=255)
                            aCGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAAc                      <  1:588301/69‑1 (MQ=255)
                                      gATTCGGGATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATc                               >  1:203314/1‑50 (MQ=37)
                                      gATTCGGGATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGc             >  1:951633/1‑68 (MQ=39)
                                      gATTCGGGATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGc             >  1:711355/1‑68 (MQ=39)
                                      gATTCGGGATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGc             >  1:2680112/1‑68 (MQ=39)
                                      gATTCGGGATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGc             >  1:1789679/1‑68 (MQ=39)
                                      gATTCGGGATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGc             >  1:2601149/1‑68 (MQ=39)
                                      gATTCGGGATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGc             >  1:1927026/1‑68 (MQ=39)
                                          cGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGaa         >  1:1374954/1‑68 (MQ=255)
                                          cGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCAGGATCAAATTGAACAGGTCTTGcag           >  1:1163404/1‑64 (MQ=255)
                                            tGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAata       <  1:2253093/68‑1 (MQ=255)
                                            tGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAata       <  1:2092741/68‑1 (MQ=255)
                                            tGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAata       <  1:2788822/68‑1 (MQ=255)
                                            tGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAata       <  1:3006337/68‑1 (MQ=255)
                                            tGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAata       <  1:1613702/68‑1 (MQ=255)
                                            tGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAata       <  1:1662084/68‑1 (MQ=255)
                                            tGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAata       <  1:50291/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:1577098/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:167477/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:1421446/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:1800124/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:1779150/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:42925/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:1904491/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:3136979/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:219469/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:2773135/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:2488490/68‑1 (MQ=255)
                                             gATTACGACCCCAGCCTGCCGGAACTGGAGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAatat      <  1:827836/68‑1 (MQ=255)
                                                 aCGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAATATtgtg  >  1:541576/1‑68 (MQ=255)
                                                 aCGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAATATtgtg  >  1:3286091/1‑68 (MQ=255)
                                                 aCGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAATATtgtg  >  1:2388473/1‑68 (MQ=255)
                                                     |                                                               
ACGCTGGTGTCGATGGAACTGCCGGACAACGTGCGGTTGATTCGTGATTACGATCCCAGCCTACCGGAACTGGCGCACGACCCGGATCAAATTGAACAGGTCTTGCTGAATATTGTG  >  minE/2342895‑2343011

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: