Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	498,522	Δ406 bp	coding (468‑873/963 nt)	hemH →	ferrochelatase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	498522	498927	406	36 [0]	[0] 36	hemH	ferrochelatase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 498521	0 (0.000)	36 (0.890)	34/482	0.3	100%	coding (467/963 nt)	hemH	ferrochelatase
	?	NC_000913	498928 =	0 (0.000)					coding (874/963 nt)	hemH	ferrochelatase

CCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGC  >  NC_000913/498277‑498762                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  ccgcTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcggg                                                                                                                                                                                                                                             >  2:157525‑M1/1‑245 (MQ=255)      gATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaa                                                                                                                                                                                                                                        <  1:104963‑M1/251‑12 (MQ=255)           tttACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCCTTCCGGGGATggcgggaaaaaatgtg                                                                                                                                                                                                                                   >  2:212626‑M1/1‑235 (MQ=255)                           aCAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgct                                                                                                                                                                                                                   <  1:157525‑M1/251‑33 (MQ=255)                             agcagGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgctt                                                                                                                                                                                                                  >  1:270090‑M1/1‑217 (MQ=255)                               cagGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaat                                                                                                                                                                                                               >  1:383091‑M1/1‑215 (MQ=255)                                       ggcccaacgtTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGTGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgcc                                                                                                                                                                                                       <  1:217154‑M1/247‑45 (MQ=255)                                          acaACGTTTACCGGAGATGCCCGTCGCGCTGGGCATGAGCTACGGCTCGCCATCCCTGGACCGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATTAACTGGCACGCATTCTGGCGCGCACACGTAGCATTCCGGGGATggcgggcaaaaatatgcatatattccggcgct                                                                                                                                                                                                                   >  2:192101‑M1/1‑204 (MQ=255)                                          acaACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGGGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgct                                                                                                                                                                                                                   <  1:192101‑M1/236‑33 (MQ=255)                                               gTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCCCGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatcttcctgcgcttaatgccacgccggaacat                                                                                                                                                                                                >  2:213151‑M1/1‑199 (MQ=255)                                               gTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTAGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggtaaaaatctgaatatattccggcgcttaatgccacgccggaccat                                                                                                                                                                                                >  2:213155‑M1/1‑199 (MQ=255)                                                   aCCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAATCGCCGTAGATGAACTCCTGGAATAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcg                                                                                                                                                                                            <  1:379802‑M1/250‑56 (MQ=255)                                                       gagaTGCCCGTAGCGCTGGGAAGGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatttgaatatattccggcgcttaatgccacgccggaacatatcgaaatg                                                                                                                                                                                       <  2:391635‑M1/251‑61 (MQ=255)                                                          aTGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatg                                                                                                                                                                                    >  1:358290‑M1/1‑188 (MQ=255)                                                                        gggAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgc                                                                                                                                                                      >  2:411734‑M1/1‑174 (MQ=255)                                                                                            cATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctg                                                                                                                                                    >  1:134097‑M1/1‑154 (MQ=255)                                                                                             aTCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgtgaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagc                                                                                                                                                 <  1:262090‑M1/251‑99 (MQ=255)                                                                                                 cTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggt                                                                                                                                              >  2:207969‑M1/1‑149 (MQ=255)                                                                                                            cGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagc                                                                                                                                  >  2:230554‑M1/1‑138 (MQ=255)                                                                                                                             cAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcct                                                                                                                 >  1:121697‑M1/1‑121 (MQ=255)                                                                                                                                   aTGTAGATCATAGTGGTGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTTCGGTATGGGATGAACTGGCACGCATTCTTGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatatgccggcgcttaatgccacgccggaacatatcgaaatgatggctaatctggttgcggcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcgg                                                                                                           <  1:251218‑M1/251‑137 (MQ=255)                                                                                                                                       aGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCCAACGTAGCATTCCGGGGATggcgggcaaaaatatgaatatattccggcgcttactgccacgccggaacatatcgacatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgccgcgaagagcctcgtcggcggg                                                                                                       >  2:299055‑M1/1‑111 (MQ=255)                                                                                                                                                tgtgGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcat                                                                                              >  1:383620‑M1/1‑102 (MQ=255)                                                                                                                                                  tggtgCTGCCGCTTTATCCGCAATTATCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggacaaaatatgcatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgccgcgaagagcctcgtcggcggttttcatcatcc                                                                                            >  2:347595‑M1/1‑100 (MQ=255)                                                                                                                                                   ggtgCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccg                                                                                           <  2:274148‑M1/251‑153 (MQ=255)                                                                                                                                                         gcCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggctttcatcatccgtgaat                                                                                      >  1:158763‑M1/1‑93 (MQ=255)                                                                                                                                                               ttATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgc                                                                                 <  1:86666‑M1/249‑163 (MQ=255)                                                                                                                                                                        aTTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtg                                                                      <  2:383620‑M1/251‑174 (MQ=255)                                                                                                                                                                                          cggtGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagag                                                    <  1:394566‑M1/251‑192 (MQ=255)                                                                                                                                                                                                        aTGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagg                                      <  2:13316‑M1/251‑206 (MQ=255)                                                                                                                                                                                                           aacgggcacgCATTCTTGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatatgccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataattcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggct                                   <  1:37752‑M1/247‑209 (MQ=255)                                                                                                                                                                                                           aaCTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgcatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgggcctggggagagtttgaactcacagggcc                                   >  2:281766‑M1/1‑43 (MQ=255)                                                                                                                                                                                                                           ggCGCGCAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgct                     <  1:411734‑M1/249‑223 (MQ=255)                                                                                                                                                                                                                                 cAAACGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctg             <  2:361980‑M1/251‑231 (MQ=255)                                                                                                                                                                                                                                   aaCGTAGCATTCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggt           <  2:221609‑M1/251‑233 (MQ=255)                                                                                                                                                                                                                                             tCCGGGGATggcgggaaaaaatatgaatatattccggcgcttaatgccacgccggaacatatcgaaatgatggctaatcttgttgccgcgtatcgctaaagctgagcggtaaagaactgagcgccgtcgcgaagagcctcgtcggcggttttcatcatccgtgaataatgcaaaaaggcgtgcagcgtgcctgggtagagtttgaactcacagggctgctgatgcgccgctaacgtctggtaaagcagac  <  2:383091‑M1/250‑242 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  CCGCTGATGGTTTACAGCCGCCAGCAACAGCAGGCGCTGGCACAACGTTTACCGGAGATGCCCGTAGCGCTGGGAATGAGCTACGGCTCGCCATCACTGGAAAGCGCCGTAGATGAACTCCTGGCAGAGCATGTAGATCATATTGTGGTGCTGCCGCTTTATCCGCAATTCTCCTGTTCTACGGTCGGTGCGGTATGGGATGAACTGGCACGCATTCTGGCGCGCAAACGTAGCATTCCGGGGATATCGTTTATACGTGATTACGCCGATAACCACGATTACATTAATGCACTGGCGAACAGCGTACGCGCTTCTTTTGCCAAACATGGCGAACCGGATCTGCTACTGCTCTCTTATCATGGCATTCCCCAGCGTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGC  >  NC_000913/498277‑498762

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