Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,968,245 C→T 11.8% intergenic (‑112/+23) glnK ← / ← PP_5723 NRII(GlnL/NtrB) phosphatase activator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,968,2450CT11.8% 83.1 / 3.6 34intergenic (‑112/+23)glnK/PP_5723NRII(GlnL/NtrB) phosphatase activator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/12);  new base T (2/2);  total (20/14)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.54e-01

AGCTTGAACGGCTTGATGATGGCTGTGACTAGCTTCATGAAACTCTCTCCCGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTC  >  NC_002947/5968096‑5968393
                                                                                                                                                     |                                                                                                                                                    
aGCTTGAACGGCTTGATGATGGCTGTGACTAGCTTCATGAAACTCTCTCCCGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCAc                                                                                                                                                      >  2:191323/1‑150 (MQ=255)
     gAACGGCTTGATGATGGCTGTGACTAGCTTCATGAAACTCTCTCCCGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCt                                                                                                                                                 >  1:309805/1‑150 (MQ=255)
       aCGGCTTGATGATGGCTGTGACTAGCTTCATGAAACTCTCTCCCGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTgg                                                                                                                                               <  2:613295/150‑1 (MQ=255)
                       tgtgACTAGCTTCATGAAACTCTCTCCCGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCtt                                                                                                                               <  1:359480/150‑1 (MQ=255)
                                     tGAAACTCTCTCCCGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCg                                                                                                                 >  2:462433/1‑150 (MQ=255)
                                              ctcCCGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCt                                                                                                        <  2:309805/150‑1 (MQ=255)
                                               tcCCGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCtt                                                                                                       >  2:145645/1‑150 (MQ=255)
                                                 ccGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTg                                                                                                     >  1:153628/1‑150 (MQ=255)
                                                   gATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCa                                                                                                   <  1:191323/150‑1 (MQ=255)
                                                         gtggACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCAt                                                                                                                       <  1:330483/124‑1 (MQ=255)
                                                         gtggACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCAt                                                                                                                       >  2:330483/1‑124 (MQ=255)
                                                                  ccccAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGtt                                                                                    >  1:127444/1‑150 (MQ=255)
                                                                    ccAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTcc                                                                                  >  1:433245/1‑150 (MQ=255)
                                                                      aGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCATTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGa                                                                                >  2:590033/1‑150 (MQ=255)
                                                                      aGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCATTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGa                                                                                >  2:589715/1‑150 (MQ=255)
                                                                              aaacCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAgcgc                                                                        >  1:566856/1‑150 (MQ=255)
                                                                              aaacCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAgcgc                                                                        <  2:566856/150‑1 (MQ=255)
                                                                                acCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCAt                                                                                                                       >  2:396135/1‑101 (MQ=255)
                                                                                acCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCAt                                                                                                                       <  1:396135/101‑1 (MQ=255)
                                                                                  ccGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAg                                                                                                                                        >  2:154152/1‑82 (MQ=255)
                                                                                  ccGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAg                                                                                                                                        <  1:154152/82‑1 (MQ=255)
                                                                                           gTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAgcgc                                                           <  1:145645/150‑1 (MQ=255)
                                                                                            tCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGcgcg                                                          >  1:259623/1‑150 (MQ=255)
                                                                                            gcTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGcgcg                                                          >  2:155307/2‑150 (MQ=255)
                                                                                                         gTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGcc                                                                       <  1:589934/124‑1 (MQ=255)
                                                                                                         gTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGcc                                                                       >  2:589934/1‑124 (MQ=255)
                                                                                                            ggCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCg                                          <  2:259623/150‑1 (MQ=255)
                                                                                                             gCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCATTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCgg                                         <  1:590033/150‑1 (MQ=255)
                                                                                                             gCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCATTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCgg                                         <  1:589715/150‑1 (MQ=255)
                                                                                                                   tAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCGGGATTgg                                   >  2:277708/1‑150 (MQ=255)
                                                                                                                     aCGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACCGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCGGGATTGGcc                                 >  2:405582/1‑150 (MQ=255)
                                                                                                                       gcgTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACCGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCGGGATTGGCCAg                               <  1:405582/150‑1 (MQ=255)
                                                                                                                                gCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCa                      >  2:565020/1‑150 (MQ=255)
                                                                                                                                            gcgTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGc          >  2:254531/1‑150 (MQ=255)
                                                                                                                                                    aCTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTc  >  2:27671/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
AGCTTGAACGGCTTGATGATGGCTGTGACTAGCTTCATGAAACTCTCTCCCGATTTGGTGGACTTGCCCCAGGAAAACAAACCCGTCTCAAGTCTAAGCGCAGCGGTTGGCTTTGTAACGCGTCGTCGGCATCCGGCTCCGCGTGCGCACTGCCTGGTCACAAGGAACTGCTTCAGTGCATGGCTCGTACAGGTCTTTGCAGAAACCTTGCCAGTTCCGACAAAGCGCCGAAAAACAGCGCGTTATGTGATTTGCTCGGGATTGGCCAGTCGCCAGCACGGGAAACATGCACAATTTC  >  NC_002947/5968096‑5968393

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