Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,264,299	Δ1 bp	9.5%	intergenic (‑53/+31)	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,299	0	A	.	9.5%	124.1 / 7.9	42	intergenic (‑53/+31)	zinU/PP_4640	P‑loop guanosine triphosphatase‑dependent zinc chaperone/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/17);  new base . (2/2);  total (23/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.23e-01

GAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACATAG  >  NC_002947/5264151‑5264447                                                                                                                                                      |                                                                                                                                                     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                                                                                                                                        >  4:103641/1‑148 (MQ=255)                cTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGtt                                                                                                                                        <  7:275409/148‑1 (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                                                                                               <  5:227062/102‑1 (MQ=255)                                                                                                       gCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTg                                                                                               >  6:227062/1‑102 (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)                                                                                                                                                    ccAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACa     >  3:52253/1‑148 (MQ=255)                                                                                                                                                     cAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAAt                                                           <  7:136595/93‑1 (MQ=255)                                                                                                                                                     cAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAAt                                                           >  8:136595/1‑93 (MQ=255)                                                                                                                                                      aGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACATAg  <  5:182986/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GAGCTTCAGGCCATGCTCGGCCTTGAGCATGTACTTGAGCAAGGTGGTCTTGCCGGCACCGAGGAAGCCTGTGAGGACGGTTACGGGAATGGGCATCTGCACGCAAAATACTCCTGAATCTGAAAATGACTGTGGGAGCC‑GGTTACCCAGGCTCCCACAGGTTTACACGAACCGCAGGCTCAACAGCACTTGGGCCCGGACTTGCCACCGTAACGTGCTTCCTGACGTTCGCGGAAGAATGCCTCGTAGCTCATCACCGGCTTGTCCGGATGCGTTTTCTGCATGTGTTCGACATAG  >  NC_002947/5264151‑5264447

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