Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 358,726 T→C 21.0% W64R (TGG→CGG) ‡ gbdR → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947358,7260TC21.0% 28.7 / 4.9 20W64R (TGG→CGG) ‡gbdRAraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base T (9/6);  new base C (2/2);  total (12/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

AACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCTGTGGTGGTGTTGGCATCCAGCGTTCCGTCACCCGTGAACACGTCACCTGGTTGCAAGCCCA  >  NC_002947/358578‑358862
                                                                                                                                                    |                                                                                                                                        
aacaCCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGtgtg                                                                                                                                         <  2:140167/150‑1 (MQ=255)
       cAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGc                                                                                                                                  <  1:367007/150‑1 (MQ=255)
             aTCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATgac                                                                                                                            <  1:70407/150‑3 (MQ=255)
                     tcttctATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAAt                                                                                                                    >  2:217750/1‑150 (MQ=255)
                              ggACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCAAGAGCGTGGTGTCCGGCGATGGCCTGCAAATCACACCGGa                                                                                                           >  2:35801/1‑150 (MQ=255)
                               gACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGCAGGCCATCGCTGGcccaca                                                                                                                       <  2:37341/137‑5 (MQ=255)
                               gACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGCAGGCCATCGCTGGcccaca                                                                                                                       >  1:37341/1‑133 (MQ=255)
                                          ccTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAg                                                                                               >  1:330616/1‑150 (MQ=255)
                                           cTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGCAGGCCATCGCTGGcccacac                                                                                                                      <  2:303564/126‑6 (MQ=255)
                                           cTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCAGTCGGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGCAGGCCATCGCTGGcccacac                                                                                                                      >  1:303564/1‑121 (MQ=255)
                                            tGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTCAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGCTGGCCTGCACATCACACCGGATGCCGCAATCAGCa                                                                                             >  2:77899/1‑150 (MQ=255)
                                                                   gAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCAGATACTGCAATCAGCAATGCTCCGCCCATCGACACCGTg                                                                      <  1:35801/150‑1 (MQ=255)
                                                                         cTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTc                                                                >  2:443402/1‑150 (MQ=255)
                                                                          tGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCt                                                               >  2:357990/1‑150 (MQ=255)
                                                                                                  gCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCTGTGGTGGTGTTGGCATCCAGCGtt                                       <  1:443402/150‑1 (MQ=255)
                                                                                                           tGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCTGTGGTGGTGTTGGCATCCAGCGTTCCGTCAccc                              >  2:372275/1‑150 (MQ=255)
                                                                                                                       acacGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCTGTGGTGGTGTTGGCATCCAGCGTTCCGTCACCCGTGAACACGTCa                  >  1:263653/1‑150 (MQ=255)
                                                                                                                             tGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCTGTGGTGGTGTTGGCATCCAGCGTTCCGTCACCCGTGAACACGTCACCTGGt            >  1:38233/1‑150 (MQ=255)
                                                                                                                                 gatcgatGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCTGTGGTGGTGTTGGCATCCAGCGTTCCGTCACCCGTGAACACGTCACCTGGTTGCa        >  1:18542/1‑150 (MQ=255)
                                                                                                                                  atcgatGGCGGCCAGGTGTTGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCTGTGGTGGTGTTGGCATCCAGCGTTCCGTCACCCGTGAACACGTCACCTGGTTGCaa       <  2:263653/150‑1 (MQ=255)
                                                                                                                                       tGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCTGTGGTGGTGTTGGCATCCAGCGTTCCGTCACCCGTGAACACGTCACCTGGTTGCAAGCCCa  <  2:368958/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
AACACCCCAGTCCATCGGGTTTCTTCTATTGGACAATTTCACCCTGATTTCCCTGGCATCGGCGGTCGAGCCGCTGCGCATGGCCAACCAGCTGTCTGGCCGTGAGCTGTATCGCTGGCACACGCTGACGATCGATGGCGGCCAGGTGTGGGCCAGCGATGGCCTGCAAATCACACCGGATGCCGCAATCAGCAATGCTCCGCCCATCGACACCGTGATTGTCTGTGGTGGTGTTGGCATCCAGCGTTCCGTCACCCGTGAACACGTCACCTGGTTGCAAGCCCA  >  NC_002947/358578‑358862

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: