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	68 [0]	[0] 69	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)	68 (0.690)	66/492	0.2	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

TGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATG  >  NC_000913/167153‑167647                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        tGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCtt                                                                                                                                                                                                                                                       <  1:302265‑M1/250‑3 (MQ=255)  gcccagtCTGCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGAATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGAGATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCTTACCGCTTGCGAACCttac                                                                                                                                                                                                                                                     <  2:571884‑M1/248‑5 (MQ=255)  gATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttac                                                                                                                                                                                                                                                     >  1:49253‑M1/1‑247 (MQ=255)       gTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTAAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCACGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGACCGCTCAGATTGCGTACAGCTTGCGAACCttactgctg                                                                                                                                                                                                                                                >  2:452651‑M1/1‑242 (MQ=255)               aGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgc                                                                                                                                                                                                                                        <  1:347701‑M1/251‑18 (MQ=255)                gCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTGCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgct                                                                                                                                                                                                                                       <  1:715635‑M1/251‑19 (MQ=255)                 ccGCAGCAACAGCCGCAACAGCAACCGGCTCNGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCGTGCGAACCttactgctggtgttcgctg                                                                                                                                                                                                                                      <  2:374261‑M1/251‑20 (MQ=255)                         aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgt                                                                                                                                                                                                                              <  1:387526‑M1/251‑28 (MQ=255)                               gCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcgg                                                                                                                                                                                                                        >  2:92900‑M1/1‑218 (MQ=255)                               gCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcgg                                                                                                                                                                                                                        <  1:387769‑M1/251‑34 (MQ=255)                                          cTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTCCATCTGGTTGTTTCTTAACCATTCAGGAACGCTCAGATTGCGAACCGCTTGCGACCCttactgctggtgttcgctgactggtgagggcggatggggccgtt                                                                                                                                                                                                             >  2:402177‑M1/1‑207 (MQ=255)                                             cTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTTGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttcc                                                                                                                                                                                                          >  1:214882‑M1/1‑204 (MQ=255)                                                 gCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattc                                                                                                                                                                                                      >  1:585014‑M1/1‑200 (MQ=255)                                                     gagCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccaga                                                                                                                                                                                                  <  2:710366‑M1/251‑56 (MQ=255)                                                         aGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatccc                                                                                                                                                                                              >  1:67556‑M1/1‑192 (MQ=255)                                                               aCAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggg                                                                                                                                                                                        <  1:236119‑M1/251‑66 (MQ=255)                                                                      cGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctg                                                                                                                                                                                 >  1:549906‑M1/1‑179 (MQ=255)                                                                       ggTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgt                                                                                                                                                                                <  1:257987‑M1/251‑74 (MQ=255)                                                                          gtAGCCGGTTGGAGCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggctgggctgctgtca                                                                                                                                                                              <  2:735058‑M1/250‑76 (MQ=255)                                                                                ggTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaccttta                                                                                                                                                                       <  1:685239‑M1/251‑83 (MQ=255)                                                                                    ggATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcgg                                                                                                                                                                   <  1:92900‑M1/251‑87 (MQ=255)                                                                                     gATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggc                                                                                                                                                                  <  1:529433‑M1/251‑88 (MQ=255)                                                                                        cAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcg                                                                                                                                                               >  1:699461‑M1/1‑161 (MQ=255)                                                                                          aGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtccaccattatcggcgcgcc                                                                                                                                                             >  2:731104‑M1/1‑159 (MQ=255)                                                                                           ggATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcatgtcggc                                                                                                                                                             >  1:490622‑M1/1‑158 (MQ=255)                                                                                           ggATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatc                                                                                                                                                                     <  2:490622‑M1/242‑85 (MQ=255)                                                                                            gATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccat                                                                                                                                                           <  2:699461‑M1/251‑95 (MQ=255)                                                                                                 gTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaccttttttggcccgccctttttt                                                                                                                                                      >  2:498754‑M1/1‑152 (MQ=255)                                                                                                 gTTTGGGAGTAATTAACATCTAATCGTGAAAGACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTACTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgacgggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctttcaacctttatcggcgcgccatattt                                                                                                                                                       <  2:67556‑M1/250‑99 (MQ=255)                                                                                                     ggTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatct                                                                                                                                                  <  2:689789‑M1/251‑104 (MQ=255)                                                                                                       tagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatctgcgcgccatattttatctat                                                                                                                                                >  1:694779‑M1/1‑146 (MQ=255)                                                                                                          taATTAACACCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctattt                                                                                                                                              <  2:178041‑M1/250‑108 (MQ=255)                                                                                                           aaTTAACATCGAAGCGTGAAACACCGGATGGGGAGGTGCCATCCGGGAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCAGCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttcaattccagatccctgcggggctgctgtcaacctttatcggcgcgccatattttatctatttgt                                                                                                                                            <  2:202379‑M1/251‑110 (MQ=255)                                                                                                            aTTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaaccnttatcggcgcgccatattttatctatttgtt                                                                                                                                           >  2:386086‑M1/1‑141 (MQ=255)                                                                                                              tAACATCTAATCGTGAAATACCGGATGGCGAGTTGACATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGCTTGCGTACCGCTTGCGAACCttactgctggcgttcgctggctggggtgggcggatggtgctgtttccattccagatcccggcgggggtgctgtcaagatttattggggcgccctatcagccgggataggggc                                                                                                                                         >  2:603047‑M1/1‑139 (MQ=255)                                                                                                              tAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccatgccagatcccggcggggctgctgtcaacctttatcggcgcgccatat                                                                                                                                                         <  1:603047‑M1/235‑97 (MQ=255)                                                                                                                 cATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaa                                                                                                                                      <  2:127281‑M1/251‑116 (MQ=255)                                                                                                                        gCGTGAACTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgcctggtgtgggcggatggtgctgtttccattccagatcccggcgggggcgccgtcaatctttatcggcgcgcaatcttttttctaatttatgagaaagcag                                                                                                                                 >  2:108231‑M1/1‑129 (MQ=255)                                                                                                                             aaaTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgtt                                                                                                                          <  1:669885‑M1/251‑128 (MQ=255)                                                                                                                                 aCCGGATGGCGCGTTGCCATCCGGTAAAATAACATCACATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgacttgtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatatgttatcgagttgttgagaaagcagggccgttaatt                                                                                                                      <  2:226573‑M1/251‑132 (MQ=255)                                                                                                                                     gATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctggtaacctttatctgggcgccatattttttcttttttttgagaaagcagagccgttaatttttt                                                                                                                  >  2:70529‑M1/1‑116 (MQ=255)                                                                                                                                        ggCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccat                                                                                                               <  1:314745‑M1/251‑139 (MQ=255)                                                                                                                                                  cATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccctattttatctatttgttgagaaagcagagccgttcatttttccctataacggccc                                                                                                     >  2:318774‑M1/1‑103 (MQ=255)                                                                                                                                                   aTCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatacgcgcgccatattttatctatttgtttagaaagcagagccgttaatttttctctataacggccct                                                                                                    >  2:688710‑M1/1‑102 (MQ=255)                                                                                                                                                            aaTAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTTGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctg                                                                                           <  1:446093‑M1/251‑159 (MQ=255)                                                                                                                                                               aaCATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcaaagccgttaatttttccatataacggccctgtacgcctggac                                                                                        >  1:482705‑M1/1‑90 (MQ=255)                                                                                                                                                                   ccccaTCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggat                                                                                    <  2:585014‑M1/250‑166 (MQ=255)                                                                                                                                                                             gATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttttcggcgcgccatattttttctattttttgagaaagcagagccgttaatttttccatataaaggccctgtacgcctggacggataaagcgtttt                                                                          >  1:694044‑M1/1‑76 (MQ=255)                                                                                                                                                                             gATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacg                                                                                                         <  1:237906‑M1/220‑145 (MQ=255)                                                                                                                                                                             gATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatccgcgcgccatattttttctatttgttgagaaagcagagccggtaatttttccatataaag                                                                                                         >  2:237906‑M1/1‑76 (MQ=255)                                                                                                                                                                                atTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgcctggtgtgggcggatggtgctgtttccattccagatcccggcggggctgccgtcaacctttatccgcgcgccctattttatctatttggtgagaaagcagcgccgttaatttttccatataacgggcctgtacccctcgacggataaaacgtttaccc                                                                       >  1:182523‑M1/1‑73 (MQ=255)                                                                                                                                                                                atTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacg                                                                        <  1:504346‑M1/250‑178 (MQ=255)                                                                                                                                                                                  tAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcg                                                                      <  1:575950‑M1/250‑180 (MQ=255)                                                                                                                                                                                      ccTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagaaccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaac                                                                 >  1:696168‑M1/1‑67 (MQ=255)                                                                                                                                                                                       ctttctttTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcatggacgggtaaagcgtttacgcgaaaacg                                                                >  2:731038‑M1/1‑66 (MQ=255)                                                                                                                                                                                               tCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccg                                                        >  1:576538‑M1/1‑58 (MQ=255)                                                                                                                                                                                                cATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccga                                                       <  1:452651‑M1/251‑195 (MQ=255)                                                                                                                                                                                                 aTCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacgggtaaagcgtttacgcgacaacggatgcccgat                                                      >  1:98859‑M1/1‑56 (MQ=255)                                                                                                                                                                                                           ttattaACCCTTCAGGAACGCACAGATTGCGTACCGCTTGCTAACCttacttctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggcccggtacgcctggacggataaagcgtttactcgacaacggatgcccgatgcgacgctg                                             <  2:694779‑M1/250‑205 (MQ=255)                                                                                                                                                                                                             attaACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggc                                           <  1:386086‑M1/250‑207 (MQ=255)                                                                                                                                                                                                               taACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatacccgatgcgacgctggcgcg                                        >  1:690408‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                              gCTCAGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctctttgttgagaaagcagagccggtaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcggcttatcggggctcc                         >  2:78925‑M1/1‑27 (MQ=255)                                                                                                                                                                                                                                  aGATTGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaac                     >  1:485048‑M1/1‑23 (MQ=255)                                                                                                                                                                                                                                     ttGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctttacgcctggacggataaagcgtttacgcgccaacggatgcccgatgcgaccctggcgcgtcttatcgggccaccaaacggc                  >  1:10560‑M1/1‑20 (MQ=255)                                                                                                                                                                                                                                     ttGCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaaccgc                  >  1:475690‑M1/1‑20 (MQ=255)                                                                                                                                                                                                                                       gCGTACCGCTTGCGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccatattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaatggatgcccgatgcgacgctggcgcgtcttatcgggcctacaaacggccc                <  1:433049‑M1/251‑234 (MQ=255)                                                                                                                                                                                                                                                   cGAACCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgccattttttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtccgcctggacggataaagcgtttacgcgacaacgggtgcccgatgcgacgctggcgcgtcttatcgggccttaaaacggcccccaatcgtaggg    >  2:757365‑M1/1‑6 (MQ=255)                                                                                                                                                                                                                                                     aaCCttactgctggtgttcgctgactggtgtgggcggatggtgctgtttccattccagatcccggcggggctgctgtcaacctttatcggcgcgcccaattttatctatttgttgagaaagcagagccgttaatttttccatataacggccctgtacgcctggacggataaagcgtttacgcgacaacgggtggccgatgcgacgctggcgcgtcttaacgggcccacaaaaggccccgaatcgtaggtcg  >  1:97261‑M1/1‑4 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                        TGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTGCGTAAAATCGCAGTTGTAGTAGCCACAGCGGTTAGCGGCATGTCTGTTTATGCACAGGCAGCGGTTGAACCGAAAGAAGACACTATCACCGTTACCGCTGCACCTGCGCCGCAAGAAAGCGCATG  >  NC_000913/167153‑167647

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