Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	167,401	Δ5,915 bp		fhuA–[fhuB]	fhuA, fhuC, fhuD, [fhuB]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	167401	173315	5915	56 [0]	[0] 56	fhuA–[fhuB]	fhuA,fhuC,fhuD,[fhuB]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 167400	0 (0.000)	56 (0.860)	52/478	0.1	100%	intergenic (+136/‑84)	mrcB/fhuA	peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB/ferrichrome outer membrane transporter/phage receptor
	?	NC_000913	173316 =	0 (0.000)					coding (1855/1983 nt)	fhuB	iron(III) hydroxamate ABC transporter membrane subunit

TCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT  >  NC_000913/167156‑167646                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       tCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactg                                                                                                                                                                                                                                                  >  1:308450‑M1/1‑245 (MQ=255)            aGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgc                                                                                                                                                                                                                                       >  1:76242‑M1/1‑234 (MQ=255)               cGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctg                                                                                                                                                                                                                                     >  1:418087‑M1/1‑231 (MQ=255)                gcagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgac                                                                                                                                                                                                                                   >  2:440025‑M1/1‑230 (MQ=255)                  agcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactg                                                                                                                                                                                                                                 <  1:477959‑M1/251‑24 (MQ=255)                     aaCAGCCGCAACAGCAACCTGCTCAGCAAGAGAAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtg                                                                                                                                                                                                                              >  2:159856‑M1/1‑225 (MQ=255)                        aGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgg                                                                                                                                                                                                                           >  2:190485‑M1/1‑222 (MQ=255)                              aaCAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggat                                                                                                                                                                                                                     >  2:427175‑M1/1‑216 (MQ=255)                               aCAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatg                                                                                                                                                                                                                    <  1:495122‑M1/251‑37 (MQ=255)                                      ccTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgt                                                                                                                                                                                                             <  1:16154‑M1/251‑44 (MQ=255)                                                       gAAAGACAGCGACGGTGTAGCCGGTTGGATCNAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTGTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccatgccagatcccg                                                                                                                                                                                            <  2:252631‑M1/251‑61 (MQ=255)                                                        aaaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccgg                                                                                                                                                                                           >  2:341883‑M1/1‑190 (MQ=255)                                                          aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCCTCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcg                                                                                                                                                                                         >  1:471920‑M1/1‑188 (MQ=255)                                                               gCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggct                                                                                                                                                                                    <  2:85186‑M1/251‑69 (MQ=255)                                                               gCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggct                                                                                                                                                                                    <  1:14390‑M1/251‑69 (MQ=255)                                                                   cGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctcctg                                                                                                                                                                                >  2:108271‑M1/1‑179 (MQ=255)                                                                            cGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaccttt                                                                                                                                                                       <  2:313520‑M1/251‑82 (MQ=255)                                                                                  gATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggc                                                                                                                                                                 <  1:341883‑M1/251‑88 (MQ=255)                                                                                  gATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggc                                                                                                                                                                 <  2:7767‑M1/251‑88 (MQ=255)                                                                                      aaGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgccgtcaaccttttttggcgcgc                                                                                                                                                             >  2:163274‑M1/1‑160 (MQ=255)                                                                                            atGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaactttatcggcgcgccatatt                                                                                                                                                       >  1:479722‑M1/1‑154 (MQ=255)                                                                                              gTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgtcatatttt                                                                                                                                                     >  1:234841‑M1/1‑152 (MQ=255)                                                                                                            aaCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTACCCCTTCCTTTTCATCTTGTTGTTTATTAACCCTTCAGGAACGCCCAGATTGCGTTCCGCTTGCGAACCttactgctggtgttcgctgcctggtgtgggcggattttgctgtttccgttccagatcccggcggagctactgtcaacctttatcggcgcgtcatatttcatctatttgttga                                                                                                                                        <  2:442921‑M1/250‑113 (MQ=255)                                                                                                              cATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctttttgttgaga                                                                                                                                      >  2:120015‑M1/1‑136 (MQ=255)                                                                                                                  tAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttttcggcgcgccatattttatctatttgttgaggaagc                                                                                                                                  >  2:320493‑M1/1‑132 (MQ=255)                                                                                                                       gTGAAATACCGGATGGCGAGTTGCCCTCCGGTAAAATAACATCCCATCTATGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctttttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagatagcagagcc                                                                                                                            <  2:479722‑M1/251‑125 (MQ=255)                                                                                                                          aaaTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgtt                                                                                                                         <  1:159856‑M1/251‑128 (MQ=255)                                                                                                                                cGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcaaagccgttaatttt                                                                                                                   >  1:99291‑M1/1‑118 (MQ=255)                                                                                                                                                  ccGGTAAAATAACATCCCATCTAAGATATTACCCCTTTCTTTTCATCGGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTTCCGCTGGCGAACCtttcggctggggttcgctgcctggtgtgggcggtgggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaagttttccatataacggccctg                                                                                                  <  2:169024‑M1/250‑151 (MQ=255)                                                                                                                                                    gttaaAATAACATCTCATCTAAGATTTTAACCCTTTCTTTTCATCTGGTTGGTTATTAACCCTTGAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgcggactggtgtgggcggatggtgctgtttccattccagatccaggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtac                                                                                               <  2:480441‑M1/249‑154 (MQ=255)                                                                                                                                                      tAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctattttttgagaaagcagacccgttaatttttccatataactgccctgtccgc                                                                                             >  2:138242‑M1/1‑96 (MQ=255)                                                                                                                                                                    aTCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGATCGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaactttacaggcgcgccatattttatctatttgttgagaaagaagagtcgttaatttttccatataacggcactgaccgccaggacggataa                                                                                 >  2:441163‑M1/1‑82 (MQ=255)                                                                                                                                                                       tacgataTTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTAAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttcccttccagatcccggcggggctgctgtagacacttattggcgcgccatattttatctttttgttgcgaaaacaaagccgttaatttttccctataaaggccctgtacaccatgacgagcaaagcgt                                                                            >  2:477934‑M1/4‑79 (MQ=255)                                                                                                                                                                          gATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttcatttttccatataacggcccttttctcctggacgggtaaagcgtttt                                                                         >  2:56040‑M1/1‑76 (MQ=255)                                                                                                                                                                               tAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttggtgagaaagcagagccgttaatttttccatataacggccctgtaagcctggacggataaaacgtttacgagc                                                                    >  1:287948‑M1/1‑71 (MQ=255)                                                                                                                                                                                aaCCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgac                                                                   <  1:258207‑M1/251‑182 (MQ=255)                                                                                                                                                                                aaCCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcga                                                                    <  1:120015‑M1/250‑181 (MQ=255)                                                                                                                                                                                   ccTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatccgcgcgccatattttatctttttgtttagaaaacagaaccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgagaaca                                                                 >  2:105554‑M1/1‑67 (MQ=255)                                                                                                                                                                                    ctttctttTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGAGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaaacagagccgttaatttttccatataacggccctgtacgcccggacggataaagcgtttacgcgacaacg                                                               >  1:151207‑M1/1‑66 (MQ=255)                                                                                                                                                                                          tttCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatg                                                           <  2:117622‑M1/249‑190 (MQ=255)                                                                                                                                                                                           ttCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGCACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccattttttttctatttgttgagaaagcaggaccgttaatttttccatataaaggccctgtacccctggacgggtaaagcgttttcgcgacaacggattccc                                                        >  1:218108‑M1/1‑59 (MQ=255)                                                                                                                                                                                           ttCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggcggctgtcaacctttatcggcgcgccatattttttctattttttgtgaaaacagagccgttaatttttccatataacggccctgtacgcctgggcggataaatggtttacgcgaccaacgatgccc                                                        >  2:90547‑M1/1‑59 (MQ=255)                                                                                                                                                                                             cATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcatagccgttaatttttccatataacggcccggtacgcctggacggataaagcgtttactcgacaacggaagcccg                                                       >  1:345878‑M1/1‑57 (MQ=255)                                                                                                                                                                                                  ggTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgat                                                     <  1:440025‑M1/247‑196 (MQ=255)                                                                                                                                                                                                   gttgttTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgac                                                <  1:56040‑M1/251‑201 (MQ=255)                                                                                                                                                                                                       tttATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccctattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtaagcctggacggataaagcgtttacgcgacaacggatgcccgctgcgacgctg                                            >  1:71471‑M1/1‑47 (MQ=255)                                                                                                                                                                                                        ttattaACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctg                                            <  1:11531‑M1/250‑205 (MQ=255)                                                                                                                                                                                                           ttaACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaaattttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggctgcccgatgcgacgcttgcgc                                        >  1:215783‑M1/1‑43 (MQ=255)                                                                                                                                                                                                            taACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgc                                        <  1:284724‑M1/250‑209 (MQ=255)                                                                                                                                                                                                                   tCAGGAACGCTCAGATTGCGTACCGCTTGCTAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggcgctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagacgttaatttttccatataacggccctgtacgcctgtcaggataaagcgattacgcgaaaacggattcccgcagcgccgctgtcgcgtcttatc                                >  2:498278‑M1/1‑35 (MQ=255)                                                                                                                                                                                                                    cAGGAACGCTCAGATTGCGTACCGCTTGCTAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcgtcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcggttacgcgaccaaggatgcccgctgcgccgctggcgcgtcttatcg                               >  2:330523‑M1/1‑34 (MQ=255)                                                                                                                                                                                                                      ggAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccggtaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatccgg                             >  2:379205‑M1/1‑32 (MQ=255)                                                                                                                                                                                                                        aaCGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggaaggataaagcgtttacgcgacaacggatgcccgatgcgccgctggcgcgtcttatcgggcc                           >  1:86909‑M1/1‑30 (MQ=255)                                                                                                                                                                                                                          cGCTCCGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggcttctgtcaacctttatcggcgcgccatattttatctttttgttgagaaatcagagccgttaatttttccatataacggccctgtacccctgggccgataaaaagtttacgccacaacggatgcccgatgcgacgctggcgcgtcttttcgggact                          >  1:93178‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                                 aTTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccctataacggccctgtacgcctggacggaaaaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcggcttatcgggcctacaaacgg                  >  1:106064‑M1/1‑21 (MQ=255)                                                                                                                                                                                                                                                 gAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcagcctttctcggcgcgccatattttatctatttgttgagaaagcagagccgttaattttttcatataacggccctgtacgcctgggcggataaagcgtttacgcgacaaccgatgcccgatgcgacgctggcccgtcttatcgggcctacaaacggcccccccacgtaaggc  >  1:3573‑M1/1‑5 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       TCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCAT  >  NC_000913/167156‑167646

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