Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,291,973	T→C	G146G (GGA→GGG)	ccmH ←	heme lyase, CcmH subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,291,973	0	T	C	100.0%	124.5 / NA	35	G146G (GGA→GGG)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (16/19);  total (16/19)

CCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTG  >  NC_000913/2291841‑2292088                                                                                                                                     |                                                                                                                    ccTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCa                                                                                                               <  2:344713/139‑1 (MQ=255)               ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCgg                                                                                                 >  2:290508/1‑139 (MQ=255)               ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCga                                                                                                 >  2:765392/1‑138 (MQ=255)               ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTACTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                 >  1:1218429/1‑123 (MQ=255)               ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTACTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                 <  2:1218429/123‑1 (MQ=255)                cGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCGCAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                <  1:799698/139‑1 (MQ=255)                     ccTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGACACGCt                                                                                           <  1:765392/139‑1 (MQ=255)                                  ctgctgCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAAc                                                                              <  1:771815/139‑1 (MQ=255)                                            aTTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGca                                                                    <  1:362194/139‑1 (MQ=255)                                                      gctgATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAc                                                          >  2:878243/1‑139 (MQ=255)                                                          aTAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGc                                                                     >  2:490874/1‑124 (MQ=255)                                                          aTAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGc                                                                     <  1:490874/124‑1 (MQ=255)                                                                 ccAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCa                                                                                                               <  1:1191404/75‑1 (MQ=255)                                                                 ccAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCa                                                                                                               >  2:1191404/1‑75 (MQ=255)                                                                 ccAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCtt                                                                                                          <  1:1178975/80‑1 (MQ=255)                                                                 ccAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCtt                                                                                                          >  2:1178975/1‑80 (MQ=255)                                                                  cAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                           <  2:182538/94‑1 (MQ=255)                                                                  cAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                           >  1:182538/1‑94 (MQ=255)                                                                   aGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                             <  1:878243/139‑1 (MQ=255)                                                                    gTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                   <  1:804985/116‑1 (MQ=255)                                                                    gTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                   >  2:804985/1‑116 (MQ=255)                                                                         gtagtaGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                 >  2:1221625/1‑65 (MQ=255)                                                                         gtagtaGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                 <  1:1221625/65‑1 (MQ=255)                                                                               gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                             <  1:1162401/127‑1 (MQ=255)                                                                               gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                             >  2:1162401/1‑127 (MQ=255)                                                                               gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATg                                 >  1:841498/1‑139 (MQ=255)                                                                               gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATg                                 >  2:728689/1‑139 (MQ=255)                                                                                  gACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGAcc                              <  2:437718/139‑1 (MQ=255)                                                                                   aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGAccc                             >  2:793815/1‑139 (MQ=255)                                                                                       cagccaCAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                         >  1:501352/1‑123 (MQ=255)                                                                                       cagccaCAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                         <  2:501352/123‑1 (MQ=255)                                                                                                 tAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgcca               <  1:290508/139‑1 (MQ=255)                                                                                                              aTACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTg  <  1:99237/139‑1 (MQ=255)                                                                                                                       aGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGAcgccg                                                 <  2:386282/83‑1 (MQ=255)                                                                                                                       aGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGAcgccg                                                 >  1:386282/1‑83 (MQ=255)                                                                                                                                     |                                                                                                                    CCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTG  >  NC_000913/2291841‑2292088

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑