Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,052,834	C→T	Q264* (CAG→TAG)	hemN →	coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,052,834	0	C	T	100.0%	215.5 / NA	55	Q264* (CAG→TAG)	hemN	coproporphyrinogen III oxidase, SAM and NAD(P)H dependent, oxygen‑independent
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/24);  total (31/24)

AAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTCC  >  NC_000913/4052690‑4052975                                                                                                                                                 |                                                                                                                                              aaaCAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCaaa                                                                                                                                         <  1:219628/151‑1 (MQ=255)    cAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAAc                                                                                                                                      >  1:305089/1‑151 (MQ=255)     aGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACt                                                                                                                                     >  1:636197/1‑151 (MQ=255)     aGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACt                                                                                                                                     <  2:474787/151‑1 (MQ=255)      gACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACt                                                                                                                                     <  1:282299/150‑1 (MQ=255)          ccGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGata                                                                                                                                <  1:565009/151‑1 (MQ=255)           cGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGatat                                                                                                                               <  1:110753/151‑1 (MQ=255)              agagTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCt                                                                                                                            >  2:234829/1‑151 (MQ=255)                  tttCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGc                                                                                                                          >  1:650304/1‑149 (MQ=255)                     cGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGaa                                                                                                                     <  1:95358/151‑1 (MQ=255)                     cGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGaa                                                                                                                     <  1:59941/151‑1 (MQ=255)                       ccTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGAACTGCCGAGTCCGTAGCAAAAACTAGATATCCTGAAGGaa                                                                                                                     >  1:44922/1‑149 (MQ=255)                           tACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATc                                                                                                               >  1:298662/1‑151 (MQ=255)                             cccTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGGCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAAcc                                                                                                                  <  1:538395/146‑1 (MQ=255)                               cTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCt                                                                                                           >  2:448004/1‑151 (MQ=255)                               cTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCt                                                                                                           >  1:413732/1‑151 (MQ=255)                               cTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCt                                                                                                           >  2:138247/1‑151 (MQ=255)                                tGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCtt                                                                                                          >  2:591304/1‑151 (MQ=255)                                     cGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGa                                                                                                     <  1:218527/151‑1 (MQ=255)                                      gtgtGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGAc                                                                                                    >  1:108934/1‑151 (MQ=255)                                           gCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAAt                                                                                               >  1:336407/1‑151 (MQ=255)                                             ggAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCg                                                                                             >  2:406083/1‑151 (MQ=255)                                              gAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCgg                                                                                            <  1:572396/151‑1 (MQ=255)                                               aGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCggg                                                                                           <  1:411297/151‑1 (MQ=255)                                                  tGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCt                                                                                         >  2:85579/1‑150 (MQ=255)                                                    aaCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATc                                                                                      <  2:302258/151‑1 (MQ=255)                                                    aaCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATc                                                                                      >  1:109518/1‑151 (MQ=255)                                                     aCCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATc                                                                                      >  1:551173/1‑150 (MQ=255)                                                        ccGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGtt                                                                                  >  2:599785/1‑151 (MQ=255)                                                                      gtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGAt                                                                     <  1:495287/150‑1 (MQ=255)                                                                      gtCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATc                                                                    <  1:26487/151‑1 (MQ=255)                                                                       tCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGCCCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGAt                                                                     <  2:287294/149‑1 (MQ=255)                                                                        cTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCAc                                                                  >  2:251721/1‑151 (MQ=255)                                                                          ttAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCa                                                                   <  1:50541/148‑1 (MQ=255)                                                                            aaCTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACttt                                                               <  2:447476/150‑1 (MQ=255)                                                                             aCTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGc                                                             <  1:515826/151‑1 (MQ=255)                                                                               tACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGccc                                                           >  1:451641/1‑151 (MQ=255)                                                                                     cATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCgg                                                     <  2:20223/151‑1 (MQ=255)                                                                                           ccGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACg                                               >  2:409140/1‑151 (MQ=255)                                                                                            cGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACg                                               >  2:202294/1‑150 (MQ=255)                                                                                             gACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACg                                               >  2:255260/1‑149 (MQ=255)                                                                                              aCCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGc                                            >  1:49653/1‑151 (MQ=255)                                                                                                cATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTg                                          >  2:254871/1‑151 (MQ=255)                                                                                                cATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTg                                          >  1:392053/1‑151 (MQ=255)                                                                                                    tttGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCgg                                      >  1:580637/1‑151 (MQ=255)                                                                                                       gctgctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTgg                                   >  2:99229/1‑151 (MQ=255)                                                                                                       gctgctCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTg                                    <  1:56719/150‑1 (MQ=255)                                                                                                              aGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGc                             <  2:575221/150‑1 (MQ=255)                                                                                                                          aaGATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtgct                <  1:264190/151‑1 (MQ=255)                                                                                                                            gATGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGtgctgc              >  1:331898/1‑151 (MQ=255)                                                                                                                             aTGCTGACCTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCa             <  1:528683/151‑1 (MQ=255)                                                                                                                                     cTGCCGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACt     >  2:443046/1‑151 (MQ=255)                                                                                                                                        ccGAGTCCGTAGCAAAAACTCGATATCCTGCTGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTcc  <  1:338707/151‑1 (MQ=255)                                                                                                                                         cGAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTcc  >  2:79137/1‑150 (MQ=255)                                                                                                                                          gAGTCCGTAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTcc  >  2:389205/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                              AAACAGACGCCGGAGAGTTTCGCCTTTACCCTGAAACGTGTGGCGGAGCTGAACCCCGATCGTCTGAGTGTCTTTAACTACGCGCATCTGCCGACCATTTTTGCTGCTCAGCGCAAAATCAAAGATGCTGACCTGCCGAGTCCGCAGCAAAAACTCGATATCCTGCAGGAAACCATCGCCTTCCTGACGCAATCGGGCTATCAGTTTATCGGTATGGATCACTTTGCCCGTCCGGATGACGAGCTGGCGGTGGCCCAGCGTGAAGGCGTGCTGCATCGTAACTTCC  >  NC_000913/4052690‑4052975

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