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

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

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

GCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCG  >  NC_000913/173067‑173562                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       ttgatcagtcttcccagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaacct                                                                                                                                                                                                                                                        <  2:394478‑M2/1‑1 (MQ=255)           ttctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGgtg                                                                                                                                                                                                                                             >  1:185435‑M2/240‑251 (MQ=255)            tctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGgtgt                                                                                                                                                                                                                                            >  2:492557‑M2/239‑251 (MQ=255)                 gccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCt                                                                                                                                                                                                                                       <  1:521070‑M2/18‑1 (MQ=255)                       gcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTgg                                                                                                                                                                                                                                 >  1:131550‑M2/228‑251 (MQ=255)                       gcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTgg                                                                                                                                                                                                                                 >  1:89022‑M2/228‑251 (MQ=255)                          ccagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGgtgt                                                                                                                                                                                                                              <  2:89022‑M2/27‑1 (MQ=255)                              ccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGCGGGc                                                                                                                                                                                                                          <  1:504225‑M2/31‑1 (MQ=255)                                             gctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGtt                                                                                                                                                                                                             <  1:233923‑M2/44‑1 (MQ=255)                                                  gcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGttt                                                                                                                                                                                                            <  1:213115‑M2/45‑1 (MQ=255)                                                  gcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGttt                                                                                                                                                                                                            >  2:213115‑M2/201‑245 (MQ=255)                                                     agagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAg                                                                                                                                                                                                   <  1:151490‑M2/54‑1 (MQ=255)                                                         cagaaagacagcgacggtgtagccggttggatcaaggatgtgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaagccttcaggaacgctcagattgcggaccgcttgcgaaccTTACTGATGGTGTTCGCTGACTGGTGTGGGAGGGTGGTGCTGTTTCCATTCCAGATcc                                                                                                                                                                                               >  1:195579‑M2/194‑251 (MQ=255)                                                              agacagcgacggtgtagccggttggctcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatacggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaaccctccaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGTTGTTCGCTGACTGGTGTGGGCGGATGGTGATGTTTCCATTCCAGCTCCcggcg                                                                                                                                                                                          >  1:459916‑M2/189‑251 (MQ=255)                                                                  agcgacggtgtagccggttggatcaaggctatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGAGGGCGGATGGTGCTGTTTCCCTTCCAGATCCCGGcggggg                                                                                                                                                                                      >  2:298139‑M2/185‑250 (MQ=255)                                                                  agcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgc                                                                                                                                                                                      >  1:533785‑M2/185‑251 (MQ=255)                                                                         gtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTc                                                                                                                                                                               <  1:155726‑M2/74‑1 (MQ=255)                                                                                cggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCttt                                                                                                                                                                        <  1:17122‑M2/81‑1 (MQ=255)                                                                                     ggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCg                                                                                                                                                                    <  1:492557‑M2/85‑1 (MQ=255)                                                                                        tcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGgcgc                                                                                                                                                                >  2:300727‑M2/163‑251 (MQ=255)                                                                                           aggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttggttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGcc                                                                                                                                                             <  1:66104‑M2/92‑1 (MQ=255)                                                                                                atgtttggtagtaattaacatctaagagtgaaataccggatggcgagttgccatccggaaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGCGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTTGCCACCGTTATCGGCGCGTCATAtt                                                                                                                                                        >  2:344768‑M2/155‑251 (MQ=255)                                                                                                  gttgggtagtaattaacatccaatcgtgaantgccggatggcgggttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttttttattaacccttcaggaacgctcagattgcgtacctcttgcgaaccTTACTGCTGGTGTGCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCTGCGCGCCATAttt                                                                                                                                                       <  2:239141‑M2/98‑1 (MQ=255)                                                                                                    ttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcgtgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTa                                                                                                                                                     <  1:49002‑M2/100‑1 (MQ=255)                                                                                                             attaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTc                                                                                                                                                                               <  1:295352‑M2/74‑1 (MQ=255)                                                                                                             attaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTc                                                                                                                                                                               >  2:295352‑M2/142‑215 (MQ=255)                                                                                                                                    cggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTTGGAAAGAAGAGCCGGTAAtttt                                                                                                                    >  2:167659‑M2/119‑251 (MQ=255)                                                                                                                                     ggatggcgagttgccatccggtaaaataacctcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAttttt                                                                                                                   <  1:192346‑M2/134‑1 (MQ=255)                                                                                                                                       atggcgagttgccatccgttaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTcc                                                                                                                 <  2:47134‑M2/136‑1 (MQ=255)                                                                                                                                       atggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTcc                                                                                                                 <  1:140017‑M2/136‑1 (MQ=255)                                                                                                                                                    atccggtaaaataacatcccatcaaatatatnaaccctttcttttcatctggttgtttattatcccttcaggaacgctcagattgcgtaccgctttcgaaccTTACTGCTTGTGTTCGCTGACTGGTGTTGGCGGATGGTTCTGTTTCCATTCCAGATCCCGGCGGGTCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTGTTTGTTGATAAAGCAGAGCCGTTAATTGTTCCATATAACGGCCCt                                                                                                    <  2:250500‑M2/149‑1 (MQ=255)                                                                                                                                                         gtaaaataacatcccatctaagatattaaccctttcttttcatctggctgtttataaacccttaaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACATTTCTCCGCGCGCCCTATTTTATCTCTTTTATGAGGAAACAGAGCCGTTAATTTTTCCATATAACGGGCCCGTACg                                                                                               >  2:482301‑M2/98‑251 (MQ=255)                                                                                                                                                           aaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTTTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGcc                                                                                             >  1:338234‑M2/96‑251 (MQ=255)                                                                                                                                                            aaataccatcccatctaagatattaaccctttcttttaatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgcaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttccattc                                                                                                                                                                                                      >  2:421616‑M2/95‑145 (MQ=255)                                                                                                                                                            aaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTttccattc                                                                                                                                                                                                      <  1:421616‑M2/51‑1 (MQ=255)                                                                                                                                                            aaaaaacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCt                                                                                            <  2:338234‑M2/157‑1 (MQ=255)                                                                                                                                                              ataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGGGAAAGCAGAGCcgttaatttttccatataagatcgggagagcacacc                                                                                          >  1:350860‑M2/93‑234 (MQ=255)                                                                                                                                                              ataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCAtata                                                                                                            <  2:350860‑M2/141‑1 (MQ=255)                                                                                                                                                                    tcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGGTGCTGTCAACCTTTTACGGCGCGCCATATTTTAACTTTTTGTTGAGAAAGCAGAGGCGTTAATTTTTCCATATAACGGACCGGTACGCCTGGacgggt                                                                                    >  1:205411‑M2/87‑249 (MQ=255)                                                                                                                                                                            aagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGtt                                                                            <  1:468254‑M2/173‑1 (MQ=255)                                                                                                                                                                            aaaatatgaaccctttctattcatctggttgtttattaacccttcaggaacgctcagcttgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGCGCTGTTTCCATTCCAGATCCCGGCGGGGGTGCTGTCAACCTTTTTTGGCGCGCCATATTTTCTCTATTTGTTGAGAAAGCAGAGCCGTTAatttttccctataacggcactgttgacatggccggcgaacacatt                                                                            >  2:451345‑M2/79‑229 (MQ=255)                                                                                                                                                                               atattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc                                                                         <  1:446505‑M2/176‑1 (MQ=255)                                                                                                                                                                                  ttaaccctttcttttcatctggttgtttnttaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc                                                                         <  2:251066‑M2/176‑1 (MQ=255)                                                                                                                                                                                  ttaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTAc                                                                         >  1:251066‑M2/73‑248 (MQ=255)                                                                                                                                                                                       cctttcttttcatctgtttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAAc                                                                 <  1:193706‑M2/184‑1 (MQ=255)                                                                                                                                                                                                tcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGccc                                                         <  1:330469‑M2/192‑1 (MQ=255)                                                                                                                                                                                                  atctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGAt                                                      >  1:110519‑M2/57‑251 (MQ=255)                                                                                                                                                                                                       gttgtttattaacccttcaggaacactcagattgcgtaccgcttgcgatccTTACTTCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAATCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGAc                                                 <  1:296436‑M2/200‑1 (MQ=255)                                                                                                                                                                                                           tttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGGCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTAAGCCCGGACGGATAAAGCGTTTACGCGACAACAGGTGACcgatgcggagctg                                             >  1:4513‑M2/48‑245 (MQ=255)                                                                                                                                                                                                              attaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGTCGCGCCATATTTTATCTATTTGTTGGGAAAGCAGAGACGGTAATTTTTCCATAAAACGGCCCTGTACGCCTGGACGGATAAAACGTTTACGCGACAACGGACGCCCGGTGCGACGCTGgcg                                          >  2:101515‑M2/45‑251 (MQ=255)                                                                                                                                                                                                              attaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGgcg                                          >  1:225755‑M2/45‑251 (MQ=255)                                                                                                                                                                                                                  acccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTCTCTATTTGTTGGGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGGCAACGGATGCCCGATGCGACGCTggcgcggc                                      >  1:216602‑M2/41‑249 (MQ=255)                                                                                                                                                                                                                   cccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTTTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGACCGATGCGACGCTGGCGCGTCt                                     >  2:414095‑M2/40‑251 (MQ=255)                                                                                                                                                                                                                   cccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAACAGAGCCGGTAATTTTTCCATATAACGGCCCTGTAAGCCTGGACGGATCAAGCGTTTACGCGACAACGGATACCCGAGGCGACGCAGGCGCGTCt                                     >  1:273643‑M2/40‑251 (MQ=255)                                                                                                                                                                                                                                 tcagattgggtaccgattgcgaaccTTACTGCTGGGGTTCGCTGATTGGTGTGGGCGGATGGTGCTGTTTCCATTCAAGATCCCGGCGGGGATGCTGTCAACCTTTATCGGCGAGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTTATTTTTACATATAAAGGCCCTGTACGCCTGGAAGGATAAAGCGTTTAAGCGAAAACGGATGCCCGATGCGACGCAGGAGCGTGTTATTGGGCctaaa                        >  2:180641‑M2/26‑248 (MQ=255)                                                                                                                                                                                                                                  cagattgcgtaccgcttgagaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCCACCTTTATCGGCGCGCCATATTTTTTCTATTTGTTTAGAAAGCATCGCCGTTAATTTTTCCATCTAACGGGCCTGTAAGCCTGGCCGGGTAAAGCGTTTACCCGACCACGGATGCCCGATTCGACGCTGGCGCGTCCTTtcgggccgccca                       >  2:442186‑M2/25‑245 (MQ=255)                                                                                                                                                                                                                                          gtacctcttgcgaatcGTACTGCGGGTGTTCGCTGACTGGGGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCTTACAAACGGCCCCg              <  2:349709‑M2/235‑1 (MQ=255)                                                                                                                                                                                                                                                   gcgaaccTTTCTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAg      <  1:70508‑M2/243‑1 (MQ=255)                                                                                                                                                                                                                                                     gacccTTCCTGCGGGTGTGCGCTGGCTGGTGGGGGCGGATGGTGCGGTTTCCATTCCAGATCCCGGCGGGGCTGCGGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTc   <  2:216602‑M2/246‑1 (MQ=255)                                                                                                                                                                                                                                                      atccTTACTGCGGGTGTTCTCTGTTTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCATCCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCATAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCg  <  1:524877‑M2/247‑1 (MQ=255)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       GCGATTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGACAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTCG  >  NC_000913/173067‑173562

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