Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,264,293	T→G	9.8%	intergenic (‑47/+37)	zinU ← / ← PP_4640	P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,264,293	0	T	G	9.8%	112.4 / 4.8	41	intergenic (‑47/+37)	zinU/PP_4640	P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (19/18);  new base G (2/2);  total (21/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.76e-01

CACGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCG  >  NC_002947/5264145‑5264441                                                                                                                                                      |                                                                                                                                                     cACGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGt                                                                                                                                                      >  2:152242/1‑149 (MQ=255)   cGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTa                                                                                                                                                    <  2:85170/149‑1 (MQ=255)       gaccctCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCa                                                                                                                                                <  4:46226/144‑1 (MQ=255)       gAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCa                                                                                                                                                <  3:84642/149‑1 (MQ=255)           ttCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGATTGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCt                                                                                                                                            >  8:275409/1‑149 (MQ=255)                 ccATGCTCGGCCTTGTGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCca                                                                                                                                        >  8:71479/1‑147 (MQ=255)                  cttGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGCTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTAGGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAg                                                                                                                                     <  7:269576/147‑1 (MQ=255)                   aTGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAgg                                                                                                                                    >  1:165248/1‑149 (MQ=255)                      cTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGtt                                                                                                                                  <  7:275409/148‑1 (MQ=255)                      cTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGtt                                                                                                                                  >  4:103641/1‑148 (MQ=255)                       tCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACC‑GGCTCCCACAg                                                                                                                                     <  3:200912/144‑1 (MQ=255)                       tCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACC‑GGCTCCCACAg                                                                                                                                     >  4:205000/1‑144 (MQ=255)                       tCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACC‑GGCTCCCACAg                                                                                                                                     >  4:200912/1‑144 (MQ=255)                       tCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCCTGGGTAACC‑GGCTCCCACAg                                                                                                                                     <  3:205000/144‑1 (MQ=255)                           ccTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTacac                                                                                                                             >  5:269895/1‑148 (MQ=255)                             ttGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGaa                                                                                                                          <  6:55664/149‑1 (MQ=255)                                 gCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCg                                                                                                                       >  4:73127/1‑148 (MQ=255)                                   aTGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAg                                                                                                                    <  7:71479/149‑1 (MQ=255)                                    tGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTaca                                                                                                                              >  3:311269/1‑138 (MQ=255)                                    tGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTaca                                                                                                                              <  4:311269/138‑1 (MQ=255)                                        cTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCa                                                                                                               <  2:196673/149‑1 (MQ=255)                                         ttGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCaa                                                                                                              >  3:228417/1‑149 (MQ=255)                                                aGGTGGTCTTGCCTGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCATGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACt                                                                                                       <  5:323617/149‑1 (MQ=255)                                                                    aGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCg                                                                                   <  2:165248/149‑1 (MQ=255)                                                                           cTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTg                                                                            >  6:271027/1‑149 (MQ=255)                                                                             gtgAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGc                                                                                        >  1:156138/1‑135 (MQ=255)                                                                             gtgAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGc                                                                                        <  2:156138/135‑1 (MQ=255)                                                                                aGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTc                                                                        <  6:269895/148‑1 (MQ=255)                                                                                             gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCGAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCgg                                                           >  2:106600/1‑148 (MQ=255)                                                                                             gAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCgg                                                           >  7:2988/1‑148 (MQ=255)                                                                                              aaTGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGgaa                                                         <  5:271027/149‑1 (MQ=255)                                                                                                   gCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACCTTCGCGGAAGAATg                                                    >  5:322097/1‑149 (MQ=255)                                                                                                            cGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGc                                           <  3:73127/149‑1 (MQ=255)                                                                                                             gCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTg                                                                                         >  6:227062/1‑102 (MQ=255)                                                                                                             gCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTg                                                                                         <  5:227062/102‑1 (MQ=255)                                                                                                                        cTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCgg                                 <  3:103641/147‑1 (MQ=255)                                                                                                                         tGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCtt                              >  8:156569/1‑149 (MQ=255)                                                                                                                            aTCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTc                           <  4:228417/149‑1 (MQ=255)                                                                                                                                        cTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGttt                >  6:182986/1‑148 (MQ=255)                                                                                                                                          gtgGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCt             >  5:100499/1‑149 (MQ=255)                                                                                                                                             ggAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCa          >  6:323031/1‑149 (MQ=255)                                                                                                                                              gAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCAt         >  7:24629/1‑149 (MQ=255)                                                                                                                                                      ttACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCg  <  5:323031/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CACGGCGAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCG  >  NC_002947/5264145‑5264441

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: ‑