Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 571,621 T→C 23.5% intergenic (+44/+47) ssb → / ← PP_0486 single‑stranded DNA‑binding protein/GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947571,6210TC23.5% 32.0 / 5.6 21intergenic (+44/+47)ssb/PP_0486single‑stranded DNA‑binding protein/GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (7/9);  new base C (2/3);  total (9/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 = 9.08e-01

CAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGC  >  NC_002947/571479‑571760
                                                                                                                                              |                                                                                                                                           
cAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTg                                                                                                                                      <  2:199664/150‑1 (MQ=255)
 aGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGATCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTgg                                                                                                                                     <  1:449225/150‑1 (MQ=255)
   cgcCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTGGGt                                                                                                                                   <  1:426244/150‑1 (MQ=255)
       cGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGttttt                                                                                                                               <  2:496995/150‑1 (MQ=255)
                       gcgcccagcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACcgc                                                                                                               <  2:425000/150‑1 (MQ=255)
                         gcccagcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGt                                                                                                             >  2:73569/1‑150 (MQ=255)
                         gcccagcgccGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATAGGCAGCACCGCGt                                                                                                             >  1:412199/1‑150 (MQ=255)
                                     aaCAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGc                                                                                                 <  1:73569/150‑1 (MQ=255)
                                         gCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTc                                                                                             <  2:311571/150‑1 (MQ=255)
                                                  gcAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGaaca                                                                                    >  1:488189/1‑150 (MQ=255)
                                                       cggcggCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCgg                                                                               <  1:230610/150‑1 (MQ=255)
                                                       cggcggCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATAGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCgg                                                                               <  2:412199/150‑1 (MQ=255)
                                                        ggcggcTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTGGGtttct                                                                                                                               <  1:162348/101‑3 (MQ=255)
                                                        ggcggcTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAACTGCCCTGGGtttct                                                                                                                               >  2:162348/1‑99 (MQ=255)
                                                             gggACTTCGACAGCTTCGATGACGAGATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACACAGGGCAGTTGCTTTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTTTGTAGAACAACCGGGGCtgc                                                                         <  2:196779/148‑3 (MQ=255)
                                                                                                    cAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAg                                  <  1:184248/150‑1 (MQ=255)
                                                                                                               ggCGGTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTg                       >  1:411886/1‑150 (MQ=255)
                                                                                                                      aGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAGCCTGGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTg                >  1:405208/1‑150 (MQ=255)
                                                                                                                              cAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAAACGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCAt        >  1:520659/1‑150 (MQ=255)
                                                                                                                               aGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCGAGAGGGATTGTCAGCTTGGTCAGGGCAAAGCCCAGGTCTTGGCGGCTGGGTGCCATc       >  2:14697/1‑150 (MQ=255)
                                                                                                                                    cccAGGGCAACTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGc  >  2:547787/1‑150 (MQ=255)
                                                                                                                                              |                                                                                                                                           
CAGCGCCCGCAACAGGCCCCTCAGCGCCCAGCGCCGCAACAGCCTGCGCCGCAGCCGGCGGCTGACTTCGACAGCTTCGATGACGATATTCCGTTCTGATCAGAGCGGTAGGGCGTTCAGCAAGACCAGAAACCCAGGGCAGTTGCCCTGGGTTTTTTTATGGGCAGCACCGCGTATCAGGCATCGCTTTCAGTAGAACAACCGGGGCTCCGGCGAATCCAGCAGGCTTGCCAGCTTGGTCAGGGCAAAGCCCAGTTCTTGGCGGCTGGGTGCCATCAGGGC  >  NC_002947/571479‑571760

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