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

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

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

GCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAT  >  NC_000913/173067‑173565                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            ttgatcagtcttctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaacctt                                                                                                                                                                                                                                                            >  1:318722‑M2/250‑251 (MQ=255)                 gccgcggcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcgtgcgaaccTTACTGCTGGTGTTCGCt                                                                                                                                                                                                                                            <  1:560622‑M2/18‑1 (MQ=255)                 gccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCt                                                                                                                                                                                                                                            <  1:74268‑M2/18‑1 (MQ=255)                   cgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTg                                                                                                                                                                                                                                           >  1:31145‑M2/232‑250 (MQ=255)                    gcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGAc                                                                                                                                                                                                                                         <  1:118987‑M2/21‑1 (MQ=255)                    gcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGAc                                                                                                                                                                                                                                         >  1:429612‑M2/231‑251 (MQ=255)                              ccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGc                                                                                                                                                                                                                               <  2:508392‑M2/31‑1 (MQ=255)                                     agcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcttgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattacccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACGGCTGGTGTTCGCTGACTGGTGTGGGCGGAtggc                                                                                                                                                                                                                        <  2:142728‑M2/38‑2 (MQ=255)                                             gctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttc                                                                                                                                                                                                                <  2:481978‑M2/46‑1 (MQ=255)                                                 agcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttccatt                                                                                                                                                                                                            <  1:547397‑M2/50‑1 (MQ=255)                                                  gcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttccattc                                                                                                                                                                                                           <  2:539835‑M2/51‑1 (MQ=255)                                                         cagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATcc                                                                                                                                                                                                    >  1:525806‑M2/194‑251 (MQ=255)                                                          agaaagacagcgacggggtagccggttggatcaaggagatgtttggtagtaattaacatcgaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATccc                                                                                                                                                                                                   <  2:140200‑M2/59‑1 (MQ=255)                                                                  agcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgc                                                                                                                                                                                           <  1:317933‑M2/67‑1 (MQ=255)                                                                  agcgacggtgaagacggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggctggcgagttgccatccggtaaactaacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggcacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGTAGCTGTTTCCATTCCAGATCCCGGcggggg                                                                                                                                                                                           >  2:101941‑M2/185‑250 (MQ=255)                                                                   gcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgct                                                                                                                                                                                          <  1:310493‑M2/68‑1 (MQ=255)                                                                   gcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgct                                                                                                                                                                                          <  1:208707‑M2/68‑1 (MQ=255)                                                                     gacggtgtagccggttggatcaaggatatgtttggtagtaattaccatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCggggcctc                                                                                                                                                                                        >  2:259307‑M2/182‑248 (MQ=255)                                                                      acggtgtagccggttggataaaggatatgtttggtagtaattaacatctaagcgtgacataccggatggcgagttgccatccggtaaaataacatcccatctaagatcttaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgctgct                                                                                                                                                                                       >  1:96652‑M2/181‑251 (MQ=255)                                                                       cggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGctgctg                                                                                                                                                                                      <  2:429612‑M2/72‑1 (MQ=255)                                                                        ggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgctgctgg                                                                                                                                                                                     >  2:153397‑M2/179‑250 (MQ=255)                                                                                 ggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTTAACctttt                                                                                                                                                                            >  1:141672‑M2/170‑250 (MQ=255)                                                                                 ggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGGGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACctttt                                                                                                                                                                            >  1:570273‑M2/170‑250 (MQ=255)                                                                                  gttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagctattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAAcctttct                                                                                                                                                                           >  1:655181‑M2/169‑249 (MQ=255)                                                                                  gttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTAAAACTTTAt                                                                                                                                                                           >  1:107017‑M2/169‑251 (MQ=255)                                                                                         caaggatatgtttggcagtaattaacatctaagcgtgaaataccggatggcgagttgcaatccggtaaaaaaacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttaaggaacgctcagattgcgtaccccttgcgaaccTTACTGCTGGTGTTCGCTGCCTGGTGGGGGCGATTGGTGCTGTTTCCTTTCAAGATCCGGCCGGGGCTGTTGTCCACCTTTATcggcggc                                                                                                                                                                    >  2:68763‑M2/162‑249 (MQ=255)                                                                                           aggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATTGGCTCGcc                                                                                                                                                                  >  1:692592‑M2/160‑251 (MQ=255)                                                                                                         agtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttattttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTACAGATCCCGGCGGGGCTGGTGGCCACCTTTATCTGCCCGCCCTATTTTatctttt                                                                                                                                                    >  1:291044‑M2/146‑248 (MQ=255)                                                                                                           taattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttc                                                                                                                                                                                                                >  2:453575‑M2/144‑189 (MQ=255)                                                                                                           taattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttc                                                                                                                                                                                                                <  1:453575‑M2/46‑1 (MQ=255)                                                                                                              ttaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCCTATTTTATCTAtttgttt                                                                                                                                               >  1:197898‑M2/141‑250 (MQ=255)                                                                                                                       aagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCag                                                                                                                                      <  1:615773‑M2/120‑1 (MQ=255)                                                                                                                         gcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCagag                                                                                                                                    <  2:379311‑M2/122‑1 (MQ=255)                                                                                                                          cgtgaaataccggatggcgtgttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaccgctcagattgcgttccgcttgcgaaccTTACTGCTGGTGTTCTCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGc                                                                                                                                   <  1:481218‑M2/123‑1 (MQ=255)                                                                                                                             gaaataccggatggccagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCGGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑Gt                                                                                                                               <  2:394940‑M2/126‑1 (MQ=255)                                                                                                                                    cggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGCGCTGCTGTCAACACTTATCGGCGCGCCATATTTTATCTAATTTTTGGGAAAGCAGAGCC‑GTTAAtttt                                                                                                                        >  2:625239‑M2/119‑251 (MQ=255)                                                                                                                                     ggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttggttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGCGCC‑GTTAAttttt                                                                                                                       <  1:495412‑M2/134‑1 (MQ=255)                                                                                                                                         ggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaaccattcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGAGGGTGCTGTTTCCACTCCAGATCCCGGCGGGGCTGCTGTGCACCTTTATAGTCGCCCCGTATTTTATCTCTTTGTTGAGAAAGCCAAGCC‑GTAAATCTTTCCat                                                                                                                   >  2:724191‑M2/114‑251 (MQ=255)                                                                                                                                           cgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcagctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGATAAAGCAGGGCC‑GTTAATTTTTCCatat                                                                                                                 <  2:635745‑M2/140‑1 (MQ=255)                                                                                                                                            gagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCAtata                                                                                                                <  1:466694‑M2/141‑1 (MQ=255)                                                                                                                                             agtggccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATaa                                                                                                               <  2:306040‑M2/142‑1 (MQ=255)                                                                                                                                                tgccatccggtaaaataatatcccatctaagatattaacccttccttttcatctggttgtttattaacccttcaggaacgataagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGCGCTGTTTCCATTCCAGTTCCCGGCGGGGCTGCTGTCAACCTTTCTCGGCGCCCCACATTTTATCTATTTTTTTAGAAAGCCGAGCC‑GGTAAtttttccattaaaatg                                                                                                            >  2:615252‑M2/107‑244 (MQ=255)                                                                                                                                                tgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACg                                                                                                             <  1:567440‑M2/144‑1 (MQ=255)                                                                                                                                                 gccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCATCCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGc                                                                                                           <  1:568721‑M2/146‑1 (MQ=255)                                                                                                                                                  ccatcccgttaacgatcatcccatctaagaaattaaccctttcttttcaactggttttttattaacccttcaggatccctcagattgcgtcccgcgtgcgaaccTTACTGCTGGTGTTCGCGGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCACCCTGTATCGGCGCGCCAGATTTTATCTGTTTGTTGAGAAAGCAGCGCC‑GTTAATTTTTCCATATGACGGcc                                                                                                          <  2:291044‑M2/147‑1 (MQ=255)                                                                                                                                                     tccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCCTTTTTTTTCTATTTTTTGGGAAAGCAGAACC‑GTTAATTTTTCCTTTTAAAGGCCc‑tt                                                                                                      >  2:180699‑M2/102‑250 (MQ=255)                                                                                                                                                        ggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTAc                                                                                                   <  1:152275‑M2/153‑1 (MQ=255)                                                                                                                                                                aacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGTGAAAGCAGAGCC‑GTTAATTTTTCCATATTACGGCCC‑TGCAAGCCTGGAc                                                                                           >  2:134599‑M2/91‑251 (MQ=255)                                                                                                                                                                  catcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGAGCGCCATATTTTATCTATTTGGTGGGAAAGCAGAGCC‑GTTAATTTTTTCCTATAACGGCCc‑tttacgcccggcagg                                                                                         >  2:434045‑M2/89‑244 (MQ=255)                                                                                                                                                                   atcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCGTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GGTAATTTTTCCATATAAAGGCCC‑TGTACGCCTGGAcggg                                                                                        >  1:472058‑M2/88‑250 (MQ=255)                                                                                                                                                                          ctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAACC‑GGTAATTTTTCCATATAAAGGCCC‑TGTACGCCTGGACGGATAAAg                                                                                   >  1:222912‑M2/81‑249 (MQ=255)                                                                                                                                                                            aagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGt                                                                                <  1:386290‑M2/172‑1 (MQ=255)                                                                                                                                                                              cccctcttccgatctttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTa                                                                             <  2:138291‑M2/175‑1 (MQ=255)                                                                                                                                                                               atattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTAc                                                                            <  1:153397‑M2/176‑1 (MQ=255)                                                                                                                                                                                 attaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTTAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTAcgc                                                                          >  1:629764‑M2/74‑251 (MQ=255)                                                                                                                                                                                      ccctttcttttcatctggttgtttattaacccttcaggaacgctcaggttgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTTCGCGACa                                                                      >  1:128214‑M2/69‑250 (MQ=255)                                                                                                                                                                                       cctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCCGCTGTCACCCTTTATCGGCGCGGCATATTTTTTCTATTTGTTGAGAATACAGAGCCGGTTAA‑TTTTCCATATAAAGGCCCTTCTA‑gcctggaccgataaaggtgttctgcgacaac                                                                    >  1:393242‑M2/68‑236 (MQ=255)                                                                                                                                                                                          ttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTTTCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGCCGGATAAAGCGTTTACGCGACAAcggc                                                                 >  2:77758‑M2/65‑250 (MQ=255)                                                                                                                                                                                             ttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGGGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTa                                                                             >  1:138291‑M2/62‑236 (MQ=255)                                                                                                                                                                                              tttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGCAAACGGATGcc                                                             >  1:391712‑M2/61‑251 (MQ=255)                                                                                                                                                                                                 catctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCg                                                           >  1:615313‑M2/58‑250 (MQ=255)                                                                                                                                                                                                 catcgggttgtttattaaccctcaaggaacgctcagatggcgtaccgcttgcgacccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGTATGCCCg                                                           <  2:699857‑M2/193‑1 (MQ=255)                                                                                                                                                                                                      ggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTAACCGACAACGGATgaccgaagccc                                                     >  2:227350‑M2/53‑246 (MQ=255)                                                                                                                                                                                                              attaacccgtcaggaacgctcagattgcgtaccgattgcgacccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGGTGCTGTCAACCTTTATCGGCGCGCCATATTTTTTCTATTTGTTGTGAAAGCCGAGCC‑GTTAATTTTTCCCTTTAACGGCCC‑TGTAGCCCTGGACGGATAAAGCGTTTACGCGAcaacggatgcccgagtcggagctggcg                                             >  2:210161‑M2/45‑238 (MQ=255)                                                                                                                                                                                                                 atcccttcatgaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCGGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGcgcg                                           <  2:75260‑M2/209‑1 (MQ=255)                                                                                                                                                                                                                 aacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGgcgcgg                                          >  1:146693‑M2/42‑250 (MQ=255)                                                                                                                                                                                                                  acccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGGGAAAGCAGAGCC‑GTTAATTTTTCCATATAACAGCCC‑TGTACGCCTGGACGGCTGACACGGTTAAGCGTCAACGGAAGCCCCATGCGACGCAGGAGCGTc                                         >  1:372344‑M2/41‑251 (MQ=255)                                                                                                                                                                                                                     cttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAAAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTcttt                                      >  1:588581‑M2/38‑250 (MQ=255)                                                                                                                                                                                                                       tcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTa                                      <  1:487035‑M2/214‑1 (MQ=255)                                                                                                                                                                                                                       gcaggaacgctcagattgcgtcccgcttgcgaaccTTACGTCTGGTGTTCGCTGACTGGTGTGGGCGGCTGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCAGATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTtag                                     <  2:43599‑M2/215‑2 (MQ=255)                                                                                                                                                                                                                        caggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTAAGCCTGGACGGATAAAGCGGTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCg                                   >  1:559867‑M2/35‑251 (MQ=255)                                                                                                                                                                                                                         aggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGGCAACCTTTATCGGCGCGCCCTATTTTATCTTTTTTTTGAGGAAGCCGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGGACGCCTGGGCGGGAAAAGCGGTTACGCGACAACGGNTGCCCGATGCGGAGCCGGCGGGtcttaaccg                                  >  2:360462‑M2/34‑247 (MQ=255)                                                                                                                                                                                                                              cgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCCTATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTTCGCGACAACGGGTGCCNGACGCGCCGCCGGCGCCTGTTATTGggcccc                             >  2:346651‑M2/29‑249 (MQ=255)                                                                                                                                                                                                                                  caatttgcgtaccgcttgcgaaccTTACTGCTGGTTTTCGCTGACGGGTGTGGGCGTATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACaaa                         <  2:671713‑M2/227‑1 (MQ=255)                                                                                                                                                                                                                                                cttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAACAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCCACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGt           >  1:241307‑M2/11‑251 (MQ=255)                                                                                                                                                                                                                                                    cgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTaggg       >  1:493284‑M2/7‑250 (MQ=255)                                                                                                                                                                                                                                                          ttACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCCTATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGGTGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGCATCGGAGGGCGGAt  >  1:312125/1‑250 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                            GCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCC‑GTTAATTTTTCCATATAACGGCCC‑TGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAT  >  NC_000913/173067‑173565

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