Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,018,201	G→A	intergenic (+144/‑36)	mvaB → / → PP_mr44	hydroxymethylglutaryl‑CoA lyase/CrcZ ncRNA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,018,201	0	G	A	100.0%	49.1 / NA	15	intergenic (+144/‑36)	mvaB/PP_mr44	hydroxymethylglutaryl‑CoA lyase/CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (8/7);  total (8/7)

TGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATT  >  NC_002947/4018079‑4018332                                                                                                                           |                                                                                                                                    tgtgGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCAGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTg                                                                                                                        <  1:280306/136‑1 (MQ=255)  gtgGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTgg                                                                                                                       >  1:27793/1‑136 (MQ=255)       tGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgc                                                                                                                  <  2:27793/136‑1 (MQ=255)       tGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacac                                                                                                                  <  1:332423/136‑3 (MQ=255)        gAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgca                                                                                                                 <  2:215813/136‑1 (MQ=255)         aaGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgcac                                                                                                                >  1:252174/1‑136 (MQ=255)             gtTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTg                                                                                                            >  2:65407/1‑136 (MQ=255)             gtTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTg                                                                                                            >  2:8186/1‑136 (MQ=255)              ttACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGc                                                                                                           <  1:8186/136‑1 (MQ=255)                                                              aTCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAaaca                                                           <  2:252174/136‑1 (MQ=255)                                                                   cAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAgg                                                      >  2:129381/1‑136 (MQ=255)                                                                                            aCCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGaaca                             >  1:211450/1‑136 (MQ=255)                                                                                            aCCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGaaca                             >  1:236026/1‑136 (MQ=255)                                                                                               aaaaaTCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAacaaca                          >  2:137227/1‑136 (MQ=255)                                                                                                                       aaGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGAtt  <  1:129381/136‑1 (MQ=255)                                                                                                                           |                                                                                                                                    TGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAGCAGATT  >  NC_002947/4018079‑4018332

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