Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 634,128:1 +C 22.2% coding (431/2016 nt) PP_0546 → sigma‑54 dependent transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947634,1281.C22.2% 41.6 / 12.4 18coding (431/2016 nt)PP_0546sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (5/9);  new base C (2/2);  total (7/11)
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.69e-01

GGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGCAAAGCCAGT  >  NC_002947/633989‑634271
                                                                                                                                            |                                                                                                                                               
ggggTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCATCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGAccatccat                                                                                                                                       <  2:368634/150‑1 (MQ=255)
     gATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGaa                                                                                                                                  <  1:31316/150‑1 (MQ=255)
                   tCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGACTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTg                                                                                                                    >  2:167101/1‑150 (MQ=255)
                      gTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGcc                                                                                                                 >  2:546085/1‑150 (MQ=255)
                      gTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGcc                                                                                                                 <  2:199026/150‑1 (MQ=255)
                                 gAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCggg                                                                                                      <  2:487382/150‑1 (MQ=255)
                                      gcgcAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCAGCGCCTGGCGCTCGAcca                                                                                                                                            >  1:404616/1‑107 (MQ=255)
                                      gcgcAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCAGCGCCTGGCGCTCGAcca                                                                                                                                            <  2:404616/107‑1 (MQ=255)
                                            cTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTc                                                                                           >  2:264471/1‑150 (MQ=255)
                                                 aGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCAGCGCCTGGCGCTCGAccagacat                                                                                                                                       <  1:15075/101‑6 (MQ=38)
                                                 aGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCAGCGCCTGGCGCTCGAccagacat                                                                                                                                       >  2:15075/1‑96 (MQ=38)
                                                                                              aGGGCACCAATGGCATCGGCAAATGTGCGTTCGAGCGCCAGGCGCTGAACATACATCAGAACTAGCATTTTCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACg                                         <  1:167101/150‑1 (MQ=255)
                                                                                               gggCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACgt                                        <  1:436867/150‑1 (MQ=255)
                                                                                                 gCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGACAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACgtgt                                      >  2:323928/1‑150 (MQ=255)
                                                                                                  cACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTc                                     >  2:217383/1‑150 (MQ=255)
                                                                                                  cACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTACTCGACGTGTc                                     <  1:523748/150‑1 (MQ=255)
                                                                                                                            tCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGc           <  2:401789/150‑1 (MQ=255)
                                                                                                                                     aGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGCAAAGCCAGt  <  2:53320/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                               
GGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGCAAAGCCAGT  >  NC_002947/633989‑634271

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: