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	77 [0]	[0] 77	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)	76 (0.850)	69/492	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

GTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA  >  NC_000913/167150‑167645                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      gTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAcc                                                                                                                                                                                                                                                       <  1:210820/251‑1 (MQ=255)   ttGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCtt                                                                                                                                                                                                                                                     >  1:318722‑M1/1‑249 (MQ=255)                   gCCGCGGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCGTGCGAACCttactgctggtgttcgct                                                                                                                                                                                                                                     <  1:560622‑M1/251‑19 (MQ=255)                   gCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgct                                                                                                                                                                                                                                     <  1:74268‑M1/251‑19 (MQ=255)                     cGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctg                                                                                                                                                                                                                                    >  1:31145‑M1/1‑231 (MQ=255)                      gcagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgac                                                                                                                                                                                                                                  >  1:429612‑M1/1‑230 (MQ=255)                      gcagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgac                                                                                                                                                                                                                                  <  1:118987‑M1/251‑22 (MQ=255)                                ccGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggc                                                                                                                                                                                                                        <  2:508392‑M1/251‑32 (MQ=255)                                       aGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCTTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTACCCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttacggctggtgttcgctgactggtgtgggcggatggc                                                                                                                                                                                                                 <  2:142728‑M1/251‑39 (MQ=255)                                               gCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttc                                                                                                                                                                                                         <  2:481978‑M1/251‑47 (MQ=255)                                                   aGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccatt                                                                                                                                                                                                     <  1:547397‑M1/251‑51 (MQ=255)                                                    gCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattc                                                                                                                                                                                                    <  2:539835‑M1/251‑52 (MQ=255)                                                           cAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcc                                                                                                                                                                                             >  1:525806‑M1/1‑193 (MQ=255)                                                            aGAAAGACAGCGACGGGGTAGCCGGTTGGATCAAGGAGATGTTTGGTAGTAATTAACATCGAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatccc                                                                                                                                                                                            <  2:140200‑M1/251‑60 (MQ=255)                                                                    aGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggc                                                                                                                                                                                    <  1:317933‑M1/251‑68 (MQ=255)                                                                    aGCGACGGTGAAGACGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGCTGGCGAGTTGCCATCCGGTAAACTAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGCACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatgtagctgtttccattccagatcccggcggggg                                                                                                                                                                                    >  2:101941‑M1/1‑184 (MQ=255)                                                                     gCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggct                                                                                                                                                                                   <  1:208707‑M1/251‑69 (MQ=255)                                                                     gCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggct                                                                                                                                                                                   <  1:310493‑M1/251‑69 (MQ=255)                                                                       gACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTACCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggcctc                                                                                                                                                                                 >  2:259307‑M1/1‑181 (MQ=255)                                                                        aCGGTGTAGCCGGTTGGATAAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGACATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATCTTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgct                                                                                                                                                                                >  1:96652‑M1/1‑180 (MQ=255)                                                                         cGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctg                                                                                                                                                                               <  2:429612‑M1/251‑73 (MQ=255)                                                                          ggTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgg                                                                                                                                                                              >  2:153397‑M1/1‑178 (MQ=255)                                                                                   ggTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgttaacctttt                                                                                                                                                                     >  1:141672‑M1/1‑169 (MQ=255)                                                                                   ggTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgggggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttt                                                                                                                                                                     >  1:570273‑M1/1‑169 (MQ=255)                                                                                    gTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGCTATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttct                                                                                                                                                                    >  1:655181‑M1/1‑168 (MQ=255)                                                                                    gTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtaaaactttat                                                                                                                                                                    >  1:107017‑M1/1‑168 (MQ=255)                                                                                           cAAGGATATGTTTGGCAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCAATCCGGTAAAAAAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTAAGGAACGCTCAGATTGCGTACCCCTTGCGAACCttactgctggtgttcgctgcctggtgggggcgattggtgctgtttcctttcaagatccggccggggctgttgtccacctttatcggcggc                                                                                                                                                             >  2:68763‑M1/1‑161 (MQ=255)                                                                                             aGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttattggctcgcc                                                                                                                                                           >  1:692592‑M1/1‑159 (MQ=255)                                                                                                           agtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTATTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattacagatcccggcggggctggtggccacctttatctgcccgccctattttatctttt                                                                                                                                             >  1:291044‑M1/1‑145 (MQ=255)                                                                                                             taATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttc                                                                                                                                                                                                         <  1:453575‑M1/189‑47 (MQ=255)                                                                                                             taATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttc                                                                                                                                                                                                         >  2:453575‑M1/1‑143 (MQ=255)                                                                                                                ttAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccctattttatctatttgttt                                                                                                                                        >  1:197898‑M1/1‑140 (MQ=255)                                                                                                                         aaGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcag                                                                                                                               <  1:615773‑M1/251‑121 (MQ=255)                                                                                                                           gCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagag                                                                                                                             <  2:379311‑M1/251‑123 (MQ=255)                                                                                                                            cGTGAAATACCGGATGGCGTGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGACCGCTCAGATTGCGTTCCGCTTGCGAACCttactgctggtgttctctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagc                                                                                                                            <  1:481218‑M1/251‑124 (MQ=255)                                                                                                                               gAAATACCGGATGGCCAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgcggtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgt                                                                                                                         <  2:394940‑M1/251‑127 (MQ=255)                                                                                                                                      cGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggcgctgctgtcaacacttatcggcgcgccatattttatctaatttttgggaaagcagagccgttaatttt                                                                                                                  >  2:625239‑M1/1‑118 (MQ=255)                                                                                                                                       ggATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGGTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagcgccgttaattttt                                                                                                                 <  1:495412‑M1/251‑135 (MQ=255)                                                                                                                                           ggCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCATTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggagggtgctgtttccactccagatcccggcggggctgctgtgcacctttatagtcgccccgtattttatctctttgttgagaaagccaagccgtaaatctttccat                                                                                                             >  2:724191‑M1/1‑113 (MQ=255)                                                                                                                                             cGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCAGCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgataaagcagggccgttaatttttccatat                                                                                                           <  2:635745‑M1/251‑141 (MQ=255)                                                                                                                                              gAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatata                                                                                                          <  1:466694‑M1/251‑142 (MQ=255)                                                                                                                                               agtggccatCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataa                                                                                                         <  2:306040‑M1/247‑143 (MQ=255)                                                                                                                                                  tGCCATCCGGTAAAATAATATCCCATCTAAGATATTAACCCTTCCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGATAAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggcgctgtttccattccagttcccggcggggctgctgtcaacctttctcggcgccccacattttatctatttttttagaaagccgagccggtaatttttccattaaaatg                                                                                                      >  2:615252‑M1/1‑106 (MQ=255)                                                                                                                                                  tGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacg                                                                                                       <  1:567440‑M1/250‑145 (MQ=255)                                                                                                                                                   gCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcatcctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggc                                                                                                     <  1:568721‑M1/251‑147 (MQ=255)                                                                                                                                                    ccatcccgttaacgatcatcccatctaagaaatTAACCCTTTCTTTTCAACTGGTTTTTTATTAACCCTTCAGGATCCCTCAGATTGCGTCCCGCGTGCGAACCttactgctggtgttcgcggactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaccctgtatcggcgcgccagattttatctgtttgttgagaaagcagcgccgttaatttttccatatgacggcc                                                                                                    <  2:291044‑M1/235‑148 (MQ=255)                                                                                                                                                       tCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgcccttttttttctattttttgggaaagcagaaccgttaatttttccttttaaaggccctt                                                                                                 >  2:180699‑M1/1‑101 (MQ=255)                                                                                                                                                          ggTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtac                                                                                              <  1:152275‑M1/251‑154 (MQ=255)                                                                                                                                                                  aaCATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgtgaaagcagagccgttaatttttccatattacggccctgcaagcctggac                                                                                      >  2:134599‑M1/1‑90 (MQ=255)                                                                                                                                                                    cATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggagcgccatattttatctatttggtgggaaagcagagccgttaattttttcctataacggccctttacgcccggcagg                                                                                    >  2:434045‑M1/1‑88 (MQ=255)                                                                                                                                                                     aTCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaccgttatcggcgcgccatattttatctatttgttgagaaagcagagccggtaatttttccatataaaggccctgtacgcctggacggg                                                                                   >  1:472058‑M1/1‑87 (MQ=255)                                                                                                                                                                            cTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagaaccggtaatttttccatataaaggccctgtacgcctggacggataaag                                                                              >  1:222912‑M1/1‑80 (MQ=255)                                                                                                                                                                              aaGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgt                                                                           <  1:386290‑M1/250‑173 (MQ=255)                                                                                                                                                                                cccctcttccgatctttttcatctggttgtttatTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgttta                                                                        <  2:138291‑M1/236‑176 (MQ=255)                                                                                                                                                                                 atatTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttac                                                                       <  1:153397‑M1/251‑177 (MQ=255)                                                                                                                                                                                   atTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgtttagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgc                                                                     >  1:629764‑M1/1‑73 (MQ=255)                                                                                                                                                                                        cccTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGGTTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgttttcgcgaca                                                                 >  1:128214‑M1/1‑68 (MQ=255)                                                                                                                                                                                         ccTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggccgctgtcaccctttatcggcgcggcatattttttctatttgttgagaatacagagccggttaattttccatataaaggcccttctagcctggaccgataaaggtgttctgcgacaac                                                               >  1:393242‑M1/1‑67 (MQ=255)                                                                                                                                                                                            ttctttTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttttctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggccggataaagcgtttacgcgacaacggc                                                            >  2:77758‑M1/1‑64 (MQ=255)                                                                                                                                                                                               ttttCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgggaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgttta                                                                        >  1:138291‑M1/1‑61 (MQ=255)                                                                                                                                                                                                tttCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgcaaacggatgcc                                                        >  1:391712‑M1/1‑60 (MQ=255)                                                                                                                                                                                                   cATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccg                                                      >  1:615313‑M1/1‑57 (MQ=255)                                                                                                                                                                                                   cATCGGGTTGTTTATTAACCCTCAAGGAACGCTCAGATGGCGTACCGCTTGCGACCCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacgtatgcccg                                                      <  2:699857‑M1/250‑194 (MQ=255)                                                                                                                                                                                                        ggTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttaaccgacaacggatgaccgaagccc                                                >  2:227350‑M1/1‑52 (MQ=255)                                                                                                                                                                                                                attaACCCGTCAGGAACGCTCAGATTGCGTACCGATTGCGACCCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcgggggtgctgtcaacctttatcggcgcgccatattttttctatttgttgtgaaagccgagccgttaatttttccctttaacggccctgtagccctggacggataaagcgtttacgcgacaacggatgcccgagtcggagctggcg                                        >  2:210161‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                   atcccttCATGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgcggactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcg                                      <  2:75260‑M1/248‑210 (MQ=255)                                                                                                                                                                                                                   aaCCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgg                                     >  1:146693‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                    aCCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgggaaagcagagccgttaatttttccatataacagccctgtacgcctggacggctgacacggttaagcgtcaacggaagccccatgcgacgcaggagcgtc                                    >  1:372344‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                       cTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcaaagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttt                                 >  1:588581‑M1/1‑37 (MQ=255)                                                                                                                                                                                                                         tCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtctta                                 <  1:487035‑M1/249‑215 (MQ=255)                                                                                                                                                                                                                         gcaGGAACGCTCAGATTGCGTCCCGCTTGCGAACCttacgtctggtgttcgctgactggtgtgggcggctggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccagataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttag                                <  2:43599‑M1/249‑216 (MQ=255)                                                                                                                                                                                                                          cAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtaagcctggacggataaagcggttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcg                              >  1:559867‑M1/1‑34 (MQ=255)                                                                                                                                                                                                                           aGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctggcaacctttatcggcgcgccctattttatcttttttttgaggaagccgagccgttaatttttccatataacggccctggacgcctgggcgggaaaagcggttacgcgacaacggntgcccgatgcggagccggcgggtcttaaccg                             >  2:360462‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                                cGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccctattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgttttcgcgacaacgggtgccngacgcgccgccggcgcctgttattgggcccc                        >  2:346651‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                                    caatttgcgtACCGCTTGCGAACCttactgctggttttcgctgacgggtgtgggcgtatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaa                    <  2:671713‑M1/247‑228 (MQ=255)                                                                                                                                                                                                                                                  cTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaaacagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgccacgctggcgcgtcttatcgggcctacaaacggccccgaatcgt      >  1:241307‑M1/1‑10 (MQ=255)                                                                                                                                                                                                                                                      cGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccccgaatcgtaggg  >  1:493284‑M1/1‑6 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                      GTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCA  >  NC_000913/167150‑167645

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