Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 634,115 Δ1 bp 21.1% coding (418/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,1150G.21.1% 43.7 / 9.6 19coding (418/2016 nt)PP_0546sigma‑54 dependent transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (5/10);  new base . (2/2);  total (7/12)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

GGTGCTGCTGACCGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGC  >  NC_002947/633967‑634262
                                                                                                                                                    |                                                                                                                                                    
ggTGCTGCTGACCGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGa                                                                                                                                                     <  2:26765/150‑1 (MQ=255)
  tgctgcTGACCGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAgc                                                                                                                                                   <  2:516343/150‑1 (MQ=255)
                      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)
                                                                                                                                                    |                                                                                                                                                    
GGTGCTGCTGACCGACGCCCGCGGGGTGATCCTCAATTGCGTCAGTGCCCCCGCCGAGCGGCGCAGCTTCGAGCGTGCCGGGCTGTGGCTGGGCGCTGACTGGAGCGAGGCGCGCGAGGGCACCAATGGCATCGGCACATGTCTGGTCGAGCGCCAGGCGCTGACCATCCATCAGAACGAGCACTTCCGTGGCCGTCACACCGGGCTGACCTGCTCGGCCAGCCCGGTGTTCGACCCGCATGGCGACCTGCTGGCGGTGCTCGACGTGTCCTCTGCGCGCCCCGACGTATCGCGGC  >  NC_002947/633967‑634262

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