Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
MC JC	NC_000913	4,001,645	Δ5 bp	100%	coding (220‑224/951 nt)	corA →	magnesium/nickel/cobalt transporter

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4001645	4001649	5	75 [0]	[0] 75	corA	magnesium/nickel/cobalt transporter

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4001644	0 (0.000)	75 (1.500)	60/488	NT	100%	coding (219/951 nt)	corA	magnesium/nickel/cobalt transporter
	?	NC_000913	4001650 =	0 (0.000)					coding (225/951 nt)	corA	magnesium/nickel/cobalt transporter

GTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCG  >  NC_000913/4001380‑4001930                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                         tctgtgtgcccaccgaactgtccgatattttaagcattgggagtcccggtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGaa                                                                                                                                                                                                                                                                        >  2:76142‑M2/271‑289 (MQ=255)  ctgtgtgcccaccgaactgtccgatattttaagcattgggagtcccggtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAg                                                                                                                                                                                                                                                                       <  2:422312‑M2/20‑1 (MQ=255)              cgaactgtccgatattttaagcattgggagtcccggtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATc                                                                                                                                                                                                                                                           <  2:287455‑M2/32‑1 (MQ=255)                 actgtccgatattttaagcattgggagtcccggtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                                                                                                                                                                        <  2:390426‑M2/35‑1 (MQ=255)                 actgtccgatattttaagcattgggagtcccggtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                                                                                                                                                                        <  2:6988‑M2/35‑1 (MQ=255)                             tttaagcattgggagtcccggtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCa                                                                                                                                                                                                                                            <  1:321256‑M2/47‑1 (MQ=255)                               taagcattgggagtcccggtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTcc                                                                                                                                                                                                                                             >  2:481075‑M2/241‑286 (MQ=255)                               taagcattgggagtcccggtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCCTTGAAGATGCGGAAGATCACGCCGGTAACTcc                                                                                                                                                                                                                                             >  2:280112‑M2/241‑286 (MQ=255)                                      tgggagtcccggtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCAt                                                                                                                                                                                                                                   >  1:136195‑M2/234‑289 (MQ=255)                                                 gtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                                                                                                                                                                        <  2:38731‑M2/35‑1 (MQ=255)                                                 gtcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACg                                                                                                                                                                                                                                                        >  1:38731‑M2/223‑257 (MQ=255)                                                  tcatgctgagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATTACGCCGGTAACTCCACTGTGGCATTTACCATCCGTg                                                                                                                                                                                                                       >  2:142172‑M2/222‑289 (MQ=255)                                                      gctgagcgcatttgaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGt                                                                                                                                                                                                                                                   >  1:479315‑M2/218‑257 (MQ=255)                                                      gctgagcgcatttaaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGt                                                                                                                                                                                                                                                   <  2:479315‑M2/40‑1 (MQ=255)                                                         gagcgcatttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcg                                                                                                                                                                                                                <  1:207840‑M2/75‑1 (MQ=255)                                                               atttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGtt                                                                                                                                                                                                          <  2:136195‑M2/81‑1 (MQ=255)                                                                 ttcaactggaaaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGttt                                                                                                                                                                                                         >  1:391190‑M2/207‑288 (MQ=255)                                                                           aaataaccgactgacccggctggaagtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTgcgtga                                                                                                                                                                                              >  1:9675‑M2/197‑289 (MQ=255)                                                                                                   agtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCg                                                                                                                                                                         >  1:259067‑M2/173‑286 (MQ=255)                                                                                                   agtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTc                                                                                                                                                                       <  2:391190‑M2/116‑1 (MQ=255)                                                                                                   agtcgaagagtcacaaccccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCt                                                                                                                                                                      <  1:76142‑M2/117‑1 (MQ=255)                                                                                                   agtcgaagaggcgcaaccccttgtaaatgcagtatggtttgttgttgtgggaccgggcgaggacgggcgactgcgcttacaaccggaacttggccagggcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGt                                                                                                                                                                        <  2:448429‑M2/115‑1 (MQ=255)                                                                                                            ggcacaaacccttgtaaatgcagtagggattgatcgtgtcgaaccggacgacgacgagcggctgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGt                                                                                                                                                               <  2:287692‑M2/124‑1 (MQ=255)                                                                                                               acgacaccttgttaaaacagtatggatagatcgtgtcgaaccggacgacgacgagcaactgcgcgtacaatctgaacttggccagggcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttcattgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCg                                                                                                                                                          <  2:308711‑M2/129‑1 (MQ=255)                                                                                                                  accccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTGTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcg                                                                                                                                                                                                                <  2:265255‑M2/75‑1 (MQ=255)                                                                                                                  accccttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcg                                                                                                                                                                                                                >  1:265255‑M2/158‑232 (MQ=255)                                                                                                                      cttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGt                                                                                                                                                   >  1:364755‑M2/154‑289 (MQ=255)                                                                                                                       ttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTa                                                                                                                                                  >  1:308711‑M2/153‑289 (MQ=255)                                                                                                                       ttgtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTa                                                                                                                                                  >  2:37161‑M2/153‑289 (MQ=255)                                                                                                                        tgtaaatgcagtatggattgatcttggcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAg                                                                                                                                                 <  1:8020‑M2/138‑1 (MQ=255)                                                                                                                         gtaaatgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATTACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGc                                                                                                                                                <  1:142172‑M2/139‑1 (MQ=255)                                                                                                                             atgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGt                                                                                                                                            <  1:218434‑M2/143‑1 (MQ=255)                                                                                                                              tgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTc                                                                                                                                           >  1:351033‑M2/146‑289 (MQ=255)                                                                                                                              tgcagtatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTc                                                                                                                                           >  1:296915‑M2/146‑289 (MQ=255)                                                                                                                                 agtatggattgatattgtcgaaccggacgacgacgagcgacggcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGGGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGa                                                                                                                                         <  1:233797‑M2/146‑1 (MQ=255)                                                                                                                                   tatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcg                                                                                                                                                                                                                <  1:162418‑M2/75‑1 (MQ=255)                                                                                                                                   tatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCCTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACCGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCCTCTGTATCGTATTCTTGCCCGTAGCCAGTCGATgg                                                                                                                                      >  2:233797‑M2/141‑289 (MQ=255)                                                                                                                                   tatggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCCCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGgtcg                                                                                                                                                                                                                >  2:162418‑M2/141‑215 (MQ=255)                                                                                                                                    atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtt                                                                                                                                                                             <  1:144196‑M2/110‑1 (MQ=255)                                                                                                                                    atggattgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCtt                                                                                                                                                                             >  2:144196‑M2/140‑249 (MQ=255)                                                                                                                                        attgatcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGAc                                                                                                                                 >  1:206362‑M2/136‑289 (MQ=255)                                                                                                                                           gatcttgtcgaacaggacgacgacggtcgactgcgcgtacaatcggaacttggccagggcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGt                                                                                                                              <  2:259067‑M2/157‑1 (MQ=255)                                                                                                                                            atcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTa                                                                                                                             <  2:296915‑M2/158‑1 (MQ=255)                                                                                                                                             tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTa                                                                                                                             <  1:223089‑M2/158‑1 (MQ=255)                                                                                                                                             tcttgtcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTa                                                                                                                             >  2:223089‑M2/131‑288 (MQ=255)                                                                                                                                                  tcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTa                                                                                                                                                                                                        <  2:39374‑M2/83‑1 (MQ=255)                                                                                                                                                  tcgaaccggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTa                                                                                                                                                                                                        >  1:39374‑M2/126‑208 (MQ=255)                                                                                                                                                        cggacgacgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGt                                                                                                                 <  1:91020‑M2/170‑1 (MQ=255)                                                                                                                                                          gacgacgacgggcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGtt                                                                                                                <  2:289005‑M2/171‑1 (MQ=255)                                                                                                                                                               cgacgagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTtgctgc                                                                                                           <  2:351033‑M2/176‑1 (MQ=255)                                                                                                                                                                   gagcgactgcgcgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGAt                                                                                                      <  1:481075‑M2/181‑1 (MQ=255)                                                                                                                                                                          tgcgcgtacaatctgaacttggccagagcctggcaacccgcgcggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCg                                                                                               >  2:116115‑M2/102‑289 (MQ=255)                                                                                                                                                                              cgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGtt                                                                                                                >  1:173674‑M2/98‑268 (MQ=255)                                                                                                                                                                              cgtacaatctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGtt                                                                                                                <  2:173674‑M2/171‑1 (MQ=255)                                                                                                                                                                                 acaatctgaacttggccagagcctggcaagccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCa                                                                                         <  1:110463‑M2/194‑1 (MQ=255)                                                                                                                                                                                      ctgaacttggccagagcctggcaacccgcgcggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCg                                                                                   <  1:116115‑M2/200‑1 (MQ=255)                                                                                                                                                                                      ctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                                                                                                                                                                  >  2:87664‑M2/90‑194 (MQ=255)                                                                                                                                                                                      ctgaacttggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGcc                                                                                                                                                                                  <  1:87664‑M2/105‑1 (MQ=255)                                                                                                                                                                                              ggccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTg                                                                           >  1:184335‑M2/82‑289 (MQ=255)                                                                                                                                                                                               gccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGt                                                                                                                                                                                                                                                   <  1:491814‑M2/40‑1 (MQ=255)                                                                                                                                                                                               gccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGt                                                                                                                                                                                                                                                   >  2:491814‑M2/81‑120 (MQ=255)                                                                                                                                                                                               gccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTgg                                                                          >  2:286365‑M2/81‑289 (MQ=255)                                                                                                                                                                                               gccagagcctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTg                                                                           >  1:96520‑M2/81‑288 (MQ=255)                                                                                                                                                                                                      cctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATc                                                                                                     >  1:463657‑M2/74‑255 (MQ=255)                                                                                                                                                                                                      cctggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATc                                                                                                     <  2:463657‑M2/182‑1 (MQ=255)                                                                                                                                                                                                        tggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGaaa                                                                 <  1:37161‑M2/218‑1 (MQ=255)                                                                                                                                                                                                         ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGCAGGTACCGCCTACGAGTTGCTGCTGGATCTGTTCGAACCCAAACTCGAACAGTTTGATAATGAAAt                                                                >  2:381699‑M2/71‑289 (MQ=255)                                                                                                                                                                                                         ggcaacccgcccggaactggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAAt                                                                <  1:280112‑M2/219‑1 (MQ=255)                                                                                                                                                                                                                        actggaagacatcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCtata                                                  <  2:206362‑M2/233‑1 (MQ=255)                                                                                                                                                                                                                                  atcgaagcatcggcacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGa                                        >  1:75253‑M2/46‑288 (MQ=255)                                                                                                                                                                                                                                               cacgtttctttgaagaggacgacggcctgcatCGCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGGGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCgg                           <  2:211983‑M2/256‑1 (MQ=255)                                                                                                                                                                                                                                               cacgtttctttgaagacgacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCgg                           <  1:286365‑M2/256‑1 (MQ=255)                                                                                                                                                                                                                                                                gacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCtata                                                  >  1:271632‑M2/16‑248 (MQ=255)                                                                                                                                                                                                                                                                gacgacggcctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCtata                                                  <  2:271632‑M2/233‑1 (MQ=255)                                                                                                                                                                                                                                                                        cctgcatCTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCg  <  2:184335‑M2/281‑1 (MQ=255)                                                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                                                         GTGCCCACCGAACTGTCCGATATTTTAAGCATTGGGAGTCCCGGTCATGCTGAGCGCATTTCAACTGGAAAATAACCGACTGACCCGGCTGGAAGTCGAAGAGTCACAACCCCTTGTAAATGCAGTATGGATTGATCTTGTCGAACCGGACGACGACGAGCGACTGCGCGTACAATCTGAACTTGGCCAGAGCCTGGCAACCCGCCCGGAACTGGAAGACATCGAAGCATCGGCACGTTTCTTTGAAGACGACGACGGCCTGCATATTCACTCCTTCTTCTTCTTTGAAGATGCGGAAGATCACGCCGGTAACTCCACTGTGGCATTTACCATCCGTGATGGTCGTCTGTTTACTCTGCGTGAGCGTGAACTGCCCGCTTTTCGTCTGTATCGTATGCGTGCCCGTAGCCAGTCGATGGTAGACGGTAACGCCTACGAGTTGCTGCTGGATCTGTTCGAAACCAAAATCGAACAGTTGGCAGATGAAATTGAAAATATCTATAGCGACCTGGAGCAGTTGAGCCGGGTGATTATGGAAGGGCATCAGGGCG  >  NC_000913/4001380‑4001930

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