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	56 [0]	[0] 56	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)	56 (0.860)	52/478	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

TTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTC  >  NC_000913/173071‑173561                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       tcagtcttctcagccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTActg                                                                                                                                                                                                                                                   >  1:308450‑M2/246‑251 (MQ=255)            agccgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGc                                                                                                                                                                                                                                        >  1:76242‑M2/235‑251 (MQ=255)               cgcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTg                                                                                                                                                                                                                                      >  1:418087‑M2/232‑250 (MQ=255)                gcagcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGAc                                                                                                                                                                                                                                    >  2:440025‑M2/231‑251 (MQ=255)                  agcaacagccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTg                                                                                                                                                                                                                                  <  1:477959‑M2/23‑1 (MQ=255)                     aacagccgcaacagcaacctgctcagcaagagaagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGgtg                                                                                                                                                                                                                               >  2:159856‑M2/226‑251 (MQ=255)                        agccgcaacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTgg                                                                                                                                                                                                                            >  2:190485‑M2/223‑251 (MQ=255)                              aacagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGAt                                                                                                                                                                                                                      >  2:427175‑M2/217‑251 (MQ=255)                               acagcaacctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATg                                                                                                                                                                                                                     <  1:495122‑M2/36‑1 (MQ=255)                                      cctgctcagcaagagcagaaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGt                                                                                                                                                                                                              <  1:16154‑M2/43‑1 (MQ=255)                                                       gaaagacagcgacggtgtagccggttggatcnaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctgtcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATGCCAGATCCcg                                                                                                                                                                                             <  2:252631‑M2/60‑1 (MQ=255)                                                        aaagacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCcgg                                                                                                                                                                                            >  2:341883‑M2/191‑251 (MQ=255)                                                          agacagcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccctctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCcggcg                                                                                                                                                                                          >  1:471920‑M2/189‑251 (MQ=255)                                                               gcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgct                                                                                                                                                                                     <  2:85186‑M2/68‑1 (MQ=255)                                                               gcgacggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGgct                                                                                                                                                                                     <  1:14390‑M2/68‑1 (MQ=255)                                                                   cggtgtagccggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGctcctg                                                                                                                                                                                 >  2:108271‑M2/180‑251 (MQ=255)                                                                            cggttggatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCttt                                                                                                                                                                        <  2:313520‑M2/81‑1 (MQ=255)                                                                                  gatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGgc                                                                                                                                                                  <  1:341883‑M2/87‑1 (MQ=255)                                                                                  gatcaaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGgc                                                                                                                                                                  <  2:7767‑M2/87‑1 (MQ=255)                                                                                      aaggatatgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCCGTCAACCTTTTTTGgcgcgc                                                                                                                                                              >  2:163274‑M2/161‑251 (MQ=255)                                                                                            atgtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAAACTTTATCGGCGCGCCATAtt                                                                                                                                                        >  1:479722‑M2/155‑251 (MQ=255)                                                                                              gtttggtagtaattaacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGTCATAtttt                                                                                                                                                      >  1:234841‑M2/153‑251 (MQ=255)                                                                                                            aacatctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaccccttccttttcatcttgttgtttattaacccttcaggaacgcccagattgcgttccgcttgcgaaccTTACTGCTGGTGTTCGCTGCCTGGTGTGGGCGGATTTTGCTGTTTCCGTTCCAGATCCCGGCGGAGCTACTGTCAACCTTTATCGGCGCGTCATATTTCATCTATTTGTTga                                                                                                                                         <  2:442921‑M2/112‑1 (MQ=255)                                                                                                              catctaagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTTTTTGTTgaga                                                                                                                                       >  2:120015‑M2/137‑250 (MQ=255)                                                                                                                  taagcgtgaaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTTTCGGCGCGCCATATTTTATCTATTTGTTGAGGAAGc                                                                                                                                   >  2:320493‑M2/133‑250 (MQ=255)                                                                                                                       gtgaaataccggatggcgagttgccctccggtaaaataacatcccatctatgatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTTTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGATAGCAGAGcc                                                                                                                             <  2:479722‑M2/124‑1 (MQ=255)                                                                                                                          aaataccggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGtt                                                                                                                          <  1:159856‑M2/127‑1 (MQ=255)                                                                                                                                cggatggcgagttgccatccggtaaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAAAGCCGTTAAtttt                                                                                                                    >  1:99291‑M2/119‑251 (MQ=255)                                                                                                                                                  ccggtaaaataacatcccatctaagatattacccctttcttttcatcgggttgtttattaacccttcaggaacgctcagattgcgttccgctggcgaaccTTTCGGCTGGGGTTCGCTGCCTGGTGTGGGCGGTGGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAGTTTTCCATATAACGGCCCTg                                                                                                   <  2:169024‑M2/150‑1 (MQ=255)                                                                                                                                                    gttaaaataacatctcatctaagattttaaccctttcttttcatctggttggttattaacccttgaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCGGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCAGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTAc                                                                                                <  2:480441‑M2/153‑1 (MQ=255)                                                                                                                                                      taaaataacatcccatctaagatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTTTTGAGAAAGCAGACCCGTTAATTTTTCCATATAACTGCCCTGTCCGc                                                                                              >  2:138242‑M2/97‑251 (MQ=255)                                                                                                                                                                    atctaagatattaaccctttcttttcatctggttgtttattaacccttcaggatcgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAAACTTTACAGGCGCGCCATATTTTATCTATTTGTTGAGAAAGAAGAGTCGTTAATTTTTCCATATAACGGCACTGACCGCCAGGACGGATaa                                                                                  >  2:441163‑M2/83‑249 (MQ=255)                                                                                                                                                                       tacgatattaaccctttcttttcatctggttgtttattaacccttaaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCCTTCCAGATCCCGGCGGGGCTGCTGTAGACACTTATTGGCGCGCCATATTTTATCTTTTTGTTGCGAAAACAAAGCCGTTAATTTTtccctataaaggccctgtacaccatgacgagcaaagcgt                                                                             >  2:477934‑M2/80‑232 (MQ=255)                                                                                                                                                                          gatattaaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTCATTTTTCCATATAACGGCCCTTTTCTCCTGGACGGGTAAAGCgtttt                                                                          >  2:56040‑M2/77‑250 (MQ=255)                                                                                                                                                                               taaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGGTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTAAGCCTGGACGGATAAAACGTttacgagc                                                                     >  1:287948‑M2/72‑248 (MQ=255)                                                                                                                                                                                aaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGa                                                                     <  1:120015‑M2/180‑1 (MQ=255)                                                                                                                                                                                aaccctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑c                                                                   <  1:258207‑M2/181‑1 (MQ=255)                                                                                                                                                                                   cctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCCGCGCGCCATATTTTATCTTTTTGTTTAGAAAACAGAACCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGcgtttacgagaaca                                                                  >  2:105554‑M2/68‑244 (MQ=255)                                                                                                                                                                                    ctttcttttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgagaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAACAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCCGGACGGATAAAGCGTTTACGCGA‑CAACg                                                               >  1:151207‑M2/67‑251 (MQ=255)                                                                                                                                                                                          tttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGATg                                                           <  2:117622‑M2/189‑1 (MQ=255)                                                                                                                                                                                           ttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgcaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATTTTTTTTCTATTTGTTGAGAAAGCAGGACCGTTAATTTTTCCATATAAAGGCCCTGTACCCCTGGACGGGTAAAGCGTTTTCGCGA‑CAACGGATTccc                                                        >  1:218108‑M2/60‑251 (MQ=255)                                                                                                                                                                                           ttcatctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCGGCTGTCAACCTTTATCGGCGCGCCATATTTTTTCTATTTTTTGTGAAAACAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGGCGGATAAATGGTTTACGCGACCAAC‑GATGccc                                                        >  2:90547‑M2/60‑251 (MQ=255)                                                                                                                                                                                             catctggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCATAGCCGTTAATTTTTCCATATAACGGCCCGGTACGCCTGGACGGATAAAGCGTTTACTCGA‑CAACGGAAGCCCg                                                       >  1:345878‑M2/58‑250 (MQ=255)                                                                                                                                                                                                  ggttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGATGCCCGAt                                                     <  1:440025‑M2/195‑1 (MQ=255)                                                                                                                                                                                                   gttgtttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGATGCCCGATGCGAc                                                <  1:56040‑M2/200‑1 (MQ=255)                                                                                                                                                                                                       tttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCCTATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTAAGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGATGCCCGCTGCGACGCTg                                            >  1:71471‑M2/48‑251 (MQ=255)                                                                                                                                                                                                        ttattaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGATGCCCGATGCGACGCTg                                            <  1:11531‑M2/204‑1 (MQ=255)                                                                                                                                                                                                           ttaacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAAATTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGCTGCCCGATGCGACGCTTgcgc                                        >  1:215783‑M2/44‑251 (MQ=255)                                                                                                                                                                                                            taacccttcaggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGATGCCCGATGCGACGCTGgcgc                                        <  1:284724‑M2/208‑1 (MQ=255)                                                                                                                                                                                                                   tcaggaacgctcagattgcgtaccgcttgctaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGCGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGACGTTAATTTTTCCATATAACGGCCCTGTACGCCTGTCAGGATAAAGCGATTACGCGA‑AAACGGATTCCCGCAGCGCCGCTGTCGCGTCTTATc                                >  2:498278‑M2/36‑251 (MQ=255)                                                                                                                                                                                                                    caggaacgctcagattgcgtaccgcttgctaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGTCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGGTTACGCGA‑CCAAGGATGCCCGCTGCGCCGCTGGCGCGTCTTATCg                               >  2:330523‑M2/35‑251 (MQ=255)                                                                                                                                                                                                                      ggaacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGGTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGATGCCCGATGCGACGCTGGCGCGTCTtatccgg                             >  2:379205‑M2/33‑248 (MQ=255)                                                                                                                                                                                                                        aacgctcagattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGAAGGATAAAGCGTTTACGCGA‑CAACGGATGCCCGATGCGCCGCTGGCGCGTCTTATCGGGcc                           >  1:86909‑M2/31‑251 (MQ=255)                                                                                                                                                                                                                          cgctccgattgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTTCTGTCAACCTTTATCGGCGCGCCATATTTTATCTTTTTGTTGAGAAATCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACCCCTGGGCCGATAAAAAGTTTACGCCA‑CAACGGATGCCCGATGCGACGCTGGCGCGTCTTttcgggact                          >  1:93178‑M2/29‑247 (MQ=255)                                                                                                                                                                                                                                 attgcgtaccgcttgcgaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCCTATAACGGCCCTGTACGCCTGGACGGAAAAAGCGTTTACGCGA‑CAACGGATGCCCGATGCGACGCTGGCGCGGCTTATCGGGCCTACAAACgg                  >  1:106064‑M2/22‑251 (MQ=255)                                                                                                                                                                                                                                                 gaaccTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAGCCTTTCTCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTTCATATAACGGCCCTGTACGCCTGGGCGGATAAAGCGTTTACGCGA‑CAACCGATGCCCGATGCGACGCTGGCCCGTCTTATCGGgcctacaaacggcccccccacgtaaggc  >  1:3573‑M2/6‑239 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       TTACCCCGCTGTGCCGCCGCTGGCTGACCATTTTACCGCTGGGTGGTGATACCGCCCGAGCCGTAGGAATGGCGCTGACGCCGACGCGAATTGCGCTGCTGCTGTTAGCGGCTTGCCTGACGGCGACCGCGACGATGACTATTGGACCGTTGAGTTTTGTTGGTTTAATGGCACCGCATATTGCGCGGATGATGGGCTTTCGACGGACGATGCCACACATCGTAATTTCGGCGCTGGTGGGTGGTTTACTGCTGGTGTTCGCTGACTGGTGTGGGCGGATGGTGCTGTTTCCATTCCAGATCCCGGCGGGGCTGCTGTCAACCTTTATCGGCGCGCCATATTTTATCTATTTGTTGAGAAAGCAGAGCCGTTAATTTTTCCATATAACGGCCCTGTACGCCTGGACGGATAAAGCGTTTACGCGA‑CAACGGATGCCCGATGCGACGCTGGCGCGTCTTATCGGGCCTACAAACGGCCCCGAATCGTAGGTC  >  NC_000913/173071‑173561

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