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

GCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGC  >  NC_000913/2291849‑2292102                                                                                                                             |                                                                                                                                  gCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTAcc                                                                                                                     <  1:612318/139‑1 (MQ=255)       ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCtt                                                                                                                   >  2:370065/1‑135 (MQ=255)       ccGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCtt                                                                                                                   <  1:370065/135‑1 (MQ=255)                     gTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAg                                                                                                 >  1:951946/1‑139 (MQ=255)                      tGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTcggg                                                                                                >  2:568001/1‑136 (MQ=255)                      tGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAgg                                                                                                >  1:1589748/1‑139 (MQ=255)                       ggCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGa                                                                                               <  2:951946/139‑1 (MQ=255)                          ctgctgCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCGGGAAAc                                                                                            <  1:568001/139‑1 (MQ=255)                                          aCCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCa                                                                                                                             <  1:117484/90‑1 (MQ=255)                                          aCCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCa                                                                                                                             >  2:117484/1‑90 (MQ=255)                                          aCCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                         <  2:1278693/110‑1 (MQ=255)                                          aCCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                         >  1:1278693/1‑110 (MQ=255)                                               ctgATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATAccc                                                                                                                                 <  1:1067146/81‑1 (MQ=255)                                                      ttGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAAc                                                                                            >  2:1293054/1‑111 (MQ=255)                                                      ttGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAAc                                                                                            <  1:1293054/111‑1 (MQ=255)                                                         ccAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCcgcg                                                             >  2:929954/1‑139 (MQ=255)                                                         ccAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCcgcg                                                             >  1:1000975/1‑139 (MQ=255)                                                         ccAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCcgcg                                                             >  2:238785/1‑139 (MQ=255)                                                          cAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGa                                                            <  1:359390/139‑1 (MQ=255)                                                           aGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                           <  1:238785/139‑1 (MQ=255)                                                           aGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                           <  2:1000975/139‑1 (MQ=255)                                                            gTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACACTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGc                                                          >  1:1510098/1‑139 (MQ=255)                                                                 gtagtaGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCg                                                     >  1:1642021/1‑139 (MQ=255)                                                                  tagtagCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAAc                                                                                            >  1:1595604/1‑99 (MQ=255)                                                                  tagtagCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAAc                                                                                            <  2:1595604/99‑1 (MQ=255)                                                                       gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                         <  1:491926/81‑1 (MQ=255)                                                                       gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                         >  2:491926/1‑81 (MQ=255)                                                                       gCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATg                                               >  2:514668/1‑139 (MQ=255)                                                                           aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACg                                                                                                           >  2:1629002/1‑75 (MQ=255)                                                                           aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACg                                                                                                           <  1:1629002/75‑1 (MQ=255)                                                                           aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                       >  1:1215792/1‑127 (MQ=255)                                                                           aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCggg                                                       <  2:1215792/127‑1 (MQ=255)                                                                           aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGta                                                   <  1:552702/131‑1 (MQ=255)                                                                           aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGta                                                   >  2:552702/1‑131 (MQ=255)                                                                           aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGAccc                                           >  1:893440/1‑139 (MQ=255)                                                                           aCGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGAccc                                           >  2:1274670/1‑139 (MQ=255)                                                                                   cacaATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgcca                                   <  2:1642021/139‑1 (MQ=255)                                                                                   cacaATTAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgcca                                   <  1:1653733/139‑1 (MQ=255)                                                                                        ttAACGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgcc                              <  2:1589748/139‑1 (MQ=255)                                                                                           aCGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCg                                                     <  1:1564057/113‑1 (MQ=255)                                                                                           aCGCCACCACAATACCCGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCg                                                     >  2:1564057/1‑113 (MQ=255)                                                                                                         cccGGCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGgca             >  2:1553959/1‑139 (MQ=255)                                                                                                            ggCAGATAAACAACATACGCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgccaa                            >  2:1099652/1‑121 (MQ=255)                                                                                                            ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCtt                                                                                                        <  1:827995/45‑1 (MQ=255)                                                                                                            ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCtt                                                                                                        >  2:827995/1‑45 (MQ=255)                                                                                                            ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCTCCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCgccaatgccaa                            <  1:1099652/121‑1 (MQ=255)                                                                                                            ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                           >  1:625932/1‑90 (MQ=255)                                                                                                            ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAg                                                           <  2:625932/90‑1 (MQ=255)                                                                                                            ggCAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGgtagca          >  1:690631/1‑139 (MQ=255)                                                                                                              cAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaa                                                                                              >  2:1006096/1‑53 (MQ=255)                                                                                                              cAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGaa                                                                                              <  1:1006096/53‑1 (MQ=255)                                                                                                              cAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCAcc        <  2:893440/139‑1 (MQ=255)                                                                                                              cAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCAcc        >  1:927795/1‑139 (MQ=255)                                                                                                              cAGATAAACAACATACCCGGCACGCTTACCTTCCGGCACACTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCAcc        <  2:1510098/139‑1 (MQ=255)                                                                                                               aGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATg                                               >  2:207894/1‑99 (MQ=255)                                                                                                               aGATAAACAACATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATg                                               <  1:207894/99‑1 (MQ=255)                                                                                                                 aTAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                         >  1:1311987/1‑39 (MQ=255)                                                                                                                 aTAAACAACATACCCGGCACGCTTACCTTCCGGCACGCt                                                                                                         <  2:1311987/39‑1 (MQ=255)                                                                                                                    aacaacATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCAcagc  <  1:1553959/139‑1 (MQ=255)                                                                                                                    aacaacATACCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCAcagc  >  2:1507695/1‑139 (MQ=255)                                                                                                                             |                                                                                                                                  GCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAGTAGCTGACGCCAGCCACAATTAACGCCACCACAATACCCGGCAGATAAACAACATATCCGGCACGCTTACCTTCCGGCACGCTTTGTTCAGGAAACGCTTCCGGCACCACGCGTACCCGACGCCGCGAGCGGGCGTATATGACCCAACCGCCAATGCCAATAGCCACTACTGGCAGCACCCACAGC  >  NC_000913/2291849‑2292102

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