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	60 [0]	[0] 60	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)	60 (0.920)	52/490	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

TTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATG  >  NC_000913/167152‑167647                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        ttGATCAGTCTTCCCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCt                                                                                                                                                                                                                                                        <  2:394478‑M1/250‑2 (MQ=255)           ttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtg                                                                                                                                                                                                                                             >  1:185435‑M1/1‑239 (MQ=255)            tctcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgt                                                                                                                                                                                                                                            >  2:492557‑M1/1‑238 (MQ=255)                 gCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgct                                                                                                                                                                                                                                       <  1:521070‑M1/251‑19 (MQ=255)                       gcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactgg                                                                                                                                                                                                                                 >  1:131550‑M1/1‑227 (MQ=255)                       gcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactgg                                                                                                                                                                                                                                 >  1:89022‑M1/1‑227 (MQ=255)                          ccaGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgt                                                                                                                                                                                                                              <  2:89022‑M1/250‑28 (MQ=255)                              ccGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgcgggc                                                                                                                                                                                                                          <  1:504225‑M1/251‑32 (MQ=255)                                             gCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtt                                                                                                                                                                                                             <  1:233923‑M1/249‑45 (MQ=255)                                                  gCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgttt                                                                                                                                                                                                            <  1:213115‑M1/245‑46 (MQ=255)                                                  gCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgttt                                                                                                                                                                                                            >  2:213115‑M1/1‑200 (MQ=255)                                                     agagCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccag                                                                                                                                                                                                   <  1:151490‑M1/251‑55 (MQ=255)                                                         cAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATGTGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAAGCCTTCAGGAACGCTCAGATTGCGGACCGCTTGCGAACCttactgatggtgttcgctgactggtgtgggagggtggtgctgtttccattccagatcc                                                                                                                                                                                               >  1:195579‑M1/1‑193 (MQ=255)                                                              aGACAGCGACGGTGTAGCCGGTTGGCTCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATACGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTCCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctgttgttcgctgactggtgtgggcggatggtgatgtttccattccagctcccggcg                                                                                                                                                                                          >  1:459916‑M1/1‑188 (MQ=255)                                                                  aGCGACGGTGTAGCCGGTTGGATCAAGGCTATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgagggcggatggtgctgtttcccttccagatcccggcggggg                                                                                                                                                                                      >  2:298139‑M1/1‑184 (MQ=255)                                                                  aGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggc                                                                                                                                                                                      >  1:533785‑M1/1‑184 (MQ=255)                                                                         gtgtAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtc                                                                                                                                                                               <  1:155726‑M1/251‑75 (MQ=255)                                                                                cGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaccttt                                                                                                                                                                        <  1:17122‑M1/251‑82 (MQ=255)                                                                                     ggATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcg                                                                                                                                                                    <  1:492557‑M1/250‑86 (MQ=255)                                                                                        tCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgc                                                                                                                                                                >  2:300727‑M1/1‑162 (MQ=255)                                                                                           aGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGGTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgcc                                                                                                                                                             <  1:66104‑M1/251‑93 (MQ=255)                                                                                                atGTTTGGTAGTAATTAACATCTAAGAGTGAAATACCGGATGGCGAGTTGCCATCCGGAAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggcgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgcttgccaccgttatcggcgcgtcatatt                                                                                                                                                        >  2:344768‑M1/1‑154 (MQ=255)                                                                                                  gttgggtagtAATTAACATCCAATCGTGAANTGCCGGATGGCGGGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTTTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCTCTTGCGAACCttactgctggtgtgcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatctgcgcgccatattt                                                                                                                                                       <  2:239141‑M1/246‑99 (MQ=255)                                                                                                    ttGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCGTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatatttta                                                                                                                                                     <  1:49002‑M1/250‑101 (MQ=255)                                                                                                             aTTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtc                                                                                                                                                                               <  1:295352‑M1/215‑75 (MQ=255)                                                                                                             aTTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtc                                                                                                                                                                               >  2:295352‑M1/1‑141 (MQ=255)                                                                                                                                    cGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgtttggaaagaagagccggtaatttt                                                                                                                    >  2:167659‑M1/1‑118 (MQ=255)                                                                                                                                     ggATGGCGAGTTGCCATCCGGTAAAATAACCTCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaattttt                                                                                                                   <  1:192346‑M1/251‑135 (MQ=255)                                                                                                                                       aTGGCGAGTTGCCATCCGTTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttcc                                                                                                                 <  2:47134‑M1/251‑137 (MQ=255)                                                                                                                                       aTGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttcc                                                                                                                 <  1:140017‑M1/251‑137 (MQ=255)                                                                                                                                                    aTCCGGTAAAATAACATCCCATCAAATATATNAACCCTTTCTTTTCATCTGGTTGTTTATTATCCCTTCAGGAACGCTCAGATTGCGTACCGCTTTCGAACCttactgcttgtgttcgctgactggtgttggcggatggttctgtttccattccagatcccggcgggtctgctgtcaacctttatcggcgcgccatattttatctgtttgttgataaagcagagccgttaattgttccatataacggccct                                                                                                    <  2:250500‑M1/251‑150 (MQ=255)                                                                                                                                                         gTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGCTGTTTATAAACCCTTAAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacatttctccgcgcgccctattttatctcttttatgaggaaacagagccgttaatttttccatataacgggcccgtacg                                                                                               >  2:482301‑M1/1‑97 (MQ=255)                                                                                                                                                           aaaaTAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctattttttgagaaagcagagccgttaatttttccatataacggccctgtacgcc                                                                                             >  1:338234‑M1/1‑95 (MQ=255)                                                                                                                                                            aaaaaacatcCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcct                                                                                            <  2:338234‑M1/247‑158 (MQ=255)                                                                                                                                                            aaaTACCATCCCATCTAAGATATTAACCCTTTCTTTTAATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGCACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattc                                                                                                                                                                                                      >  2:421616‑M1/1‑94 (MQ=255)                                                                                                                                                            aaaTAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattc                                                                                                                                                                                                      <  1:421616‑M1/145‑52 (MQ=255)                                                                                                                                                              aTAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgggaaagcagagccgttaatttttccatataagatcgggagagcacacc                                                                                          >  1:350860‑M1/1‑92 (MQ=255)                                                                                                                                                              aTAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatata                                                                                                            <  2:350860‑M1/233‑142 (MQ=255)                                                                                                                                                                    tCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcgggggtgctgtcaaccttttacggcgcgccatattttaactttttgttgagaaagcagaggcgttaatttttccatataacggaccggtacgcctggacgggt                                                                                    >  1:205411‑M1/1‑86 (MQ=255)                                                                                                                                                                            aaaatatgaaCCCTTTCTATTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGCTTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggcgctgtttccattccagatcccggcgggggtgctgtcaaccttttttggcgcgccatattttctctatttgttgagaaagcagagccgttaatttttccctataacggcactgttgacatggccggcgaacacatt                                                                            >  2:451345‑M1/4‑78 (MQ=255)                                                                                                                                                                            aaGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtt                                                                            <  1:468254‑M1/251‑174 (MQ=255)                                                                                                                                                                               atatTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttac                                                                         <  1:446505‑M1/251‑177 (MQ=255)                                                                                                                                                                                  ttAACCCTTTCTTTTCATCTGGTTGTTTNTTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttac                                                                         <  2:251066‑M1/248‑177 (MQ=255)                                                                                                                                                                                  ttAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttac                                                                         >  1:251066‑M1/1‑72 (MQ=255)                                                                                                                                                                                       ccTTTCTTTTCATCTGTTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaac                                                                 <  1:193706‑M1/251‑185 (MQ=255)                                                                                                                                                                                                tCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgccc                                                         <  1:330469‑M1/250‑193 (MQ=255)                                                                                                                                                                                                  aTCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgat                                                      >  1:110519‑M1/1‑56 (MQ=255)                                                                                                                                                                                                       gttgttTATTAACCCTTCAGGAACACTCAGATTGCGTACCGCTTGCGATCCttacttctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaatcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgac                                                 <  1:296436‑M1/251‑201 (MQ=255)                                                                                                                                                                                                           tttATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctggcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtaagcccggacggataaagcgtttacgcgacaacaggtgaccgatgcggagctg                                             >  1:4513‑M1/1‑47 (MQ=255)                                                                                                                                                                                                              attaACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcgtcgcgccatattttatctatttgttgggaaagcagagacggtaatttttccataaaacggccctgtacgcctggacggataaaacgtttacgcgacaacggacgcccggtgcgacgctggcg                                          >  2:101515‑M1/1‑44 (MQ=255)                                                                                                                                                                                                              attaACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcg                                          >  1:225755‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                  aCCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttctctatttgttgggaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcggcaacggatgcccgatgcgacgctggcgcggc                                      >  1:216602‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                   cccTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctattttttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgaccgatgcgacgctggcgcgtct                                     >  2:414095‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                   cccTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaaacagagccggtaatttttccatataacggccctgtaagcctggacggatcaagcgtttacgcgacaacggatacccgaggcgacgcaggcgcgtct                                     >  1:273643‑M1/1‑39 (MQ=255)                                                                                                                                                                                                                                 tCAGATTGGGTACCGATTGCGAACCttactgctggggttcgctgattggtgtgggcggatggtgctgtttccattcaagatcccggcggggatgctgtcaacctttatcggcgagccatattttatctatttgttgagaaagcagagccgtttatttttacatataaaggccctgtacgcctggaaggataaagcgtttaagcgaaaacggatgcccgatgcgacgcaggagcgtgttattgggcctaaa                        >  2:180641‑M1/1‑25 (MQ=255)                                                                                                                                                                                                                                  cAGATTGCGTACCGCTTGAGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtccacctttatcggcgcgccatattttttctatttgtttagaaagcatcgccgttaatttttccatctaacgggcctgtaagcctggccgggtaaagcgtttacccgaccacggatgcccgattcgacgctggcgcgtcctttcgggccgccca                       >  2:442186‑M1/1‑24 (MQ=255)                                                                                                                                                                                                                                          gTACCTCTTGCGAATCgtactgcgggtgttcgctgactggggtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcttacaaacggccccg              <  2:349709‑M1/251‑236 (MQ=255)                                                                                                                                                                                                                                                   gCGAACCtttctgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccccgaatcgtag      <  1:70508‑M1/250‑244 (MQ=255)                                                                                                                                                                                                                                                     gacccttcctgcgggtgtgcgctggctggtgggggcggatggtgcggtttccattccagatcccggcggggctgcggtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccccgaatcgtaggtc   <  2:216602‑M1/248‑247 (MQ=255)                                                                                                                                                                                                                                                      atccttactgcgggtgttctctgtttggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcatcctttatcggcgcgccatattttatctatttgttgagaaagcatagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccccgaatcgtaggtcg  <  1:524877‑M1/249‑248 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                        TTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATG  >  NC_000913/167152‑167647

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: ‑