Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,537,079	C→T	25.8%	intergenic (+56/+47)	nadE → / ← PP_4870	NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,537,079	0	C	T	25.8%	54.9 / 13.6	31	intergenic (+56/+47)	nadE/PP_4870	NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base C (18/5);  new base T (3/5);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 7.42e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01

GCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCG  >  NC_002947/5536939‑5537222                                                                                                                                             |                                                                                                                                                 gCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑gg                                                                                                                                         <  2:392863/149‑1 (MQ=25) aacGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑gg                                                                                                                                         <  2:420479/147‑1 (MQ=25)   cGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGt                                                                                                                                       <  2:288619/149‑1 (MQ=255)     gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGa                                                                                                                                     <  5:234790/149‑1 (MQ=255)                                                           aaGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCTGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGc                                                                               >  5:329927/1‑149 (MQ=255)                                                                                                    gAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                      >  7:293510/1‑149 (MQ=255)                                                                                                    gAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                      >  3:174401/1‑149 (MQ=255)                                                                                                        gggTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGCCGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt                                  >  3:430088/1‑149 (MQ=255)                                                                                                           ttCACCCTCGTATTCGCGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATAGAGATGTGGACCGGGAAGGTTCTGAAGATGCTGTAg                                >  8:208436/1‑148 (MQ=255)                                                                                                           ttCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                               >  6:398464/1‑149 (MQ=255)                                                                                                            tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GTGTTAACCCGCTCCTGCAGGGTTCGCCGTCATCCTTTGATTACTTGACGATTACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAATAAGCt                                    <  5:264230/143‑1 (MQ=18)                                                                                                            tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg                                                  <  3:100732/129‑1 (MQ=18)                                                                                                            tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg                                                  >  4:100732/1‑129 (MQ=18)                                                                                                            tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                                    <  5:270122/143‑1 (MQ=21)                                                                                                            tCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCt                                    >  6:270122/1‑143 (MQ=21)                                                                                                            tCACCCGCGAATGCGTCGGTGGGTTCACAATCGTATTCGC‑GGGTGAACCAGCTCCTGCAGGGTTAGCCGTCAGCCTTTGATTTCTTGACGATGACCTTGCCTTTCATCATCGAGATGTGTCACGGGAAGGTGCAGAAGAAGATgatatc                              >  6:264230/1‑144 (MQ=18)                                                                                                                ccGCGAATACGAGGGTGAACCCTCCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                    <  6:44511/107‑1 (MQ=255)                                                                                                                ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCatcatc                                                                    >  5:44511/1‑107 (MQ=255)                                                                                                                ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                          >  3:300219/1‑149 (MQ=255)                                                                                                                ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                          >  1:230628/1‑149 (MQ=255)                                                                                                                ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                          >  7:201674/1‑149 (MQ=255)                                                                                                                ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCATTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                          >  5:34267/1‑149 (MQ=255)                                                                                                                ccGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGAAGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacca                          >  2:378736/1‑149 (MQ=255)                                                                                                                       tACGAGGGTGAACCCACCGACGCATTCGCAGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCt                    >  1:34540/1‑149 (MQ=255)                                                                                                                           tcggTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg               >  4:427812/3‑149 (MQ=255)                                                                                                                               tGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg           >  5:359063/1‑149 (MQ=255)                                                                                                                                gAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg           >  4:6415/1‑148 (MQ=255)                                                                                                                                  aCCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGa         <  6:34267/148‑1 (MQ=255)                                                                                                                                      aCCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgtgt    >  2:67754/1‑149 (MQ=255)                                                                                                                                       ccGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCCTGGAAGTGTc   <  6:329927/149‑1 (MQ=255)                                                                                                                                         gACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCg  >  2:301854/1‑148 (MQ=255)                                                                                                                                             |                                                                                                                                                 GCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGC‑GGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCG  >  NC_002947/5536939‑5537222

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