Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,221,231 G→T 21.0% intergenic (+46/+216) PP_3699 → / ← PP_3700 hypothetical protein/chromosome partitioning ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,221,2310GT21.0% 53.8 / 8.3 29intergenic (+46/+216)PP_3699/PP_3700hypothetical protein/chromosome partitioning ATPase
Reads supporting (aligned to +/- strand):  ref base G (12/10);  new base T (3/3);  total (16/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

ACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACAT  >  NC_002947/4221098‑4221371
                                                                                                                                     |                                                                                                                                            
aCGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccctga                                                                                                                              <  2:618813/150‑5 (MQ=255)
  gTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctc                                                                                                                            >  1:497148/1‑150 (MQ=255)
    cGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAAcctcag                                                                                                                          <  1:561088/150‑1 (MQ=255)
          gACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCg                                                                                                                                              >  1:1083941/1‑124 (MQ=255)
          gACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCg                                                                                                                                              <  2:1083941/124‑1 (MQ=255)
               tGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAATGAATTGTCAAATGCTGACCGTCGTGCATATCGCGCTCAGTGGGCATCTTACAGATGCCCCAAccccatcctcagcctca                                                                                                               >  1:248611/1‑150 (MQ=255)
                      gCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGTGGGCATCg                                                                                                                                              <  1:365719/112‑1 (MQ=255)
                      gCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGTGGGCATCg                                                                                                                                              >  2:365719/1‑112 (MQ=255)
                      gCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCaa                                                                                                        <  1:1044250/150‑1 (MQ=255)
                            gATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTa                                                                                                  <  1:272962/150‑1 (MQ=255)
                                cTGAAGCTTTGATCAGGGCTTATATGCTTGGGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGGTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTa                                                                                              >  2:486155/1‑150 (MQ=255)
                                             cAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc                                                                                                             >  1:1031678/1‑122 (MQ=255)
                                             cAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc                                                                                                             <  2:1031678/122‑1 (MQ=255)
                                             cAGGGCTGATATGCTTGAGGAGTTTCAACAAGACTGCAAGTGAAGTGGCAAATGCTAACCGTCGTACATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCCGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAgcgcgc                                                                                 >  2:867294/1‑150 (MQ=255)
                                                ggCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccct                                                                                                                                <  1:714258/100‑3 (MQ=255)
                                                ggCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCTGTCCGATGccccct                                                                                                                                >  2:714258/1‑98 (MQ=255)
                                                ggCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGa                                                                              >  2:631057/1‑150 (MQ=255)
                                                           ttGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc                                                                                                             >  1:726474/1‑108 (MQ=255)
                                                           ttGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc                                                                                                             <  2:726474/108‑1 (MQ=255)
                                                           ttGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCCAGCAGAGGCCCCAACCTCAGCCTCAGCGTCAGCATCAAAACCTAGCGATGGCCGAGCGCGCCGAAGCCACAAtct                                                                   >  1:870391/1‑150 (MQ=255)
                                                             gaggagTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCa                                                                 <  1:631057/150‑1 (MQ=255)
                                                                  gTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCaa                                                            >  2:729855/1‑150 (MQ=255)
                                                                          aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc                                                                                                             >  1:730244/1‑93 (MQ=255)
                                                                          aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACctcagcctcagcctcagc                                                                                                             <  2:730244/93‑1 (MQ=255)
                                                                          aaGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGcc                                                    >  1:674245/1‑150 (MQ=255)
                                                                                    gtgaagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGc                                                                                                <  1:794733/96‑1 (MQ=255)
                                                                                    gtgaagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGc                                                                                                >  2:794733/1‑96 (MQ=255)
                                                                                     tgaagtgGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCgg                                         <  1:729855/150‑1 (MQ=255)
                                                                                            gCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTa                                  >  1:404463/1‑150 (MQ=255)
                                                                                                                      tgaggttggggCATCTGTCCGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATg           >  1:586478/8‑147 (MQ=255)
                                                                                                                      tgaggttggggCATCTGTCCGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATg           <  2:586478/140‑1 (MQ=255)
                                                                                                                      cgcTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGccc        <  2:497148/150‑1 (MQ=255)
                                                                                                                            gggggCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACAt  >  2:51147/1‑150 (MQ=255)
                                                                                                                                     |                                                                                                                                            
ACGTCGAGTTGACGCTGTTCACGCTGAAGATCCTGAAGCTTTGATCAGGGCTGATATGCTTGAGGAGTTTGAACAAGACTGCAAGTGAAGTGGCAAATGCTGACCGTCGTGCATAGCGCGCTCAGGGGGCATCGGACAGATGCCCCAACCTCAGCCTCAGCCTCAGCCTCAAAACCTAGCTATGGCCGAGCGCGCCGAAGCCACAATCTCATGCAACCGCTGCCGTACGAGTCGGCACAGTATTGCATCTCTAGCTCCGAGCATGCCCTGACAT  >  NC_002947/4221098‑4221371

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