Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	499,014	C→T	R320R (CGC→CGT)  *320* (TAG→TAA)	hemH → aes ←	ferrochelatase acetyl esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	499,014	0	C	T	100.0%	137.4 / NA	43	R320R (CGC→CGT)  *320* (TAG→TAA)	hemH aes	ferrochelatase acetyl esterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/20);  total (23/20)

AGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCG  >  NC_000913/498878‑499144                                                                                                                                         |                                                                                                                                   agaCGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTa                                                                                                                                  >  1:426569/1‑139 (MQ=255)     gCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGc                                                                                                                              >  2:55995/1‑139 (MQ=255)                          aaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGTCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAACGCTGAGCGGTAAAGAACTGAgcg                                                                                                         <  1:1183020/139‑1 (MQ=255)                          aaaCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAgcg                                                                                                         <  2:426569/139‑1 (MQ=255)                                 gAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTcgc                                                                                                  >  2:422173/1‑139 (MQ=255)                                                         aaaaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcgg                                                                                   <  1:860422/130‑1 (MQ=255)                                                         aaaaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcgg                                                                                   >  2:860422/1‑130 (MQ=255)                                                          aaaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         <  1:55995/139‑1 (MQ=255)                                                          aaaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         <  1:422173/139‑1 (MQ=255)                                                                  gAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         >  1:1636022/1‑131 (MQ=255)                                                                  gAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         <  2:1636022/131‑1 (MQ=255)                                                                   aaTATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAg                                                                                                           <  2:1241794/96‑1 (MQ=255)                                                                   aaTATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAg                                                                                                           >  1:1241794/1‑96 (MQ=255)                                                                        atTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaataat                                                           >  2:1789492/1‑139 (MQ=255)                                                                          tCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGc                                                         >  1:1745623/1‑139 (MQ=255)                                                                            cGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCaa                                                       >  2:674204/1‑139 (MQ=255)                                                                             ggCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCaaa                                                      <  1:674204/139‑1 (MQ=255)                                                                             ggCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCaaa                                                      <  1:1789492/139‑1 (MQ=255)                                                                              gcgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCaaaa                                                     >  2:1958589/1‑139 (MQ=255)                                                                              gcgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCaaaa                                                     >  2:1208477/1‑139 (MQ=255)                                                                                 cTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         <  1:744220/116‑1 (MQ=255)                                                                                 cTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         >  2:744220/1‑116 (MQ=255)                                                                                            gCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgt                                                                                      >  2:2052696/1‑92 (MQ=255)                                                                                            gCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgt                                                                                      <  1:2052696/92‑1 (MQ=255)                                                                                                      atCGAAATGATAGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTaga                             >  1:1596343/1‑139 (MQ=255)                                                                                                        cGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAg                                           >  1:1745842/1‑123 (MQ=255)                                                                                                        cGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAg                                           <  2:1745842/123‑1 (MQ=255)                                                                                                          aaaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTcgcg                                                                                                 >  2:890977/1‑67 (MQ=255)                                                                                                          aaaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTcgcg                                                                                                 <  1:890977/67‑1 (MQ=255)                                                                                                           aaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcgg                                                                                   >  2:1844521/1‑80 (MQ=255)                                                                                                           aaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcgg                                                                                   <  1:1844521/80‑1 (MQ=255)                                                                                                                tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         <  2:4192/85‑1 (MQ=255)                                                                                                                tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         >  1:4192/1‑85 (MQ=255)                                                                                                                tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaa                                                               >  2:1951730/1‑95 (MQ=255)                                                                                                                tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaa                                                               <  1:1951730/95‑1 (MQ=255)                                                                                                                tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTg                                       >  2:1816390/1‑119 (MQ=255)                                                                                                                tgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTg                                       <  1:1816390/119‑1 (MQ=255)                                                                                                                         ttgttgCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         >  2:1627269/1‑76 (MQ=255)                                                                                                                         ttgttgCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                         <  1:1627269/76‑1 (MQ=255)                                                                                                                              gCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGAt     <  1:1208477/139‑1 (MQ=255)                                                                                                                                 gcgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATgcg  >  1:549348/1‑139 (MQ=255)                                                                                                                                  cgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTc                                                                                                    <  2:1551648/40‑1 (MQ=255)                                                                                                                                  cgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTc                                                                                                    >  1:1551648/1‑40 (MQ=255)                                                                                                                                         |                                                                                                                                   AGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCG  >  NC_000913/498878‑499144

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