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

GATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAT  >  NC_000913/173069‑173565                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          gatcagtcttctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattacccctttcttttcatctggttgtttattaacccttcaggaangctcagattgcgtaccgcttgcgaaccTTAc                                                                                                                                                                                                                                                        >  2:254014‑M2/248‑251 (MQ=255) gatcagacttctcagccgcagcaccagccgcaaaagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggctggctcacggatatgattggtagtaattaacatctacgcgtgaaataccggatggcgagttgccatccggaaaaataacatcccatctaagatattaaccctttcttttcatctggatgtttattaccccttcaggaacgctacgattgcgtaccgctagcgaaccTTAc                                                                                                                                                                                                                                                        >  2:189988‑M2/248‑251 (MQ=255) aatcagtctcctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccctccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTAc                                                                                                                                                                                                                                                        <  2:73043‑M2/4‑1 (MQ=255)  atcagtcttctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtatcgcttgcgaaccttcct                                                                                                                                                                                                                                                       >  1:201886‑M2/247‑248 (MQ=255)          tctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGgt                                                                                                                                                                                                                                                 <  1:58102‑M2/11‑1 (MQ=255)                  tcagcaacagccgcaacagcaacctgaacagcaagagcagaaagacaacgacggagtagccggttggatccaggatatgtttggtagtaattaacatcgaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTactggtgttggttgac                                                                                                                                                                                                                                       >  2:352074‑M2/231‑244 (MQ=255)                     gcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTgg                                                                                                                                                                                                                                    >  1:240799‑M2/228‑251 (MQ=255)                       aacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGgtg                                                                                                                                                                                                                                  >  1:182717‑M2/226‑251 (MQ=255)                        acagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGgtgt                                                                                                                                                                                                                                 >  1:370249‑M2/225‑251 (MQ=255)                             cgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaccatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccatttcttttcatctggttgtttattaacccttcaggaacgctcagattgagtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGtgtggccg                                                                                                                                                                                                                            >  2:386597‑M2/220‑248 (MQ=255)                               caacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCgg                                                                                                                                                                                                                           >  1:138082‑M2/218‑250 (MQ=255)                                aacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTANTGCTGGTGTTCGCTGACTGGTGTGGGCgg                                                                                                                                                                                                                           >  2:294568‑M2/217‑249 (MQ=255)                                      aacctgctcagcaagcgcagaaagaccgcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataacggatggcgagttgccatccggtaaaataacatcacatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTTCTGGTGTTCAATGACTGGGGGGGGCGGATGGGGCt                                                                                                                                                                                                                   >  2:354657‑M2/211‑251 (MQ=255)                                           gctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGttt                                                                                                                                                                                                               <  1:265500‑M2/45‑1 (MQ=255)                                            ctcagcaagagcagaaaggcagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttcc                                                                                                                                                                                                             >  2:478129‑M2/205‑251 (MQ=255)                                              cagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaatacaggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacacttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttccat                                                                                                                                                                                                           >  1:39183‑M2/203‑251 (MQ=255)                                                 caagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttccattcc                                                                                                                                                                                                        <  2:370249‑M2/52‑1 (MQ=255)                                                         gaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCcg                                                                                                                                                                                                >  1:133395‑M2/192‑251 (MQ=255)                                                         gaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCcg                                                                                                                                                                                                <  1:41742‑M2/60‑1 (MQ=255)                                                          aaagacagcgacggagtagccggttggatcacggatatgtttggtagtacttaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagctattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCcgg                                                                                                                                                                                               >  2:366072‑M2/191‑251 (MQ=255)                                                            agacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaacccttttttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGGGCTGTTTCCATTCCAGATCCcggcg                                                                                                                                                                                             >  1:215112‑M2/189‑251 (MQ=255)                                                                agcgacggtgtagccggttggatcaaggatatgtttggcagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttaaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCTCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgc                                                                                                                                                                                         >  2:171677‑M2/185‑251 (MQ=255)                                                                   gacggtgtagccggttggatcaaggctatgtttggtagtaattaacatctaagcgtgaaataccggaaggcgagttgccatccggtaaaataacatcacatctaagacattaaccctttcttttcatctggttgtatattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGCTGGTGCTGTTTCAATTCCAGATACCGGCGGCgctgc                                                                                                                                                                                      >  1:356465‑M2/182‑251 (MQ=255)                                                                           agccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAAcc                                                                                                                                                                              <  1:528387‑M2/78‑1 (MQ=255)                                                                             ccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCtt                                                                                                                                                                            <  1:91282‑M2/80‑1 (MQ=255)                                                                              cgggtggatcaaggatatgtttgggagtaattaacatctaagcgtggaatcccggatggcgagttgccagccggtaaaataacatcccatctaagatattagccctttcttttcatcttgttgcttattaaccctccaggaacgctcagattgggtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGTATCGTGCTGTTTCCTTTCCAGATCCCGGCGGGGCTGCTGTcaaccgtg                                                                                                                                                                           <  2:20131‑M2/81‑4 (MQ=255)                                                                                  tggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGTTGTGGGCGGTTGGTGCTGTTTCCATTCCAGATCCCGGCGGAGCGGCTGTCAGCCTTtattt                                                                                                                                                                       >  2:383739‑M2/167‑249 (MQ=255)                                                                                        aaggatatgtttggtagtaattaacctctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCNGCGGGGCTGCCGGCAACCTTTATCGgctcgc                                                                                                                                                                 >  2:287875‑M2/161‑251 (MQ=255)                                                                                         aggatatggttggtagtaattaacatctaagcgtgaaatcccgtatggcgagttgcaatccggtaaaataacatcccatctaagatattaaccctttcgtttcatctggttgtttatttacccttcaggaacgctcagattgcgtcccgcgtgcgaaccTTACTGCGGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGcc                                                                                                                                                                <  2:21680‑M2/92‑1 (MQ=255)                                                                                              atgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCatat                                                                                                                                                            <  1:530020‑M2/96‑1 (MQ=255)                                                                                                    ggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATc                                                                                                                                                      <  1:539238‑M2/102‑1 (MQ=255)                                                                                                          aagtaacatctaagcgtgaaataccggatggctagttgccatccggtaaaataacatcacatctaagatattaacccttacttttcatctggttgtggattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGGTCGCTGACTGGTGTGGGCGGAGGGTGCTGTTTCCATTCCAGATCCCGTCGGGGCTGCTGTCAACCTGTATCGGCTCGCAATATTTTAACTATttg                                                                                                                                                <  2:39183‑M2/108‑1 (MQ=255)                                                                                                                           gaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCg                                                                                                                               <  1:285658‑M2/125‑1 (MQ=255)                                                                                                                           gaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGt                                                                                                                              <  1:440189‑M2/126‑1 (MQ=255)                                                                                                                            aaataccggagggcgagttgccatccggtaaaataacatcccatctaagatattaacccgtgcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtacctcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGGTTCCATTCCAGATCCCGGAGGGGCTGCTGTCAACCGTGTGCGGCGCGCCATATTTTAGCTATTTGTTGAGAAAGAAGAGCCGtt                                                                                                                             <  2:92802‑M2/127‑1 (MQ=255)                                                                                                                              ataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTaa                                                                                                                           <  1:254014‑M2/129‑1 (MQ=255)                                                                                                                                   ggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAttttt                                                                                                                      <  1:376310‑M2/134‑1 (MQ=255)                                                                                                                                     atggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCCTATTTTATCTATTTGTTGAGAAAGAAGAGCCGTTAATTTTTcc                                                                                                                    >  1:562651‑M2/116‑251 (MQ=255)                                                                                                                                     atggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTcc                                                                                                                    >  1:3746‑M2/116‑251 (MQ=255)                                                                                                                                         cgagttgccatccggtaaaataacagcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCatat                                                                                                                <  1:248126‑M2/140‑1 (MQ=255)                                                                                                                                               gccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAATCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACgg                                                                                                           <  1:386597‑M2/145‑1 (MQ=255)                                                                                                                                                           aataccatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCCCGCCCTATTTTATCTATTTGTTGAGAAAGCAAAACCGTTAATTTTTTCATATAAACGCCCTGTAcgcccg                                                                                              >  2:231181‑M2/94‑249 (MQ=255)                                                                                                                                                               acatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACg                                                                                          <  1:389969‑M2/162‑1 (MQ=255)                                                                                                                                                                catcccatctaagctattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAAACTTTATCTGCCCGCCCTATTTTATCTATTTGTTGAGTAAGACGAGCAGCTAATTTTTCCATATAACGGCCCTGTACGCCTAGACgg                                                                                         >  2:542841‑M2/89‑251 (MQ=255)                                                                                                                                                                 atcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGAcggg                                                                                        >  1:443992‑M2/88‑250 (MQ=255)                                                                                                                                                                     catctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATaaa                                                                                    <  1:478129‑M2/168‑1 (MQ=255)                                                                                                                                                                     catctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATaa                                                                                     >  1:397983‑M2/84‑250 (MQ=255)                                                                                                                                                                             atattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagatttcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCGGCTGTCAAACTTTATCGGCGCGCCATATTTTTTCTCTTTGTTGAGAAAGCAGAGCCGTTCATTTTTCCATATAACAGCCCTGTACGCCTGGACGGGTAAAGCGTTTAc                                                                            >  2:392495‑M2/76‑251 (MQ=255)                                                                                                                                                                                         tcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAtata                                                                                                               >  1:246673‑M2/64‑204 (MQ=255)                                                                                                                                                                                         tcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAGTTTTCCAtata                                                                                                               <  2:246673‑M2/141‑1 (MQ=255)                                                                                                                                                                                            tttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGAt                                                                <  1:242790‑M2/188‑1 (MQ=255)                                                                                                                                                                                              tcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTAAACCTTTATCGGCGCGCCATCTTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCCTATAACGGACCTGTACGCCTGGACGGATAAAGCGTTTACGCGATCACGGATGCCCg                                                           >  1:423585‑M2/59‑251 (MQ=255)                                                                                                                                                                                               catctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCg                                                           <  2:138082‑M2/193‑1 (MQ=255)                                                                                                                                                                                                    tgttgtttattaacccttcaggaacgctcagattgcttaccgcttgcgaaccTTACCGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTCCGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGc                                                       <  2:517897‑M2/197‑1 (MQ=255)                                                                                                                                                                                                    ggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGc                                                       <  1:403856‑M2/197‑1 (MQ=255)                                                                                                                                                                                                    ggttgtttattaacccttcaggaaagctcagattgcgtaccgcttgcgaaccGTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGAt                                                         <  1:256583‑M2/195‑1 (MQ=255)                                                                                                                                                                                                     gttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGTTGTGGGCGGATGGTGCTGTGTCCATTCCAGTTCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGAc                                                    <  2:312505‑M2/200‑1 (MQ=255)                                                                                                                                                                                                           tattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTTTTGAGAAAGCAGAACCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGAGACAACGGATGCCCGATGCGACGCCGgc                                              >  1:560955‑M2/46‑251 (MQ=255)                                                                                                                                                                                                            attaacccttcaggaacgctcagattgcgtaccgcttacgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTTTCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTAAGACTGGACGGCTAAAGCGTTTACGCGACCACCGATGCCCGATGCGACGCTGgcg                                             >  1:444702‑M2/45‑251 (MQ=255)                                                                                                                                                                                                                acccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGGTTACGCGACAACGGATGCCCGATGCGAAGCTggcgcggc                                         >  1:555279‑M2/41‑249 (MQ=255)                                                                                                                                                                                                                      caggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATc                                    <  1:287875‑M2/216‑1 (MQ=255)                                                                                                                                                                                                                       aggatcgcgcagattgagttccgcttgcgaaccTTCCTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCgg                                  <  2:17694‑M2/218‑1 (MQ=255)                                                                                                                                                                                                                             gctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTAc                            >  1:414916‑M2/28‑251 (MQ=255)                                                                                                                                                                                                                              ctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGCCGCTGGCGCGTCTTATCGGGCCTAc                            >  1:101137‑M2/27‑250 (MQ=255)                                                                                                                                                                                                                                     tgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGGCTTTTCGGGCCTACAAACGGcc                    >  2:463616‑M2/20‑251 (MQ=255)                                                                                                                                                                                                                                       cgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGGGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGTTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTGTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGcccc                  <  2:182717‑M2/234‑1 (MQ=255)                                                                                                                                                                                                                                         taccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCCGCGCGCCATATTTTATCTATTTTTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTTGCGCGTCTTATCGGGCCTACAAACGGCCCCg                 >  1:339102‑M2/16‑250 (MQ=255)                                                                                                                                                                                                                                                tgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTa          <  1:266340‑M2/242‑1 (MQ=255)                                                                                                                                                                                                                                                ggcgatccTTACTGCTGGTGTTCGCTGACTGGGGGGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGTCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAATCGGCCCCGAATCGt           <  2:53302‑M2/241‑1 (MQ=255)                                                                                                                                                                                                                                                   gaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCTGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGTGTCTTTTCGGGCCTACAAACGGCCCCGCATCGtagggc      >  1:400593‑M2/6‑249 (MQ=255)                                                                                                                                                                                                                                                    aaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTTTTTGTTGAGAAAACAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCCTAGGTCg     >  1:216370‑M2/5‑251 (MQ=255)                                                                                                                                                                                                                                                     accTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCgg    <  1:503059‑M2/248‑1 (MQ=255)                                                                                                                                                                                                                                                        ttACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGATGCTGTCAACCTTTATAGGCGCGCCATATTTTTTCTCTTTGTTGATAAAGCAGAGCCGTTAATTTTTCCCTATAACGGCCCTTTACGCCCGGACGGATAAAACGTTTACGCGACAACGGATGCCCGATGCCCCGCTGGCGCCTCTTATCATGACCACAAACGGccccccaacctagggctgac  >  2:430486/1‑234 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                          GATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCGGAT  >  NC_000913/173069‑173565

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