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%	199.4 / NA	55	G146G (GGA→GGG)	ccmH	heme lyase, CcmH subunit
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/26);  total (29/26)

TCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCAC  >  NC_000913/2291845‑2292099                                                                                                                                 |                                                                                                                               tCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCt                                                                                                                      >  2:833853/1‑139 (MQ=255)         cGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                            <  1:1178659/125‑1 (MQ=255)         cGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                            >  2:1178659/1‑125 (MQ=255)            cGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGc                                                                                                           <  2:210769/139‑1 (MQ=255)               aGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACg                                                                                                        <  1:833853/139‑1 (MQ=255)                   tgtgCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCtt                                                                                                                >  2:931461/1‑127 (MQ=255)                   tgtgCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCtt                                                                                                                <  1:931461/127‑1 (MQ=255)                           ggCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTTACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACACACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGa                                                                                            <  1:1587249/139‑1 (MQ=255)                             ccTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaaa                                                                                          >  2:534679/1‑139 (MQ=255)                              ctgctgCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAAc                                                                                         <  1:534679/139‑1 (MQ=255)                                       gATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGc                                                                                                                           <  1:638079/96‑1 (MQ=255)                                       gATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGc                                                                                                                           >  2:638079/1‑96 (MQ=255)                                                  gctgATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAc                                                                     >  2:906704/1‑139 (MQ=255)                                                   ctgATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAcc                                                                    >  2:1392769/1‑139 (MQ=255)                                                         aTTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGAcgc                                                              >  1:1401337/1‑139 (MQ=255)                                                             ccAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTAcc                                                                                                                  >  2:516568/1‑83 (MQ=255)                                                             ccAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTAcc                                                                                                                  <  1:516568/83‑1 (MQ=255)                                                               aGTCTGGTAGTAGCTGATGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTACAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                        <  1:1065856/139‑1 (MQ=255)                                                                 tCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaa                                                                                           >  1:1140216/1‑102 (MQ=255)                                                                 tCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaa                                                                                           <  2:1140216/102‑1 (MQ=255)                                                                           gCTGACGCCAGCCACAATTACCGCCACCACAATACCCGCCAGATAAATAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCTCCACGCGTACCCGCCGCCGCGAGCGGTCGTATATg                                            >  1:1276232/1‑139 (MQ=255)                                                                               aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGAccc                                        >  2:274348/1‑139 (MQ=255)                                                                               aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGAccc                                        >  2:1642319/1‑139 (MQ=255)                                                                                cGCCAGCCACAGTTAACGCCACCACAATACCCGGCAGATAAACAACATAcccag                                                                                                                            >  2:1004476/1‑52 (MQ=255)                                                                                cGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCgg                                                                                                                            <  1:1004476/54‑1 (MQ=255)                                                                                cGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                    <  1:728090/126‑1 (MQ=255)                                                                                cGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                    >  2:728090/1‑126 (MQ=255)                                                                                cGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCa                                       <  1:1642319/139‑1 (MQ=255)                                                                                 gccagccaCAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCaa                                      <  1:274348/139‑1 (MQ=255)                                                                                     gccaCAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCAGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgc                                  >  2:997106/1‑139 (MQ=255)                                                                                      ccaCAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgcc                                 >  2:1332073/1‑139 (MQ=255)                                                                                             tAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgcca                          <  1:1392769/139‑1 (MQ=255)                                                                                                    accacAATACCCGGCAGATAAACAACATACCCGGCa                                                                                                                          >  1:681428/1‑36 (MQ=255)                                                                                                    accacAATACCCGGCAGATAAACAACATACCCGGCa                                                                                                                          <  2:681428/36‑1 (MQ=255)                                                                                                    accacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                              <  1:1541100/80‑1 (MQ=255)                                                                                                    accacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcac                                                                              >  2:1541100/1‑80 (MQ=255)                                                                                                    accacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                        <  2:627383/102‑1 (MQ=255)                                                                                                    accacAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                        >  1:627383/1‑102 (MQ=255)                                                                                                      cacaATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCact                 >  1:99265/1‑139 (MQ=255)                                                                                                        caATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCactac               <  1:1332073/139‑1 (MQ=255)                                                                                                                ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                        <  1:1362607/90‑1 (MQ=255)                                                                                                                ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                        >  2:1362607/1‑90 (MQ=255)                                                                                                                ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCgg                                                     >  2:1731496/1‑93 (MQ=255)                                                                                                                ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCgg                                                     <  1:1731496/93‑1 (MQ=255)                                                                                                                  cAGATAAACAACATTCCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcaccac                                                                           <  2:353857/69‑1 (MQ=255)                                                                                                                  cAGATAAACAACATTCCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGcaccac                                                                           >  1:353857/1‑69 (MQ=255)                                                                                                                    gATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCca   >  2:1527051/1‑139 (MQ=255)                                                                                                                     aTAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCcac  >  1:1104005/1‑139 (MQ=255)                                                                                                                       aaaCAACATACCCGGTACGCTTACCTTCCGGCACGCTTTGt                                                                                                  >  2:1474759/1‑41 (MQ=255)                                                                                                                       aaaCAACATACCCGGTACGCTTACCTTCCGGCACGCTTTGt                                                                                                  <  1:1474759/41‑1 (MQ=255)                                                                                                                       aaaCAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                        >  1:830114/1‑83 (MQ=255)                                                                                                                       aaaCAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                        <  2:830114/83‑1 (MQ=255)                                                                                                                        aacaacATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                        >  2:5456/1‑82 (MQ=255)                                                                                                                        aacaacATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                        <  1:5456/82‑1 (MQ=255)                                                                                                                            acaTACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAccc                                                                   >  2:398698/1‑67 (MQ=255)                                                                                                                            acaTACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTAccc                                                                   <  1:398698/67‑1 (MQ=255)                                                                                                                                 |                                                                                                                               TCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCAC  >  NC_000913/2291845‑2292099

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