Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,053,227	C→T	H395Y (CAC→TAC)	hemN →	coproporphyrinogen III dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,053,227	0	C	T	100.0%	115.6 / NA	32	H395Y (CAC→TAC)	hemN	coproporphyrinogen III dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (16/16);  total (16/16)

ATTCGCCGCGATGTGATTAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGCACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGATGGGCTGGTGGATGTGGAT  >  NC_000913/4053149‑4053301                                                                               |                                                                           aTTCGCCGCGATGTGATTAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCTCCTATTGAGAAACAGTGGGATTTGTAc                                                                          >  2:798673/1‑81 (MQ=255)    cgccgcGATGTGATTAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACtt                                                                        <  2:798617/80‑1 (MQ=255)    cgccgcGATGTGATTAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACtt                                                                        <  1:798616/80‑1 (MQ=255)      ccgcGATGTGATTAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGc                                                                     <  2:798618/81‑1 (MQ=255)            tgtgATTAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTa                                                               <  1:798619/81‑1 (MQ=255)              tgATTAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTAc                                                              <  2:798620/80‑1 (MQ=255)              tgATTAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACt                                                             <  2:798621/81‑1 (MQ=255)                 ttAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTg                                                          <  1:798622/81‑1 (MQ=255)                         cTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCgg                                                       >  2:798664/1‑76 (MQ=255)                           cATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGAtct                                                >  1:798680/1‑81 (MQ=255)                               tGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAg                                            >  1:798676/1‑81 (MQ=255)                                 cAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAgc                                           <  1:798623/80‑1 (MQ=255)                                  aaCTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAgctg                                         <  1:798624/81‑1 (MQ=255)                                     ttCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTc                                      >  2:798667/1‑81 (MQ=255)                                       ccGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCg                                     <  1:798625/80‑1 (MQ=255)                                        cGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGcc                                   >  2:798668/1‑81 (MQ=255)                                          tCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGcccc                                 <  2:798626/81‑1 (MQ=255)                                              gATTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTa                             <  2:798627/81‑1 (MQ=255)                                               aTTACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAg                            >  1:798696/1‑81 (MQ=255)                                                ttACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAg                            >  1:798670/1‑80 (MQ=255)                                                ttACGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGc                           <  2:798628/81‑1 (MQ=255)                                                  aCGCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCa                          >  1:798684/1‑80 (MQ=255)                                                    gCCCCTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCaa                         >  1:798669/1‑79 (MQ=255)                                                      cccTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGa                     <  1:798629/81‑1 (MQ=255)                                                       ccTATTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGAt                    <  2:798630/81‑1 (MQ=255)                                                          aTTGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGATggg                 >  1:798675/1‑81 (MQ=255)                                                           ttGAGAAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGATGGGc                >  1:798674/1‑81 (MQ=255)                                                             gagaAACAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGATGGGCtg              >  2:798634/1‑81 (MQ=255)                                                                aaaCAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGATGGGCtggtg           >  1:798671/1‑81 (MQ=255)                                                                 aaCAGTGGGATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGATGGGCtggtg           >  2:798685/1‑80 (MQ=255)                                                                         gATTTGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGATGGGCTGgtggatgtggat  <  1:798631/81‑1 (MQ=255)                                                                            ttGTACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGATGGGCTGgtggatgtggat  >  2:798672/1‑78 (MQ=255)                                                                               |                                                                           ATTCGCCGCGATGTGATTAAGTCGCTCATCTGCAACTTCCGTCTGGATTACGCCCCTATTGAGAAACAGTGGGATTTGCACTTCGCTGATTACTTTGCGGAAGATCTCAAGCTGCTCGCCCCGTTAGCAAAAGATGGGCTGGTGGATGTGGAT  >  NC_000913/4053149‑4053301

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