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	72 [0]	[0] 73	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)	72 (0.990)	67/492	0.0	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

GATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGG  >  NC_000913/167154‑167648                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         gATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTACCCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAANGCTCAGATTGCGTACCGCTTGCGAACCttac                                                                                                                                                                                                                                                      >  2:254014‑M1/1‑247 (MQ=255) gATCAGACTTCTCAGCCGCAGCACCAGCCGCAAAAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGCTGGCTCACGGATATGATTGGTAGTAATTAACATCTACGCGTGAAATACCGGATGGCGAGTTGCCATCCGGAAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGATGTTTATTACCCCTTCAGGAACGCTACGATTGCGTACCGCTAGCGAACCttac                                                                                                                                                                                                                                                      >  2:189988‑M1/1‑247 (MQ=255) aatCAGTCTCCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCCTCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttac                                                                                                                                                                                                                                                      <  2:73043‑M1/250‑5 (MQ=255)  aTCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTATCGCTTGCGAACCttcct                                                                                                                                                                                                                                                     >  1:201886‑M1/1‑246 (MQ=255)          tctcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggt                                                                                                                                                                                                                                               <  1:58102‑M1/249‑12 (MQ=255)                  tcagcaaCAGCCGCAACAGCAACCTGAACAGCAAGAGCAGAAAGACAACGACGGAGTAGCCGGTTGGATCCAGGATATGTTTGGTAGTAATTAACATCGAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactactggtgttggttgac                                                                                                                                                                                                                                     >  2:352074‑M1/2‑230 (MQ=255)                     gcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactgg                                                                                                                                                                                                                                  >  1:240799‑M1/1‑227 (MQ=255)                       aaCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtg                                                                                                                                                                                                                                >  1:182717‑M1/1‑225 (MQ=255)                        aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgt                                                                                                                                                                                                                               >  1:370249‑M1/1‑224 (MQ=255)                             cGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTACCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCATTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGAGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtggccg                                                                                                                                                                                                                          >  2:386597‑M1/1‑219 (MQ=255)                               cAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcgg                                                                                                                                                                                                                         >  1:138082‑M1/1‑217 (MQ=255)                                aaCAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttantgctggtgttcgctgactggtgtgggcgg                                                                                                                                                                                                                         >  2:294568‑M1/1‑216 (MQ=255)                                      aaCCTGCTCAGCAAGCGCAGAAAGACCGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATAACGGATGGCGAGTTGCCATCCGGTAAAATAACATCACATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttacttctggtgttcaatgactggggggggcggatggggct                                                                                                                                                                                                                 >  2:354657‑M1/1‑210 (MQ=255)                                           gCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgttt                                                                                                                                                                                                             <  1:265500‑M1/250‑46 (MQ=255)                                            cTCAGCAAGAGCAGAAAGGCAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttcc                                                                                                                                                                                                           >  2:478129‑M1/1‑204 (MQ=255)                                              cAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACAGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACACTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccat                                                                                                                                                                                                         >  1:39183‑M1/1‑202 (MQ=255)                                                 cAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattcc                                                                                                                                                                                                      <  2:370249‑M1/251‑53 (MQ=255)                                                         gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccg                                                                                                                                                                                              >  1:133395‑M1/1‑191 (MQ=255)                                                         gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccg                                                                                                                                                                                              <  1:41742‑M1/251‑61 (MQ=255)                                                          aaaGACAGCGACGGAGTAGCCGGTTGGATCACGGATATGTTTGGTAGTACTTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGCTATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccgg                                                                                                                                                                                             >  2:366072‑M1/1‑190 (MQ=255)                                                            aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTTTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggggctgtttccattccagatcccggcg                                                                                                                                                                                           >  1:215112‑M1/1‑188 (MQ=255)                                                                aGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGCAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTAAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttctctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggc                                                                                                                                                                                       >  2:171677‑M1/1‑184 (MQ=255)                                                                   gACGGTGTAGCCGGTTGGATCAAGGCTATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGAAGGCGAGTTGCCATCCGGTAAAATAACATCACATCTAAGACATTAACCCTTTCTTTTCATCTGGTTGTATATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggctggtgctgtttcaattccagataccggcggcgctgc                                                                                                                                                                                    >  1:356465‑M1/1‑181 (MQ=255)                                                                           aGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacc                                                                                                                                                                            <  1:528387‑M1/251‑79 (MQ=255)                                                                             ccGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctt                                                                                                                                                                          <  1:91282‑M1/251‑81 (MQ=255)                                                                              cgggtggatCAAGGATATGTTTGGGAGTAATTAACATCTAAGCGTGGAATCCCGGATGGCGAGTTGCCAGCCGGTAAAATAACATCCCATCTAAGATATTAGCCCTTTCTTTTCATCTTGTTGCTTATTAACCCTCCAGGAACGCTCAGATTGGGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcgtatcgtgctgtttcctttccagatcccggcggggctgctgtcaaccgtg                                                                                                                                                                         <  2:20131‑M1/247‑82 (MQ=255)                                                                                  tGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactgttgtgggcggttggtgctgtttccattccagatcccggcggagcggctgtcagcctttattt                                                                                                                                                                     >  2:383739‑M1/1‑166 (MQ=255)                                                                                        aaGGATATGTTTGGTAGTAATTAACCTCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccngcggggctgccggcaacctttatcggctcgc                                                                                                                                                               >  2:287875‑M1/1‑160 (MQ=255)                                                                                         aGGATATGGTTGGTAGTAATTAACATCTAAGCGTGAAATCCCGTATGGCGAGTTGCAATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCGTTTCATCTGGTTGTTTATTTACCCTTCAGGAACGCTCAGATTGCGTCCCGCGTGCGAACCttactgcgggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgcc                                                                                                                                                              <  2:21680‑M1/251‑93 (MQ=255)                                                                                              atGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatat                                                                                                                                                          <  1:530020‑M1/250‑97 (MQ=255)                                                                                                    ggTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatc                                                                                                                                                    <  1:539238‑M1/250‑103 (MQ=255)                                                                                                          aagtaacATCTAAGCGTGAAATACCGGATGGCTAGTTGCCATCCGGTAAAATAACATCACATCTAAGATATTAACCCTTACTTTTCATCTGGTTGTGGATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtggtcgctgactggtgtgggcggagggtgctgtttccattccagatcccgtcggggctgctgtcaacctgtatcggctcgcaatattttaactatttg                                                                                                                                              <  2:39183‑M1/247‑109 (MQ=255)                                                                                                                           gAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgt                                                                                                                            <  1:440189‑M1/251‑127 (MQ=255)                                                                                                                           gAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccg                                                                                                                             <  1:285658‑M1/250‑126 (MQ=255)                                                                                                                            aaaTACCGGAGGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCGTGCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCTCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctggttccattccagatcccggaggggctgctgtcaaccgtgtgcggcgcgccatattttagctatttgttgagaaagaagagccgtt                                                                                                                           <  2:92802‑M1/251‑128 (MQ=255)                                                                                                                              aTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaa                                                                                                                         <  1:254014‑M1/251‑130 (MQ=255)                                                                                                                                   ggATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaattttt                                                                                                                    <  1:376310‑M1/251‑135 (MQ=255)                                                                                                                                     aTGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccctattttatctatttgttgagaaagaagagccgttaatttttcc                                                                                                                  >  1:562651‑M1/1‑115 (MQ=255)                                                                                                                                     aTGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttcc                                                                                                                  >  1:3746‑M1/1‑115 (MQ=255)                                                                                                                                         cGAGTTGCCATCCGGTAAAATAACAGCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatat                                                                                                              <  1:248126‑M1/251‑141 (MQ=255)                                                                                                                                               gCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaatctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacgg                                                                                                         <  1:386597‑M1/250‑146 (MQ=255)                                                                                                                                                           aaTACCATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcccgccctattttatctatttgttgagaaagcaaaaccgttaattttttcatataaacgccctgtacgcccg                                                                                            >  2:231181‑M1/1‑93 (MQ=255)                                                                                                                                                               aCATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacg                                                                                        <  1:389969‑M1/251‑163 (MQ=255)                                                                                                                                                                cATCCCATCTAAGCTATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaactttatctgcccgccctattttatctatttgttgagtaagacgagcagctaatttttccatataacggccctgtacgcctagacgg                                                                                       >  2:542841‑M1/1‑88 (MQ=255)                                                                                                                                                                 aTCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggg                                                                                      >  1:443992‑M1/1‑87 (MQ=255)                                                                                                                                                                     cATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaa                                                                                  <  1:478129‑M1/251‑169 (MQ=255)                                                                                                                                                                     cATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataa                                                                                   >  1:397983‑M1/1‑83 (MQ=255)                                                                                                                                                                             atatTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTTCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggcggctgtcaaactttatcggcgcgccatattttttctctttgttgagaaagcagagccgttcatttttccatataacagccctgtacgcctggacgggtaaagcgtttac                                                                          >  2:392495‑M1/1‑75 (MQ=255)                                                                                                                                                                                         tctttTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatata                                                                                                             >  1:246673‑M1/1‑63 (MQ=255)                                                                                                                                                                                         tctttTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaagttttccatata                                                                                                             <  2:246673‑M1/204‑142 (MQ=255)                                                                                                                                                                                            tttCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggat                                                              <  1:242790‑M1/248‑189 (MQ=255)                                                                                                                                                                                              tCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtaaacctttatcggcgcgccatcttttatctatttgttgagaaagcagagccgttaatttttccctataacggacctgtacgcctggacggataaagcgtttacgcgatcacggatgcccg                                                         >  1:423585‑M1/1‑58 (MQ=255)                                                                                                                                                                                               cATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccg                                                         <  2:138082‑M1/250‑194 (MQ=255)                                                                                                                                                                                                    tgttgtttATTAACCCTTCAGGAACGCTCAGATTGCTTACCGCTTGCGAACCttaccgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtccgcctggacggataaagcgtttacgcgacaacggatgcccgatgc                                                     <  2:517897‑M1/248‑198 (MQ=255)                                                                                                                                                                                                    ggTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgc                                                     <  1:403856‑M1/249‑198 (MQ=255)                                                                                                                                                                                                    ggTTGTTTATTAACCCTTCAGGAAAGCTCAGATTGCGTACCGCTTGCGAACCgtactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgat                                                       <  1:256583‑M1/247‑196 (MQ=255)                                                                                                                                                                                                     gttgttTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactgttgtgggcggatggtgctgtgtccattccagttcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgac                                                  <  2:312505‑M1/251‑201 (MQ=255)                                                                                                                                                                                                           tattaACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctattttttgagaaagcagaaccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgagacaacggatgcccgatgcgacgccggc                                            >  1:560955‑M1/1‑45 (MQ=255)                                                                                                                                                                                                            attaACCCTTCAGGAACGCTCAGATTGCGTACCGCTTACGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttttctatttgttgagaaagcagagccgttaatttttccatataacggccctgtaagactggacggctaaagcgtttacgcgaccaccgatgcccgatgcgacgctggcg                                           >  1:444702‑M1/1‑44 (MQ=255)                                                                                                                                                                                                                aCCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcggttacgcgacaacggatgcccgatgcgaagctggcgcggc                                       >  1:555279‑M1/1‑40 (MQ=255)                                                                                                                                                                                                                      cAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatc                                  <  1:287875‑M1/250‑217 (MQ=255)                                                                                                                                                                                                                       aGGATCGCGCAGATTGAGTTCCGCTTGCGAACCttcctgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgg                                <  2:17694‑M1/251‑219 (MQ=255)                                                                                                                                                                                                                             gCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctac                          >  1:414916‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                              cTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgccgctggcgcgtcttatcgggcctac                          >  1:101137‑M1/1‑26 (MQ=255)                                                                                                                                                                                                                                     tGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcggcttttcgggcctacaaacggcc                  >  2:463616‑M1/1‑19 (MQ=255)                                                                                                                                                                                                                                       cGTACCGCTTGCGAACCttactgctggtgttcgctgactggggtgggcggatggtgctgtttccattccagatcccggcggggttgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtgtacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggcccc                <  2:182717‑M1/251‑235 (MQ=255)                                                                                                                                                                                                                                         tACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatccgcgcgccatattttatctattttttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgcttgcgcgtcttatcgggcctacaaacggccccg               >  1:339102‑M1/1‑15 (MQ=255)                                                                                                                                                                                                                                                tGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccccgaatcgta        <  1:266340‑M1/250‑243 (MQ=255)                                                                                                                                                                                                                                                ggcGATCCttactgctggtgttcgctgactggggggggcggatggtgctgtttccattccagatcccggcgggtctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaatcggccccgaatcgt         <  2:53302‑M1/248‑242 (MQ=255)                                                                                                                                                                                                                                                   gAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatctgcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgtgtcttttcgggcctacaaacggccccgcatcgtagggc    >  1:400593‑M1/1‑5 (MQ=255)                                                                                                                                                                                                                                                    aaCCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctttttgttgagaaaacagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccccgaatcctaggtcg   >  1:216370‑M1/1‑4 (MQ=255)                                                                                                                                                                                                                                                     aCCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccccgaatcgtaggtcgg  <  1:503059‑M1/251‑249 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                         GATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATGG  >  NC_000913/167154‑167648

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