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

TGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAG  >  NC_002947/4018081‑4018326                                                                                                                         |                                                                                                                               tgGTTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGTCAGAATTTCAGTCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgc                                                                                                                 <  1:211907/136‑1 (MQ=255) tgGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgc                                                                                                                 <  2:528254/136‑1 (MQ=255)   ggTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcac                                                                                                               <  2:722264/136‑1 (MQ=255)    gTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacg                                                                                                              >  1:638035/1‑136 (MQ=255)     tGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgc                                                                                                             <  1:550993/136‑1 (MQ=255)     tGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgc                                                                                                             <  1:685556/136‑1 (MQ=255)      gAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGgcacgca                                                                                                            <  1:822161/136‑1 (MQ=255)           gtTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTg                                                                                                       >  1:682549/1‑136 (MQ=255)                       acaTTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTgg                                                                                           >  2:693510/1‑136 (MQ=255)                       acaTTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTgg                                                                                           <  2:46455/136‑1 (MQ=255)                         aTTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATATTCAAAAGTTGGCACGCACCCAGCTATATCTCAGGTa                                                                                         >  2:110501/1‑136 (MQ=255)                                    aaaaaTCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑aaaa                                                                             <  2:594292/136‑1 (MQ=255)                                                          ataTCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATaaaaa                                                       <  2:378798/134‑1 (MQ=255)                                                                 cAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAgg                                                <  2:591333/136‑1 (MQ=255)                                                                 aaGATTTTCAGGCTGGTTGTTTTCGACCAAAAATCATAACTAATTAATTTTAAAGTATTTTCAAAAGTTGGCACGCACCATGATATATCTCTGGTACAACAAC‑AAAAAAAAGGCGATACCCAATAAAAACAACAgt                                                <  1:541694/135‑2 (MQ=255)                                                                            gCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAACAAAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGctct                                      <  1:218421/136‑1 (MQ=255)                                                                                   ggtgTGGAACAAAAATCTTTACTCATTTATTTTAAAGTATTTTCAAAAGTTGGCAGGCACCCTGCTATATCTCTTGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAAc                              <  1:910088/136‑1 (MQ=255)                                                                                        cGACCAAAAATCTTAACTCATTAATATTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGaa                         <  1:709902/136‑1 (MQ=255)                                                                                             aaaaaTCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAacacca                    >  1:789716/1‑136 (MQ=255)                                                                                               aaaTCTTAACTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACg                  <  1:693510/136‑1 (MQ=255)                                                                                                        cTCATTAATTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCg         >  2:453429/1‑136 (MQ=255)                                                                                                              aaTTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTagcaa   >  2:371763/1‑136 (MQ=255)                                                                                                               aTTTTAAAGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTagcaag  >  2:808886/1‑136 (MQ=255)                                                                                                                      aGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGt        >  2:950894/1‑123 (MQ=255)                                                                                                                      aGGATTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGt        <  1:950894/123‑1 (MQ=255)                                                                                                                         |                                                                                                                               TGGGTGAAGTGTTACCTCCCGCACATTTCACAGCAAAAAATCGGGTAACAGGGAAACAAATCGGCAGAATTTCAGGCTGGGTGGTGTCGACCAAAAATCTTAACTCATTAATTTTAAAGGGTTTTCAAAAGTTGGCACGCACCCTGCTATATCTCTGGTACAACAAC‑AAAAAAAACGCGATACCCAATAAAAACAACAGGTAACGGCTCTGACATAACAAGAACAACACGGCAGAGGCGTAGCAAG  >  NC_002947/4018081‑4018326

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