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%	169.0 / NA	53	R320R (CGC→CGT)  *320* (TAG→TAA)	hemH aes	ferrochelatase acetyl esterase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/27);  total (26/27)

CTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGC  >  NC_000913/498883‑499145                                                                                                                                    |                                                                                                                                    cTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCt                                                                                                                              >  1:959291/1‑139 (MQ=255) cTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCt                                                                                                                              >  1:296973/1‑139 (MQ=255)  tGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCCCGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTg                                                                                                                             >  2:217389/1‑139 (MQ=255)    gAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                           <  1:1027700/139‑1 (MQ=255)            tGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAg                                                                                                                   >  1:1180928/1‑139 (MQ=255)                  gCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGa                                                                                                             >  1:1558497/1‑139 (MQ=255)                       aCCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGcc                                                                                                        <  1:217389/139‑1 (MQ=255)                           tGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                           >  2:693624/1‑116 (MQ=255)                           tGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                           <  1:693624/116‑1 (MQ=255)                              ggTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGa                                                                                                 <  2:965408/139‑1 (MQ=255)                                          tGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgtcg                                                                                     >  1:1122064/1‑139 (MQ=255)                                           gCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcgg                                                                                    <  2:1454067/139‑1 (MQ=255)                                            ccGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcggc                                                                                   >  2:1483360/1‑139 (MQ=255)                                             cggcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTcggcg                                                                                  <  2:296973/139‑1 (MQ=255)                                               gcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTcgcg                                                                                                  >  1:1673040/1‑121 (MQ=255)                                               gcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTcgcg                                                                                                  <  2:1673040/121‑1 (MQ=255)                                               gcggGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGt                                                                                >  2:1417253/1‑139 (MQ=255)                                                     aaaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          <  2:1558497/139‑1 (MQ=255)                                                     aaaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          <  2:1122064/139‑1 (MQ=255)                                                     aaaaaTATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          <  2:1001074/139‑1 (MQ=255)                                                             gAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTg                                                                  >  1:535142/1‑139 (MQ=255)                                                                        ggCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          <  2:131828/120‑1 (MQ=255)                                                                        ggCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          >  1:131828/1‑120 (MQ=255)                                                                         gcgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTa                                                                                                                      >  1:797563/1‑75 (MQ=255)                                                                         gcgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTa                                                                                                                      <  2:797563/75‑1 (MQ=255)                                                                          cgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                           >  1:709428/1‑69 (MQ=255)                                                                          cgcTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                           <  2:709428/69‑1 (MQ=255)                                                                              tAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCg                                                 <  2:1180928/139‑1 (MQ=255)                                                                                 tGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGc                                              >  2:1218305/1‑139 (MQ=255)                                                                                   ccACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAg                                            <  2:959291/139‑1 (MQ=255)                                                                                   ccACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAg                                            <  1:1483360/139‑1 (MQ=255)                                                                                       gCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          <  2:406404/105‑1 (MQ=255)                                                                                       gCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          >  1:406404/1‑105 (MQ=255)                                                                                           gAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                           <  2:560562/52‑1 (MQ=255)                                                                                           gAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAg                                                                                                                           >  1:560562/1‑52 (MQ=255)                                                                                           gAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTg                                    <  1:1218305/139‑1 (MQ=255)                                                                                                 atCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgt                                                                                       >  1:758622/1‑82 (MQ=255)                                                                                                 atCGAAATGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgt                                                                                       <  2:758622/82‑1 (MQ=255)                                                                                                      aaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgtcg                                                                                     >  1:358326/1‑79 (MQ=255)                                                                                                      aaTGATGGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgtcg                                                                                     <  2:358326/79‑1 (MQ=255)                                                                                                          atgGCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACt                     >  2:1454873/1‑139 (MQ=255)                                                                                                          atgGCTAATCTTGTAGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          <  2:150661/86‑1 (MQ=255)                                                                                                          atgGCTAATCTTGTAGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          >  1:150661/1‑86 (MQ=255)                                                                                                            ggCTAATCTTGTTGCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTca                   <  1:1417253/139‑1 (MQ=255)                                                                                                                      gttgCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTcaca                 >  2:1289457/1‑131 (MQ=255)                                                                                                                      gttgCCGCGTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTcaca                 <  1:1289457/131‑1 (MQ=255)                                                                                                                           cgcgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          <  2:1288180/69‑1 (MQ=255)                                                                                                                           cgcgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTtcat                                                                          >  1:1288180/1‑69 (MQ=255)                                                                                                                            gcgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAg                                            <  2:896322/98‑1 (MQ=255)                                                                                                                            gcgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAg                                            >  1:896322/1‑98 (MQ=255)                                                                                                                             cgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGt                                                                   >  1:511264/1‑74 (MQ=255)                                                                                                                             cgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGt                                                                   <  2:511264/74‑1 (MQ=255)                                                                                                                             cgTATCGTTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATgcgc  >  1:1230054/1‑139 (MQ=255)                                                                                                                                    |                                                                                                                                    CTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGC  >  NC_000913/498883‑499145

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