Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,139,771	10483 bp→TG		[ptsA]–[eptC]	[ptsA], frwC, frwB, pflD, pflC, frwD, yijO, [eptC]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4139771	4150253	10483	76 [0]	[0] 74	[ptsA]–[eptC]	[ptsA],frwC,frwB,pflD,pflC,frwD,yijO,[eptC]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4139770	0 (0.000)	74 (0.760) +TG	67/488	0.1	100%	coding (2452/2502 nt)	ptsA	putative PTS multiphosphoryl transfer protein PtsA
	?	NC_000913	4150254 =	0 (0.000)					coding (12/1734 nt)	eptC	phosphoethanolamine transferase EptC

CGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACG  >  NC_000913/4139521‑4140014                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     cGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCt                                                                                                                                                                                                                                                     >  1:606005‑M1/1‑250 (MQ=255)  ggCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGAGAATGCGCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTGCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtg                                                                                                                                                                                                                                                    <  2:290303‑M1/251‑3 (MQ=255)              gAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcat                                                                                                                                                                                                                                        <  2:658572‑M1/251‑15 (MQ=255)                 cAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTTTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagt                                                                                                                                                                                                                                     <  1:738461‑M1/251‑18 (MQ=255)                 cAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagt                                                                                                                                                                                                                                     <  1:190010‑M1/251‑18 (MQ=255)                  aaaCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTCCGGTTCAGGTTCATGTTGCGGC‑GGCTTGGCGTTAGCGTCGGCAGGGGCtgtgtggaatgcatagt                                                                                                                                                                                                                                     >  2:596055‑M1/1‑233 (MQ=255)                            cccacagcaccaatcgCGTTTTGCAGACCCGGCAGCACTTCCCATATGGACCTCTTTCTGGCACCGATAATGCTCGGGTTTGTTGGCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGTGATAGCGGCGGCAGAGGCtgtgtggaatgcatagtaggaccctgtc                                                                                                                                                                                                                          <  2:590170‑M1/244‑29 (MQ=255)                            cccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtc                                                                                                                                                                                                                          >  1:96955‑M1/1‑223 (MQ=255)                            cccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtc                                                                                                                                                                                                                          <  1:437464‑M1/251‑29 (MQ=255)                             cctcctCACCAATCGCTTTTTGCAGACGCGGAAGCACTTCCCATATGGCCTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCTATCGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTATCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGAtgtgtggaatgcatagtaagaacctgtct                                                                                                                                                                                                                         <  1:336886‑M1/251‑30 (MQ=255)                             ccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtct                                                                                                                                                                                                                         <  1:602462‑M1/251‑30 (MQ=255)                                catccCCAATCGCTTTTTGCAGACGCGGCAGCTCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGCGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtctt                                                                                                                                                                                                                        <  2:560349‑M1/249‑31 (MQ=255)                                    cccaATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagatcctgtcttgaaaaa                                                                                                                                                                                                                  <  1:88738‑M1/250‑37 (MQ=255)                                     ccAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaa                                                                                                                                                                                                                 <  2:212196‑M1/251‑38 (MQ=255)                                       aaTCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaata                                                                                                                                                                                                               <  1:256339‑M1/251‑40 (MQ=255)                                          cgcgtttagcagacgcgGCAGCACGTCCCATATGGTCTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTGTTACAGTTCCAGTTCATGTTGCAGCAGGG‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcataggaagaacctggcttgaaaaaatatcg                                                                                                                                                                                                            <  1:228834‑M1/243‑43 (MQ=255)                                           gcattttgcaGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGTCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgc                                                                                                                                                                                                           <  2:221440‑M1/248‑44 (MQ=255)                                              ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccg                                                                                                                                                                                                         >  1:523543‑M1/1‑205 (MQ=255)                                              ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGGTAGTTGTCACCCAGGAAATGGGGAATCTGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGAGACAGTTCCATGATGTTTCCTTTTACAATTCCAGTTCATGTTGCAGCCGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcctngtaagaacctgtcttgaaaaaatatcgccg                                                                                                                                                                                                         >  2:395047‑M1/1‑205 (MQ=255)                                                    aGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATTCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaa                                                                                                                                                                                                  <  2:625313‑M1/251‑53 (MQ=255)                                                    aGACGCGGCAGCACTTCCAATATGGACTCCTTGCTGGCCGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggtatgcatagaaagaacctgtcttgaaaaaatatcgccgaatgta                                                                                                                                                                                                   >  2:612853‑M1/1‑199 (MQ=255)                                                     gACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaac                                                                                                                                                                                                 >  1:730571‑M1/1‑198 (MQ=255)                                                            cagcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaatt                                                                                                                                                                                          <  1:23974‑M1/251‑61 (MQ=255)                                                             agccagtcccagaTGGACCCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattc                                                                                                                                                                                         <  2:410988‑M1/245‑62 (MQ=255)                                                              gcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTGCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattcc                                                                                                                                                                                        <  1:232532‑M1/251‑63 (MQ=255)                                                                aCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacggcaattcctt                                                                                                                                                                                      <  1:217760‑M1/251‑65 (MQ=255)                                                                   tCCCATAAGGACTCCTGGCTGGCAGCGATAATGCTCGTGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGAGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtgtaatgcatagtaagaacctgtctttaaaaaatatcgccgaatgtaacgacaattccttaa                                                                                                                                                                                    <  2:562971‑M1/250‑67 (MQ=255)                                                                     ccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaagg                                                                                                                                                                                  >  2:587803‑M1/1‑182 (MQ=255)                                                                      cATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggat                                                                                                                                                                                <  2:606005‑M1/251‑71 (MQ=255)                                                                           ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctg                                                                                                                                                                           >  2:38555‑M1/1‑176 (MQ=255)                                                                           ggACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctg                                                                                                                                                                           <  2:457822‑M1/251‑76 (MQ=255)                                                                                ccTTGGTGGCAGCGATAATGCTCGGGTTAGTTTGCACCCCGGCAATGGGGAATATGCGTGCCAGACCTTCGACTTCTGCGACGGTAGCGGTGTCCAGATACCGTTCCATTATTTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaatgatatctgaaggt                                                                                                                                                                      <  1:596055‑M1/251‑81 (MQ=255)                                                                                    gCTGGCAGCGATAATGCTCGGGTTAGTTTGNACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtata                                                                                                                                                                   <  2:389865‑M1/250‑84 (MQ=255)                                                                                    gCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatat                                                                                                                                                                  <  1:673707‑M1/251‑85 (MQ=255)                                                                                    gCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtata                                                                                                                                                                   <  1:574904‑M1/250‑84 (MQ=255)                                                                                         cAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcaga                                                                                                                                                             <  1:258670‑M1/251‑90 (MQ=255)                                                                                            cgagaatgcaCGGGTTAGTTGTCCCCCCGGCTATGGGGAATATGCGTGCAAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGAGACAGTTCCATGATGTTTCCTTTTGCAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtatcgacaattccttaaggatatctgaaggtatattcagaatt                                                                                                                                                          <  2:724448‑M1/247‑93 (MQ=255)                                                                                            cGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatt                                                                                                                                                          >  1:165769‑M1/1‑159 (MQ=255)                                                                                                  ggcTCGGGTTTGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaata                                                                                                                                                    <  2:705564‑M1/250‑99 (MQ=255)                                                                                                      cGGGTTAGTTGTCAACCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaat                                                                                                                                                <  2:278896‑M1/251‑103 (MQ=255)                                                                                                            aGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagac                                                                                                                                          <  2:91754‑M1/251‑109 (MQ=255)                                                                                                                     cccGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatata                                                                                                                                 <  1:118867‑M1/251‑118 (MQ=255)                                                                                                                     cccGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatata                                                                                                                                 <  1:585327‑M1/251‑118 (MQ=255)                                                                                                                                ggAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgag                                                                                                                      <  1:457584‑M1/251‑129 (MQ=255)                                                                                                                                   atatGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggt                                                                                                                   <  1:242868‑M1/251‑132 (MQ=255)                                                                                                                                   atatGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagaaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggt                                                                                                                   <  1:612853‑M1/251‑132 (MQ=255)                                                                                                                                           gCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaac                                                                                                           <  1:38555‑M1/251‑140 (MQ=255)                                                                                                                                             cAGACGTTCGACTTCTGCGACGGTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacag                                                                                                         <  2:177571‑M1/251‑142 (MQ=255)                                                                                                                                                 cGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagt                                                                                                     <  1:565462‑M1/251‑146 (MQ=255)                                                                                                                                                  gTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagt                                                                                                     >  1:740402‑M1/1‑105 (MQ=255)                                                                                                                                                     cGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTGCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtattt                                                                                                 <  2:96955‑M1/251‑150 (MQ=255)                                                                                                                                                                  gTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcg                                                                                    >  1:272794‑M1/1‑89 (MQ=255)                                                                                                                                                                    tAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgc                                                                                   <  1:468858‑M1/250‑164 (MQ=255)                                                                                                                                                                    tAGCGGTGTCCAGATACAGTTCCAGGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagtcagaaatatattgcaaacgagggtgttagaacagaagtatttcagaaaaccctcgcg                                                                                  >  2:718966‑M1/1‑87 (MQ=255)                                                                                                                                                                         gtgtCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaa                                                                              >  1:94305‑M1/1‑82 (MQ=255)                                                                                                                                                                              cAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtgagaacctgtcttgaaaaaatatcgtcgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacg                                                                        <  1:617044‑M1/251‑175 (MQ=255)                                                                                                                                                                                   aCAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgtgagaacagaagtatttcagaaaaccctcgcgcaaaagcacgaggg                                                                    <  2:730571‑M1/250‑179 (MQ=255)                                                                                                                                                                                         ccATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgca                                                              <  1:609780‑M1/250‑185 (MQ=255)                                                                                                                                                                                             gatgTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaaga                                                         <  1:226392‑M1/251‑190 (MQ=255)                                                                                                                                                                                                  ttCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaactagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaa                                                     <  1:387702‑M1/250‑194 (MQ=255)                                                                                                                                                                                                    ccTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCCGAGGCtgtgaggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtctattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaaaagaagtttttcagaaaaccctcgcgnaaaagcacgagggtttgcagaagagga                                                      >  2:379692‑M1/1‑55 (MQ=255)                                                                                                                                                                                                      ttttACAGTTCCAGTTAATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGTGGCtgtgtggaatgcataggaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaagaatatctgaaggtatcttcagaatttgaataaactgcagacagcactatattgaaaacgagggcgttagaacagaagtatttcagaaaaccctcgcgcanaagccagagggtttgcagaagatgaagattc                                                >  2:394592‑M1/1‑53 (MQ=255)                                                                                                                                                                                                          aCAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagcc                                             <  1:306230‑M1/250‑202 (MQ=255)                                                                                                                                                                                                                 cagttaATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgacaaaatatcgccgaatgtaaagacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgcaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgaaaaagcacgggggtttgcagaagaggaagattagccgggattac                                     >  2:194528‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                  agttcATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacg                                    >  1:751858‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                   gttcATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgc                                   <  1:485796‑M1/251‑212 (MQ=255)                                                                                                                                                                                                                     tcATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcat                                 >  1:588661‑M1/1‑38 (MQ=255)                                                                                                                                                                                                                      cATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcat                                 <  2:7211‑M1/250‑214 (MQ=255)                                                                                                                                                                                                                       aTGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatac                               >  1:758192‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                         gTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaaagacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaaagagggagttagaacagaagtatttcagaaacccctcgcgcaaaagcacgcgggtttgcagaagaggaagattagccggcataccgcataccg                             >  2:398615‑M1/1‑34 (MQ=255)                                                                                                                                                                                                                              agcagGC‑TGGCGATAGCGTCGGCAGAGGCtgtgtggaacgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgatagaacagaagtatttcagaaaaccctcgagcaaaagcacgagggtttgcagaagaggaagattacccggcattacgcatacctgccgc                        >  1:60703‑M1/1‑29 (MQ=255)                                                                                                                                                                                                                                      tGGCGATAGCGTCGGCAGAGGCtgagtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgactgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccg                 >  2:520566‑M1/1‑22 (MQ=255)                                                                                                                                                                                                                                                 gcgGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatcggaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgac      <  1:94052‑M1/250‑241 (MQ=255)                                                                                                                                                                                                                                                  cGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgacc     >  1:183434‑M1/1‑10 (MQ=255)                                                                                                                                                                                                                                                        agGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgaccatt  >  2:152266‑M1/1‑4 (MQ=255)                                                                                                                                                                                                                                                        agGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaattggaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaacccctcgcgaaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgaccatt  <  1:152266‑M1/248‑245 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                     CGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGC‑TGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACG  >  NC_000913/4139521‑4140014

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