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	59 [0]	[0] 57	[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)	57 (0.790) +TG	52/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

GCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAAC  >  NC_000913/4139520‑4140013                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    gcgGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAggc                                                                                                                                                                                                                                                     >  1:561437/1‑251 (MQ=255)   ggCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtg                                                                                                                                                                                                                                                   <  1:543250‑M1/251‑3 (MQ=255)       caggggctgagcaAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCTCTTCCCAGATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCACGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCAGGATGTTTCCGTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtg                                                                                                                                                                                                                                                 <  1:274078‑M1/243‑5 (MQ=255)           gTCTGAGCAAACAGAATGCCCACATCACCAATCGCTTTTTGCAGGAGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGAGATAATGCTCGGGTTAGTTGTCACCCCGGCAAGGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTGTACAGTTCCAGTTCATGGTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaat                                                                                                                                                                                                                                           <  2:137727‑M1/251‑11 (MQ=255)           gTCTGAGCAAAAAGAACGCCCTCATCACAAATCTCTTTTTGCATACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTATCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGAtgtgtggaat                                                                                                                                                                                                                                           <  2:515941‑M1/251‑11 (MQ=255)               gAGCAAACAGAATGCCCTCATCCCCAATCTCCTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCGGGCAGCGATAATTCTCGGGTTAGTGGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGGTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgc                                                                                                                                                                                                                                         <  2:32665‑M1/249‑13 (MQ=255)               gAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcat                                                                                                                                                                                                                                       <  1:370292‑M1/251‑15 (MQ=255)                aGCAAACAGAATGCCCTCATCACCAATCGCTGTTTGCAGACGCGGCAGCACTTCTCATATGGACGCCTTGCTGGCAGCGATAATGCTCGGTGTAGTTGTCACCCCTGCAATGGGGAATATTCGTGCCAGACGTTCGACTTCTGCGACTTTAGCGGTGTGCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCATTTCATGTTGCAGCAGACTGGCGATAGCGTCGGCAGAGGCtgtgtgtattgcata                                                                                                                                                                                                                                      <  2:499121‑M1/251‑16 (MQ=255)                aGCAAACAGAATGACCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCCCCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGGGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCAAGTTCATGTTGCAGCGGGCTGGCGNTAGCGTCGGAAGAGGCtgtgtggaatgcatc                                                                                                                                                                                                                                      >  2:287417‑M1/1‑236 (MQ=255)                  cAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACTCGGCAGCCCTTCCCACATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCTTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGGCAAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagt                                                                                                                                                                                                                                    <  1:297097‑M1/251‑18 (MQ=255)                        gAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCTGCAGAGGCtgtgtggaatgcatagtaagaac                                                                                                                                                                                                                              <  1:413526‑M1/251‑24 (MQ=255)                         aaTGCCCTCATCCCCAATCGCTTTTTGCAGACGCGGCAGCTCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCCCCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTGCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacc                                                                                                                                                                                                                             <  2:105534‑M1/251‑25 (MQ=255)                              ccTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtct                                                                                                                                                                                                                        <  1:141505‑M1/251‑30 (MQ=255)                                  atcaCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGATAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGGGTCCAGATACAGTCCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCCGAGGCtgtgtgggatgcatagtaagaacctgtcttga                                                                                                                                                                                                                     >  2:564174‑M1/1‑218 (MQ=255)                                    caCCAATAGCTTTTTGACGACGCGGCAGCCCATCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGCTACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGAAGCAGGCTGGCGATCGCGTCGGCAGAGGCtgtggggaatacatagtcagaacctgtcttgaaa                                                                                                                                                                                                                   >  2:43809‑M1/1‑216 (MQ=255)                                     aCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtctttaaaaa                                                                                                                                                                                                                 <  1:125596‑M1/251‑37 (MQ=255)                                      ccAATCGCTTGTTGCAGACGCGGCAGCACTTCCCATATGGACTACTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCCCCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaa                                                                                                                                                                                                                <  2:98484‑M1/251‑38 (MQ=255)                                               ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCGCCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgcc                                                                                                                                                                                                         <  1:294487‑M1/249‑45 (MQ=255)                                                      gACGCGGCAGCACTTCCCATATGGACTCCTGGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatagcgccgaatgtaac                                                                                                                                                                                                <  2:98375‑M1/251‑54 (MQ=255)                                                       acggggcagcACTTCCCATAGGGACTCCTTGCTGGCAGAGATAATGCTCGGGTTAGTTGTCACCCCGGCAAGGGGGAATATGCGTGCAAGACGGTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCAGGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacg                                                                                                                                                                                               <  2:319698‑M1/247‑55 (MQ=255)                                                              agcaCTTCCCATATTGACTCCTTGCTTGCAGCGATCATGCTCGGGTTAGTTTTCACCCCGGCACTGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGCTACAGTTCCCTGATGTTTCCTTTTACAGTTCCAGATCATGTTGCAGCAGGCTGGCGATACCGTCGGCAGAGGCtgtgtggaatgcatggtaagaacatgtcttgaaacaatatcgccgcatgaaacgacaattc                                                                                                                                                                                        >  2:221908‑M1/1‑190 (MQ=255)                                                                 aCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattcctt                                                                                                                                                                                     >  1:140209‑M1/1‑187 (MQ=255)                                                                     cccATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaagg                                                                                                                                                                                 <  1:504266‑M1/251‑69 (MQ=255)                                                                              aCTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaa                                                                                                                                                                        <  1:205580‑M1/251‑78 (MQ=255)                                                                                 ccTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacg                                                                                                                                                                                               >  1:503177‑M1/1‑171 (MQ=255)                                                                                 ccTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacg                                                                                                                                                                                               <  2:503177‑M1/225‑55 (MQ=255)                                                                                     gCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatat                                                                                                                                                                 >  1:81450‑M1/1‑167 (MQ=255)                                                                                       tGGCAGCGATAATGCTCGCGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattc                                                                                                                                                               >  2:480008‑M1/1‑165 (MQ=255)                                                                                            gCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagact                                                                                                                                                          >  1:181910‑M1/1‑160 (MQ=255)                                                                                                 aaTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGCTAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttga                                                                                                                                                      >  1:178287‑M1/1‑155 (MQ=255)                                                                                                       cGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaat                                                                                                                                               >  1:8162‑M1/1‑149 (MQ=255)                                                                                                       cGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaa                                                                                                                                                <  1:526114‑M1/250‑102 (MQ=255)                                                                                                            tAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagac                                                                                                                                          >  1:547322‑M1/1‑144 (MQ=255)                                                                                                             aGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTACTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggaatattcagaatttgaataaaatgcagac                                                                                                                                         >  2:511379‑M1/1‑143 (MQ=255)                                                                                                                        cGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatatt                                                                                                                              >  1:491544‑M1/1‑132 (MQ=255)                                                                                                                                        gCGTGCCATGCGTTCGACTTCTGCGACTTTAGCTGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatctccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgacaacgagggtgtta                                                                                                              <  2:520420‑M1/251‑136 (MQ=255)                                                                                                                                         cGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttag                                                                                                             <  1:512490‑M1/251‑137 (MQ=255)                                                                                                                                               aGACGTGCGACTTCTGCTACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGTAGAGGCtgtgtggaatgcatagtaagaacctttcttgaaaaaatatcgccgaatgtaacgaccattccttaaggagatctgaaggtagattcagaatttgaataaaatgcagacagaaatatattgaaaacgcggttgttacaacaga                                                                                                       <  2:554882‑M1/251‑143 (MQ=255)                                                                                                                                                gACGTTCGACTTCTGCGACGTTAGCGTTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaa                                                                                                      <  1:252494‑M1/251‑144 (MQ=255)                                                                                                                                                              gCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTGGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaa                                                                                         <  1:374868‑M1/250‑157 (MQ=255)                                                                                                                                                                gACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccc                                                                                      >  2:105698‑M1/1‑92 (MQ=255)                                                                                                                                                                gACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatggaacgacaattccttaaggatatctgaaggtaaattcagaatttgaataaaatgcagacagaaatatattgaaaacganggtgttagaacagaaggatttcagaaaaccc                                                                                      >  2:292550‑M1/1‑92 (MQ=255)                                                                                                                                                                 aCGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGATCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaactgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccct                                                                                     >  1:491516‑M1/1‑91 (MQ=255)                                                                                                                                                                  cGTTAGCGGTGTCCAGATTCAGTTCCATGATGTTTCCTTTTACATTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaacccgc                                                                                    <  2:82527‑M1/251‑162 (MQ=255)                                                                                                                                                                       gCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgca                                                                               <  2:140209‑M1/251‑167 (MQ=255)                                                                                                                                                                              ccAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcac                                                                        >  1:214450‑M1/1‑78 (MQ=255)                                                                                                                                                                                            atgatgTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaa                                                          <  1:549596‑M1/251‑188 (MQ=255)                                                                                                                                                                                            atgatgTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGAGAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaa                                                          <  1:43809‑M1/251‑188 (MQ=255)                                                                                                                                                                                             tgatgTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcaga                                                           >  2:374279‑M1/1‑63 (MQ=255)                                                                                                                                                                                                 gTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttgaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagaacgagggtttgcagaagagg                                                      >  2:6992‑M1/1‑59 (MQ=255)                                                                                                                                                                                                        tttACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagatt                                                <  1:389166‑M1/249‑198 (MQ=255)                                                                                                                                                                                                          tgaagtttcagttcatgTGGCAGCAGGGTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattacttaaggatatatgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgatggtgttagaacagaagtagtacagaaaaccctcgcgcaaaagcacgatggtttgcagaagaggaagattagc                                             <  2:127512‑M1/247‑201 (MQ=255)                                                                                                                                                                                                                      tcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgca                                 <  1:287417‑M1/250‑213 (MQ=255)                                                                                                                                                                                                                       cATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcata                               <  2:432719‑M1/251‑215 (MQ=255)                                                                                                                                                                                                                                    gCTGGCGATAGCGTCGGCAGAGGCtgtggggaatgcatagtaagaacctgtcttgaaaaaatatcgccgcatgaaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacaagggtgttagaacagaaggatttcagaaaaccctctcgcaaacgcacgagggtttgcagaagaggaagataagccggtattacgcatacctaccgcaaacc                  >  2:581382‑M1/1‑24 (MQ=255)                                                                                                                                                                                                                                        gCGATAGCGTCGGCAGAGGCtgttgggaatgcagagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgtgggtgttagaacagaagtatttcaaaaacccctcgcgcaaaagcacgagggtttgcagaagaggaaaattagccggtattacgcatacctgccgcaatcccggc              <  2:514425‑M1/251‑232 (MQ=255)                                                                                                                                                                                                                                           aTAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaat           <  1:467108‑M1/251‑235 (MQ=255)                                                                                                                                                                                                                                            aagagtcggCAGAGGCagtggggaatgaatagtaagaacctggcgagaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtaaattcagaattggaagaaaatgcaaacagaaatatattgaaaacgagggtgttagaacagaagtatatcagaaaaccctcgcggaaaagcacgagggttggcagaagaggaagagtagccggtattacgcatccctgccgcaatcccggcaata          <  2:232378‑M1/247‑236 (MQ=255)                                                                                                                                                                                                                                                    gCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgaaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagcagaggaagattagccggtattacgcatacctgacgcaatcccggcaatagtgaccat  >  1:238249‑M1/1‑8 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                    GCGGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAAC  >  NC_000913/4139520‑4140013

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