Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 337,822 C→G 26.7% intergenic (‑72/+57) PP_0279 ← / ← PP_0280 hypothetical protein/amino acid ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947337,8220CG26.7% 49.2 / 19.2 30intergenic (‑72/+57)PP_0279/PP_0280hypothetical protein/amino acid ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (12/10);  new base G (6/2);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.78e-01

CTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACG  >  NC_002947/337676‑337969
                                                                                                                                                   |                                                                                                                                                   
ctgGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTc                                                                                                                                                     <  1:16788/147‑1 (MQ=255)
         gATCAGCAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTTTTTGGGCATGGCAATTATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGgc                                                                                                                                            <  1:120866/147‑1 (MQ=39)
             aGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTaca                                                                                                                                      <  6:175889/149‑1 (MQ=255)
              ggAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTacac                                                                                                                                     <  1:959/149‑1 (MQ=255)
                aaCTGGGCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCt                                                                                                                                   <  5:50533/149‑1 (MQ=255)
                  cttgTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGTCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                 <  5:377403/146‑1 (MQ=255)
                  cTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCtgt                                                                                                                                 <  4:193475/149‑1 (MQ=255)
                      tCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                       <  8:75267/139‑1 (MQ=255)
                      tCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTac                                                                                                                                       >  7:75267/1‑139 (MQ=255)
                              cTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCAGGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGt                                                                                                                     <  3:328092/149‑1 (MQ=255)
                                  aaGCGGGTGCTTAGGGGGGGGCTGTGGAGGTGGTTGGGCATGGCACGGCTGTGTGGGGGTGTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGccc                                                                                                                 <  8:154639/149‑1 (MQ=38)
                                                                                                                     gCCCGCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAg                              >  8:83742/1‑149 (MQ=21)
                                                                                                                     gCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAg                              >  5:114982/1‑149 (MQ=255)
                                                                                                                         gCTCCCACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCAt                           >  5:263767/1‑148 (MQ=21)
                                                                                                                             ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                       >  8:247107/1‑148 (MQ=25)
                                                                                                                             ccACAGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                       >  3:157372/1‑148 (MQ=25)
                                                                                                                             ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                       >  7:79304/1‑148 (MQ=255)
                                                                                                                             ccACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                       >  7:264129/1‑148 (MQ=255)
                                                                                                                              cacaGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                       >  2:61118/1‑147 (MQ=25)
                                                                                                                               acaGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc                                   >  4:376200/1‑134 (MQ=21)
                                                                                                                               acaGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGc                                   <  3:376200/134‑1 (MQ=21)
                                                                                                                               acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                       >  3:75491/1‑146 (MQ=255)
                                                                                                                               acaGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                       >  7:120851/1‑146 (MQ=255)
                                                                                                                                 aGGTGTAGCGCCAACCTTGAATGTGGCGCTATACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGAt                  <  4:157372/149‑1 (MQ=25)
                                                                                                                                 aGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                       >  3:338312/1‑144 (MQ=255)
                                                                                                                                 aGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAAc                       >  4:272894/1‑144 (MQ=255)
                                                                                                                                         cgcCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAgg           >  1:60749/1‑148 (MQ=255)
                                                                                                                                           ccACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcagtca                                                                                                     >  7:149462/1‑56 (MQ=255)
                                                                                                                                           ccACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGcagtca                                                                                                     <  8:149462/56‑1 (MQ=255)
                                                                                                                                                 ttCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAacac   >  4:1454/1‑149 (MQ=255)
                                                                                                                                                  tCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACg  >  1:238029/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTGGTGCTCGATCAGGAACTGGTCGAGGGCCTGGAAGCGGGTGCTTAGGTGGGGGCTGTGGAGGTGGTTGGGCATGGCAAGGATGTGTGGTGGTTTTGAGGGCCTCTTCGCGGGCACGCCCGCTCCCACAGGTATAGCGCCACATTCAAGGTTGGCGCTACACCTGTGGGAGCGGGCGTGCCCGCGAAGCAGTCAGTGAGGTCTCAGCGGCCTTGGCAGGCATCCACGCGCAACCAGCGCTCCAGCAGCTTGAAGCCTTGCACCAGCATGAACGAGATCACCAGGTAGAACACG  >  NC_002947/337676‑337969

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