Predicted mutation
evidence seq id position mutation freq annotation gene description
RA minE 2,342,351 C→T 27.5% N26N (AAC→AAT 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,3510CT27.5% 106.1 / 36.9 94N26N (AAC→AATglnLsensory kinase in two‑component regulatory system with GlnG
Reads supporting (aligned to +/- strand):  ref base C (56/12);  new base T (8/18);  total (64/30)
Fisher's exact test for biased strand distribution p-value = 3.86e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCCGATGCTGGGCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGCTCC  >  minE/2342292‑2342405
                                                           |                                                      
cccGATGCTGGGCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAATCTGGCGATc                                               >  1:1012041/1‑69 (MQ=255)
cccGATGCTGGGCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAATCTGGCGATc                                               >  1:2872828/1‑69 (MQ=255)
     tGCTGGGCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTa                                          <  1:742417/69‑1 (MQ=255)
         gggCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGcc                                      <  1:581096/69‑1 (MQ=255)
         gggCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGcc                                      <  1:3067895/69‑1 (MQ=255)
         gggCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGcc                                      <  1:2505721/69‑1 (MQ=255)
         gggCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGcc                                      <  1:2350512/69‑1 (MQ=255)
         gggCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGcc                                      <  1:2336222/69‑1 (MQ=255)
         gggCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGcc                                      <  1:2063917/69‑1 (MQ=255)
         gggCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGcc                                      <  1:52654/69‑1 (MQ=255)
         gggCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGcc                                      <  1:1345931/69‑1 (MQ=255)
         gggCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGcc                                      <  1:971082/69‑1 (MQ=255)
            cAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAAc                                   <  1:185090/69‑1 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCaa                                    >  1:3215133/1‑64 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCt                                >  1:848697/1‑68 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCt                                >  1:693041/1‑68 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCt                                >  1:1681656/1‑68 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCt                                >  1:1751792/1‑68 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCt                                >  1:2524199/1‑68 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTg                               >  1:3376683/1‑69 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTg                               >  1:1149286/1‑69 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTg                               >  1:2002602/1‑69 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTg                               >  1:3167323/1‑69 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTg                               >  1:626732/1‑69 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTg                               >  1:303227/1‑69 (MQ=255)
                tCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTg                               >  1:1306414/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCt                                >  1:1214559/1‑67 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTg                               >  1:840596/1‑68 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:2926176/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:2363795/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:334793/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:285443/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:150489/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:637912/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:1295397/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:586348/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:550472/1‑69 (MQ=255)
                 ccTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGc                              >  1:1996258/1‑69 (MQ=255)
                              aTTAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCATTACGCCAACCCCGCCGCGCAACAACTg                 >  1:922789/1‑69 (MQ=255)
                              aTTAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCATTACGCCAACCCCGCCGCGCAACAACTg                 >  1:733951/1‑69 (MQ=255)
                              aTTAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCATTACGCCAACCCCGCCGCGCAACAACTg                 >  1:1086488/1‑69 (MQ=255)
                              aTTAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCATTACGCCAACCCCGCCGCGCAACAACTg                 >  1:2351438/1‑69 (MQ=255)
                              aTTAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCATTACGCCAACCCCGCCGCGCAACAACTg                 >  1:2324716/1‑69 (MQ=255)
                              aTTAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCATTACGCCAACCCCGCCGCGCAACAACTg                 >  1:3395240/1‑69 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:103786/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:3123980/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:929825/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:416232/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:646507/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:1548591/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:2045277/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:2037087/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:2089114/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:1780985/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:1693941/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:2226212/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:1608255/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:2377954/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:2506720/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:1231801/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:1082531/40‑1 (MQ=255)
                                tAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCt               >  1:2618645/1‑69 (MQ=255)
                                tAACAGTATTTTGTGAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCt               >  1:963604/1‑69 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAATCTGGCGATCCAt                                            <  1:2658269/39‑1 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:82377/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:2690103/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:1523413/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:1483403/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:90974/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:1984646/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:57649/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:2051397/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:393317/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:2645483/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:3374273/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:477866/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:2595294/1‑61 (MQ=255)
                                 aaCAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGcaac                      >  1:2754151/1‑61 (MQ=255)
                                      tATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGcc           >  1:2223602/1‑67 (MQ=255)
                                        ttttGTTAATCGATGACAACCTGGCGATCCATTACGCCAAc                                   <  1:180590/41‑1 (MQ=255)
                                        ttttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCaaa       >  1:3406935/1‑69 (MQ=255)
                                          ttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCaaa       >  1:903116/1‑67 (MQ=255)
                                          ttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGc     >  1:2784371/1‑69 (MQ=255)
                                          ttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGc     >  1:860157/1‑69 (MQ=255)
                                          ttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGc     >  1:2891532/1‑69 (MQ=255)
                                          ttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGc     >  1:80888/1‑69 (MQ=255)
                                          ttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGc     >  1:2687852/1‑69 (MQ=255)
                                          ttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGc     >  1:2401743/1‑69 (MQ=255)
                                          ttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGc     >  1:1873100/1‑69 (MQ=255)
                                          ttGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGc     >  1:582710/1‑69 (MQ=255)
                                              tAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGCTcc  >  1:106772/1‑68 (MQ=255)
                                              tAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGCTcc  >  1:2535846/1‑68 (MQ=255)
                                              tAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGCTcc  >  1:1530242/1‑68 (MQ=255)
                                              tAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGCTcc  >  1:3261296/1‑68 (MQ=255)
                                                           |                                                      
CCCGATGCTGGGCAGATCCTCAACTCGCTGATTAACAGTATTTTGTTAATCGATGACAACCTGGCGATCCATTACGCCAACCCTGCCGCGCAACAACTGCTCGCCCAAAGCTCC  >  minE/2342292‑2342405

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: