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%	45.4 / NA	14	intergenic (+144/‑36)	mvaB/PP_mr44	hydroxymethylglutaryl‑CoA lyase/CrcZ ncRNA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (10/4);  total (10/4)

CATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATA  >  NC_002947/4018073‑4018296                                                                                                                                 |                                                                                                caTTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCaa                                                                                          <  2:76312/136‑1 (MQ=255)        gtgGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTgg                                                                                   >  2:372023/1‑136 (MQ=255)              gAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgcc                                                                             >  1:218438/1‑135 (MQ=255)              gAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgca                                                                             >  1:218199/1‑136 (MQ=255)                         tCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTatat                                                                  >  1:141226/1‑136 (MQ=255)                         tCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTatat                                                                  >  1:184601/1‑136 (MQ=255)                               acaTTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTgg                                                            >  1:145040/1‑136 (MQ=255)                                                                 cAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAAACCGATACCCCATtaaa                          >  2:94584/1‑136 (MQ=255)                                                                         cAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAgg                  >  1:289740/1‑136 (MQ=255)                                                                              tttCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAATGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAATCAACAGGTAACg             <  1:372023/136‑1 (MQ=255)                                                                                   ggCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGctct        <  2:158186/136‑1 (MQ=255)                                                                                      ggggTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGAc     <  2:141226/135‑1 (MQ=255)                                                                                         gtggtgTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATa  >  2:151128/1‑136 (MQ=255)                                                                                         gtggtgTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATa  >  2:6042/1‑136 (MQ=255)                                                                                                                                 |                                                                                                CATTGGTGTGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATA  >  NC_002947/4018073‑4018296

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