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	50 [0]	[0] 50	hemH	ferrochelatase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 498521	0 (0.000)	50 (0.590)	40/598	0.8	100%	coding (467/963 nt)	hemH	ferrochelatase
	?	NC_000913	498928 =	0 (0.000)					coding (874/963 nt)	hemH	ferrochelatase

GTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAG  >  NC_000913/498650‑499215                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                ctgtatgcctctgtctggatggaaggtggctcgccgctgatggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAatata                                                                                                                                                                                                                                                                            >  1:83513‑M2/279‑300 (MQ=255)   gtatgcctctgtctggatggaaggtggctcgccgctgatggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATAtt                                                                                                                                                                                                                                                                          >  1:314402‑M2/277‑300 (MQ=255)                 ggatggaaggtggctcgccgctgatggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagataaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGc                                                                                                                                                                                                                                                            <  2:206101‑M2/38‑1 (MQ=255)                            ggctcgcctctgatggtttacagccaccagcaacagcaggagcgggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctgtcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAAc                                                                                                                                                                                                                                                 <  2:256705‑M2/49‑1 (MQ=255)                                        atggtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAatga                                                                                                                                                                                                                                     >  2:289633‑M2/240‑300 (MQ=255)                                           gtttacagccgccagcaacagcaggcgctggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATgg                                                                                                                                                                                                                                  >  1:299628‑M2/237‑300 (MQ=255)                                                 ctgcgagatgtgtataagagacaggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgtggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATc                                                                                                                                                                                                                            >  1:494319‑M2/231‑300 (MQ=255)                                                        agcaacagcaggcgctggcacaacgtttaccgtagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtattggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGc                                                                                                                                                                                                                     <  2:229539‑M2/77‑1 (MQ=255)                                                                        ggcgcgacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattttggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGc                                                                                                                                                                                                     <  2:8704‑M2/93‑1 (MQ=255)                                                                        ggcacaacgtttaccgtagatgcccgtagcgctgggaatgagctacggctcgccatcacttgaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGc                                                                                                                                                                                                     <  1:165776‑M2/93‑1 (MQ=255)                                                                        ggcacaacgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGc                                                                                                                                                                                                     <  2:297166‑M2/93‑1 (MQ=255)                                                                              acgtttaccggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCg                                                                                                                                                                                               >  1:93293‑M2/202‑300 (MQ=255)                                                                                       ggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactccttgcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAAc                                                                                                                                                                                      <  2:112988‑M2/108‑1 (MQ=255)                                                                                       ggagatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAAc                                                                                                                                                                                      <  1:114976‑M2/108‑1 (MQ=255)                                                                                         agatgcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTg                                                                                                                                                                                    <  1:639123‑M2/110‑1 (MQ=255)                                                                                             gcccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGTGGTAAAGAACTGcgcg                                                                                                                                                                                >  1:280354‑M2/187‑300 (MQ=255)                                                                                               ccgtagcgctgggaatgagctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattcaggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTGATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGcc                                                                                                                                                                              >  2:566903‑M2/185‑300 (MQ=255)                                                                                                                agctacggctcgccatcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCtcgt                                                                                                                                                             >  1:252325‑M2/168‑300 (MQ=255)                                                                                                                               tcactggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGCTGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCggcggtttttctca                                                                                                                                              >  2:91887‑M2/153‑295 (MQ=255)                                                                                                                                  ctggaaagcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattcaggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATcc                                                                                                                                           >  2:412293‑M2/150‑300 (MQ=255)                                                                                                                                        agcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaat                                                                                                                                     <  1:371663‑M2/157‑1 (MQ=255)                                                                                                                                         gcgccgtagatgaactcctggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTAGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGaata                                                                                                                                    >  2:356521‑M2/143‑300 (MQ=255)                                                                                                                                                           tggcagagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAg                                                                                                                  >  1:101010‑M2/125‑300 (MQ=255)                                                                                                                                                                gagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgccttctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGc                                                                                                             >  2:509524‑M2/120‑300 (MQ=255)                                                                                                                                                                gagcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGc                                                                                                             >  1:258060‑M2/120‑300 (MQ=255)                                                                                                                                                                aaacatgtagatcatattgtggggctgccgctttatccgcaattctcctgttctacggtcggtgcggtattggaggaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATTGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGc                                                                                                             <  2:77864‑M2/181‑1 (MQ=255)                                                                                                                                                                  gcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCt                                                                                                           >  1:584810‑M2/118‑300 (MQ=255)                                                                                                                                                                  gcatgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGNTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCt                                                                                                           >  1:485637‑M2/118‑300 (MQ=255)                                                                                                                                                                   catgtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCgtgcatg                                                                                                          >  2:254941‑M2/117‑297 (MQ=255)                                                                                                                                                                      gtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcngtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTNGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGt                                                                                                       <  2:615460‑M2/187‑1 (MQ=255)                                                                                                                                                                      gtagatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGt                                                                                                       <  1:112677‑M2/187‑1 (MQ=255)                                                                                                                                                                        agatcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTag                                                                                                     <  1:132386‑M2/189‑1 (MQ=255)                                                                                                                                                                          atcatattgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTCTCGCTAAAGCTGAGCGGTAAAGAACGGAGCGCCGTCGCGAAGAGCCTCGTCGGNGGTTTTCATCATCCGTGAATAATGCAAAAAGGAGTGCAGCGTGCATGGGTagag                                                                                                   >  2:497740‑M2/110‑300 (MQ=255)                                                                                                                                                                                ttgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGNAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGaa                                                                                             <  1:457815‑M2/197‑1 (MQ=255)                                                                                                                                                                                ttgtggtgctgccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGaa                                                                                             >  1:338352‑M2/104‑300 (MQ=255)                                                                                                                                                                                          gccgctttatccgcaattctcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGgc                                                                                   <  2:607068‑M2/207‑1 (MQ=255)                                                                                                                                                                                                    ccgcaattctccttttctacggtcggtgcggtatgggatgaacttacacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATgcg                                                                         <  2:135143‑M2/217‑1 (MQ=255)                                                                                                                                                                                                    ccgcaattctcctgttctactgtcggtgcggtattggatgaactggcacgcattctggcgcgcaaacgtagcattcctgggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACNGATCGAAATAACGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATgcg                                                                         <  2:161341‑M2/217‑1 (MQ=255)                                                                                                                                                                                                             tcctgttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACg                                                                <  1:272554‑M2/226‑1 (MQ=255)                                                                                                                                                                                                                 gttctacggtcggtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTg                                                            >  1:585461‑M2/71‑300 (MQ=255)                                                                                                                                                                                                                        ggtcggtgcagtatgggatgaactggcacgcattctggcgcgtaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGc                                                     <  2:79261‑M2/237‑1 (MQ=255)                                                                                                                                                                                                                        ggtcggngcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcatnccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCCCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGCCTGCTGATGCGCCGCTAACGTCTGGTAAAGc                                                     <  2:553573‑M2/237‑1 (MQ=255)                                                                                                                                                                                                                             gtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACg                                                >  1:393597‑M2/59‑300 (MQ=255)                                                                                                                                                                                                                             gtgcggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGAGGTAAAGAACTNAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGATAACGTCTGGTAAAGCAGACg                                                >  2:161269‑M2/59‑300 (MQ=255)                                                                                                                                                                                                                                 ggtatgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTg                                            >  1:513344‑M2/55‑300 (MQ=255)                                                                                                                                                                                                                                    atgggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCNAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCa                                         <  2:574516‑M2/249‑1 (MQ=255)                                                                                                                                                                                                                                      gggatgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATc                                       <  1:506222‑M2/251‑1 (MQ=255)                                                                                                                                                                                                                                          tgaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCcagc                                   <  1:52005‑M2/255‑1 (MQ=255)                                                                                                                                                                                                                                           gaactggcacgcattctggcgcgcaaacgtagcattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCATGGGTAGAGTTTGAACTCANAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCcagca                                  >  2:496413‑M2/45‑300 (MQ=255)                                                                                                                                                                                                                                                                           cattccggggatGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAg  <  1:55850‑M2/288‑1 (MQ=255)                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                GTTATGCAGATGAAGGCGATGATTACCCGCAACGTTGCCGCACAACGACTCGTGAACTGGCTTCCGCATTGGGGATGGCACCGGAAAAAGTGATGATGACCTTTCAGTCGCGCTTTGGTCGGGAACCCTGGCTGATGCCTTATACCGACGAAACGCTGAAAATGCTCGGAGAAAAAGGCGTAGGTCATATTCAGGTGATGTGCCCGGGCTTTGCTGCGGATTGTCTGGAGACGCTGGAAGAGATTGCCGAGCAAAACCGTGAGGTCTTCCTCGGTGCCGGCGGGAAAAAATATGAATATATTCCGGCGCTTAATGCCACGCCGGAACATATCGAAATGATGGCTAATCTTGTTGCCGCGTATCGCTAAAGCTGAGCGGTAAAGAACTGAGCGCCGTCGCGAAGAGCCTCGTCGGCGGTTTTCATCATCCGTGAATAATGCAAAAAGGCGTGCAGCGTGCCTGGGTAGAGTTTGAACTCACAGGGCTGCTGATGCGCCGCTAACGTCTGGTAAAGCAGACGGCTGTCATCCAGCAGCGGATCGAACTCCGCCCCGGCAATAAAACAG  >  NC_000913/498650‑499215

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