Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,100,007 A→C 21.6% intergenic (‑118/‑87) PP_1875 ← / → PP_1876 two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,100,0070AC21.6% 59.7 / 3.0 31intergenic (‑118/‑87)PP_1875/PP_1876two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (14/9);  new base C (5/2);  total (19/12)
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

CGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAAC  >  NC_002947/2099859‑2100154
                                                                                                                                                       |                                                                                                                                                      
cGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCa                                                                                                                                                        <  5:355786/149‑1 (MQ=255)
    ttCCCGTTCCTTGCGCTGTGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGACCGGTCACCCAGGTCACGCATAGCCCCAggag                                                                                                                                                    >  7:160967/1‑147 (MQ=255)
        cGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAAc                                                                                                                                                <  8:224712/149‑1 (MQ=255)
        cGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAAc                                                                                                                                                <  2:40988/149‑1 (MQ=255)
           gccTTCCGCTGCGCAAGCATAGCCGAATCGCCCCGGCTGCGGGTGGGGCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGc                                                                                                                                              <  6:206655/147‑1 (MQ=255)
             cTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCaa                                                                                                                                            <  8:190444/148‑1 (MQ=255)
                       gcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGa                                                                                                                                <  2:74917/149‑1 (MQ=255)
                                 ccGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCgg                                                                                                                            <  6:323892/143‑1 (MQ=255)
                                 ccGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCgg                                                                                                                            >  5:323892/1‑143 (MQ=255)
                                 ccGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCAGCGGGCGTGCCCGTTCCACCAGGGTCAGCATAGCCCAAGGCTGAACGCAAGCCCTGAGGGACCGGGCacca                                                                                                                      >  4:244181/1‑146 (MQ=255)
                                                        ggTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTTCCCGGTCCCCCAGGGACAGCTTAGCCCCAGGCGGAACGCCAGCCCTGTGGGCGCGGGCAAGCCCGCGGAGGGGCGCGCCGCGGcc                                                                                             >  6:367624/1‑149 (MQ=38)
                                                             gtctgtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCATTGCCGTTACCCCACCTAAACCTTAACCCCCCGCGGAACCCACGCCCGGTGGGAGCCGGGCAGCCCCCGGCGGGGCGCGCAGCCGGccccc                                                                                          >  4:117100/1‑148 (MQ=255)
                                                              tctgtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTTCCCGTTACCCCAGGTACAGCGTTGCCCCCCGGGGCACGCCCGCCCGGTGGGGGCCGGCCAGCCCCCGGCGGGGCGCCCCGGGGcccccca                                                                                        >  5:381078/1‑146 (MQ=255)
                                                                  ttCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGGGTGGCCCGTAACCCACGGGAAGCATTGGCCCCCGGGGGAACCAAGCCCGGTGGGGGCGGGGCACCCCCCGGCGGGGCGCGGCGCGGCCCCCcaaaaa                                                                                    >  8:6738/1‑148 (MQ=255)
                                                                  ttCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTACCAAAGGTACAGCTTAGCCCCCGGCGGGACGCAAGCCCGGTGGGCGCGGGCCAGCCCGCGGAGGGGCGCGCCGCGGCCCCCAAAAAAc                                                                                   >  7:71097/1‑149 (MQ=255)
                                                                                                   gggCGCCCCCCCGCGGGGTGGCCCGTCCCCCCCGGTCCCGCTTCCCCCCCGCCGTACCCCACCCCTGGGGGGCGGGGGAACGCCCGCGGAGGGGCGCGCAGTGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGAcgctcgc                                                   <  6:381078/147‑1 (MQ=255)
                                                                                                         cccTCCGCGGGCTTGGCCGTTCCCACAGGTACAGCTTAGCCCCCGGCTGGACGCAAGCCCGGTGGGAGCGGGGCAGCCCGCGGGGGGGCGCGCCGGGGCCCCCCAAAAACTGAGGTTGTTTGGGGTCGTGATGGACGCTCGCGTGTTCg                                            >  2:250996/1‑149 (MQ=255)
                                                                                                                     cTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGAc                                 >  5:266577/1‑149 (MQ=255)
                                                                                                                      ttGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCACGCCCGCGGCGGGGCGCGCAGCGGCCCCCAAAAACCTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGAcc                                >  1:6805/1‑149 (MQ=255)
                                                                                                                           cGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGcc                           >  4:297387/1‑149 (MQ=255)
                                                                                                                                  cacaGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGccagtacc                     >  5:226916/1‑143 (MQ=255)
                                                                                                                                  cacaGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCg                     >  8:228700/1‑149 (MQ=255)
                                                                                                                                   acaGGTACAGCTTAGCCCCAGGCTGAACGCAAACCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCgg                    >  1:233472/1‑149 (MQ=255)
                                                                                                                                        tACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGAAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTg               >  5:343162/1‑149 (MQ=255)
                                                                                                                                           aGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGgcg            >  2:232672/1‑149 (MQ=255)
                                                                                                                                               tttggcccTGGCTCAAGGCAAGCCCTGGGGGTGGGGCAACGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGgcgcgc         <  5:345584/144‑1 (MQ=255)
                                                                                                                                                   ccccCGGCTGAACGCACGCCGTGGGGGAGCGGGAAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTg     <  5:367624/149‑1 (MQ=255)
                                                                                                                                                   ccccAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCtt      >  3:4992/1‑148 (MQ=255)
                                                                                                                                                     ccAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTg     >  7:273931/1‑147 (MQ=255)
                                                                                                                                                     ccAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTg     <  8:273931/147‑1 (MQ=255)
                                                                                                                                                     ccAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGaa   <  1:232672/149‑1 (MQ=255)
                                                                                                                                                      cAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAAc  <  6:266577/149‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                      
CGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTCGTCATGGACGCTCGCTTGATCGCTTTTCTGGACCGCGCCGAATCGGTACTGGCGCGCCTTGAAC  >  NC_002947/2099859‑2100154

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