Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	499,203	A→G	intergenic (+47/‑249)	uhpA → / → PP_0411	two‑component system response regulator UhpA/polyamine ABC transporter ATP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	499,203	0	A	G	97.5%	129.0 / ‑4.4	40	intergenic (+47/‑249)	uhpA/PP_0411	two‑component system response regulator UhpA/polyamine ABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (20/19);  total (21/19)
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.89e-01

AAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCAGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTG  >  NC_002947/499058‑499326                                                                                                                                                  |                                                                                                                            aaGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGTAGCggg                                                                                                                           >  4:186472/1‑148 (MQ=255) aaGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCggg                                                                                                                           >  7:24963/1‑148 (MQ=255)   gACCATCAGCACCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCAGGc                                                                                                                          >  8:298284/1‑147 (MQ=255)           gcccccACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACACGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCcgcg                                                                                                                 >  5:269556/4‑148 (MQ=255)                cacaAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACggg                                                                                                           >  6:151188/1‑149 (MQ=255)                       cACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGggcctggcc                                                                                                    <  8:24963/149‑1 (MQ=255)                           ccTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCctgct                                                                                                <  6:6373/149‑1 (MQ=255)                             tGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGTGCATGCCCGCCAACGGGCCTGGCCTGCTGa                                                                                              >  5:384224/1‑149 (MQ=255)                             tGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGCACGGGCCAGGCCTGCTGa                                                                                              >  3:110044/1‑149 (MQ=255)                               aTGCAGATGTTGAAAGTTAACTGGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACATGCTGATCTGATTGCGACTTACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGAt                                                                                             <  2:50977/148‑1 (MQ=255)                               aTGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGAt                                                                                             <  5:151188/148‑1 (MQ=255)                                 gCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGAtct                                                                                           <  3:186472/148‑1 (MQ=255)                                     aaGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCAGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATg                                                                                      >  6:293812/1‑149 (MQ=255)                                       gTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGc                                                                                    <  4:236907/149‑1 (MQ=255)                                           aaaGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGgcctggcctg                                                                                                  <  1:320364/131‑1 (MQ=255)                                           aaaGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGgcctggcctg                                                                                                  >  2:320364/1‑131 (MQ=255)                                           aaaGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAg                                                                                >  8:319747/1‑149 (MQ=255)                                              gTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACATGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCg                                                                             <  7:319747/149‑1 (MQ=255)                                                  ccgcgcgGGCGGCGCGGGTGACGGACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCct                                                                         <  6:269556/146‑1 (MQ=255)                                                      gCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTc                                                                     >  7:257979/1‑149 (MQ=255)                                                          gCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTcgc                                                                   <  8:216389/147‑1 (MQ=255)                                                           cGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCgga                                                                <  5:293812/149‑2 (MQ=255)                                                               gCTGGNGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATg                                                            >  7:364730/1‑149 (MQ=255)                                                                cTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGc                                                           <  8:364730/149‑1 (MQ=255)                                                                     gAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGcc                                                                          <  8:260223/129‑1 (MQ=255)                                                                     gAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGcc                                                                          >  7:260223/1‑129 (MQ=255)                                                                               aTGGAGGACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGc                                                           <  8:241557/134‑1 (MQ=255)                                                                               aTGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGc                                                           >  7:241557/1‑134 (MQ=255)                                                                                 gtaGCACATTCTTATCTGTTTGCTACATAGGGTGGCCCTAAATTTCCCGAAGCCTTTGTGGGTGGGGGCTTGCCCGTGAAGGGGCCTGGCCGGCTGATCTCTATGGCTGAGCGGGCCTCTTCGGGGGAATGCCCGCTCCCACAGGttt                                           <  8:271529/146‑1 (MQ=255)                                                                                 ggAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATgg                                                                                     <  3:322714/106‑1 (MQ=255)                                                                                 ggAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATgg                                                                                     >  4:322714/1‑106 (MQ=255)                                                                                 ggAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGGTGAGCCGGCCGCTTCGCGGGCATGCCCGCGCCCACACGttt                                           >  7:271529/1‑148 (MQ=255)                                                                                        actctgatctgatTCCGACATACGGGTGCCCTAAATTTGCGTATGCCTATGTGTGAGCGTGAATGCCCGCTAACGGGCCTGGCCTGATTATCTCTTTGGCTGAGCGGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCTTGGTTg                                   <  6:384224/146‑1 (MQ=25)                                                                                        aaGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTg                                   >  5:75617/1‑149 (MQ=255)                                                                                                 tgatTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCCAGCCCGCGAACCGGCCTGGCCTGGGGATCTATATGGCTGAGCCGGCCAATTAGCCGGCCTGCCCCCTACCACACGTTTAGTGGTTTGCTCGCa                            >  4:135671/1‑147 (MQ=255)                                                                                                     tGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGAtct                                                                                           >  2:324955/1‑80 (MQ=255)                                                                                                     tGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGAtct                                                                                           <  1:324955/80‑1 (MQ=255)                                                                                                     tGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCCCCAACGGGCCGGGCCGGCAGTACTCAATGGCTGAGCCCGCCGCTTCGCGGGCCAGGCCCCACCCCCAGGTTTAGGGGATGGCGCGCAAgggtg                      >  8:275050/1‑147 (MQ=255)                                                                                                        gACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTc                   >  7:256138/1‑149 (MQ=255)                                                                                                             acaggTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGCACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTAGTGGTTGGCTCGCAAGGGCATTCCAATc              >  5:86553/4‑149 (MQ=255)                                                                                                                     ccGAAACTTGCCGAAGCCTATGTGGGAGCGGGCATGCCCGCGAACGGGCCTGGCCTGCGGATGTCTATGGGTGAGCCGGGCTCTTCGCGGGCCTGCCCCCTCCCCCAGGTTTCGTGGGTGGGGCGCAAGGGCAGTGCCATCCCGCCCTc      >  5:348775/1‑149 (MQ=255)                                                                                                                         aaCTTGCCGAAGCTTTTGTGGGAGGGGGCATGCCCGCTAACGGGCCTGGCCTGCTTATCTCTATGGCTGAGCGGGCCTCTTCGCGGGCTTGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTg  <  2:131992/149‑1 (MQ=38)                                                                                                                                                  |                                                                                                                            AAGACCATCAGCACCCACAAGGCACGCCTGATGCAGAAGTTGAAAGTGAACTCGCTGGCGGAGCTGGTGAAGTACGCCATGGAGCACAAGCTGATCTGATTGCGACATACGGGTGCCCGAAACTTGCCGAAGCCTATGTGGGAGCAGGCATGCCCGCGAACGGGCCTGGCCTGCTGATCTCTATGGCTGAGCCGGCCTCTTCGCGGGCATGCCCGCTCCCACAGGTTTCGTGGTTGGCTCGCAAGGGCATTCCAATCCCGCCATCCTTG  >  NC_002947/499058‑499326

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