Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,562,170 C→G 15.4% intergenic (‑80/‑89) pctA ← / → PP_1372 methyl‑accepting chemotaxis protein PctA/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,562,1700CG15.4% 56.3 / 5.1 26intergenic (‑80/‑89)pctA/PP_1372methyl‑accepting chemotaxis protein PctA/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (10/11);  new base G (2/2);  total (12/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.84e-01

GAGTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGGC  >  NC_002947/1562028‑1562316
                                                                                                                                              |                                                                                                                                                   
gAGTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCtttt                                                                                                                                               >  2:245670/1‑149 (MQ=255)
   ccgatctaTCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCgg                                                                                                                                           <  7:162901/142‑1 (MQ=255)
     tCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGGGACTTTTCggga                                                                                                                                         <  5:105950/149‑3 (MQ=255)
           tCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCgg                                                                                                                                           >  8:162901/1‑141 (MQ=255)
                         aaaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGc                                                                                                                       <  4:300019/147‑1 (MQ=255)
                         aaaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGc                                                                                                                       >  3:300019/1‑147 (MQ=255)
                          aaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCctctc                                                                                                                                >  1:400/1‑137 (MQ=255)
                          aaaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCctctc                                                                                                                                <  2:400/137‑1 (MQ=255)
                           aaGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAAt                                                                                                                   <  8:64265/149‑1 (MQ=255)
                             gAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATaa                                                                                                                 <  1:245670/149‑1 (MQ=255)
                                        cTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAgg                                                                                                       >  6:200060/1‑148 (MQ=255)
                                         tGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGGGACTTTTCgg                                                                                                                                           >  6:63316/1‑111 (MQ=255)
                                                       ttGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTaa                                                                                        <  8:273917/148‑1 (MQ=255)
                                                          ttCATCGCGTACTCGCTTCCGCTAAGCGTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCt                                                                                     <  6:82019/148‑1 (MQ=255)
                                                                     aTCGGTTCCGCTATTGTTTTATATGAGGCTTTCGCACGAGTTGAGCCGTAGGTTGTGCAACATCCAGAAAATTACCTTTTCGGCTGTTCCTATCTATTAAGGCTAATAACACCCACAGGTATCGGCCTAATTAATCTAAAGCTTGAga                                                                          <  4:97002/148‑2 (MQ=255)
                                                                             cGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAg                                                                                                        <  1:162146/110‑1 (MQ=255)
                                                                             cGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAg                                                                                                        >  2:162146/1‑110 (MQ=255)
                                                                                              cggCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAg                                                                                 <  4:278313/115‑1 (MQ=255)
                                                                                                                      aGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCAGGc                        >  1:280648/1‑149 (MQ=255)
                                                                                                                               aaCATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCACTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGAGCGGGCCCGGaagtc               >  2:234344/1‑145 (MQ=255)
                                                                                                                                  aTCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgccag            >  7:304950/1‑149 (MQ=255)
                                                                                                                                   tCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGgccagc           <  8:304950/149‑1 (MQ=255)
                                                                                                                                    gccGAAAAGGGACTTTTCGGGATGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCATTGTAGGCGCGGGCCCGGAt                  <  1:342577/141‑1 (MQ=255)
                                                                                                                                    gccGAAAAGGGACTTTTCGGGATGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCATTGTAGGCGCGGGCCCGGAt                  >  2:342577/2‑142 (MQ=255)
                                                                                                                                            gTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCTAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAATGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGgc  <  1:234344/149‑1 (MQ=255)
                                                                                                                                             tCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGgc  >  7:316188/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                   
GAGTATCACGATCAGTGATGCCGCCAAAAGAATCTTGTGGCTGAAACGAAGGCTTTTGTTCATGGCGTAATCGGTTCCGCTAAGGTTTTATATCAGGCGTTCGCACGAGTGGAGCCGTAGGTTGTGCAACATCCCGAAAAGTCCCTTTTCGGCTGTTCCTCTCTATTAAGGCGAATAACACCCACAGGTATCGGCCGAATTAAGCTAAAGCTTGAGCAGAGGATTAAAGAGATGTCGGAATTTTCGACCATCGTTGTAGGCGCGGGCCCGGATGGCCAGCCAGTCGGGC  >  NC_002947/1562028‑1562316

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