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	66 [0]	[0] 67	[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)	65 (0.690) +TG	60/484	0.4	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

GGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACG  >  NC_000913/4139522‑4140014                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     ggCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtg                                                                                                                                                                                                                                                    >  2:838532‑M1/1‑249 (MQ=255)  gCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgt                                                                                                                                                                                                                                                   >  2:671187‑M1/1‑248 (MQ=255)     cATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGCGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtg                                                                                                                                                                                                                                                <  1:547759‑M1/251‑7 (MQ=255)         gTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaat                                                                                                                                                                                                                                            <  2:533160‑M1/251‑11 (MQ=255)             gAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcat                                                                                                                                                                                                                                        <  2:445538‑M1/251‑15 (MQ=255)             gAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgca                                                                                                                                                                                                                                         <  1:140109‑M1/250‑14 (MQ=255)             gAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATACGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTATTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcat                                                                                                                                                                                                                                        <  2:661597‑M1/251‑15 (MQ=255)                cAAACAGCATGCCGTCATCACCAATCGCTTTTTGCAGACGCGTCAGCACTTCCCATATGTACTCCTTGCGGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCAATAGAGTCGGCAGAGTCtgtgtggaatgcatagt                                                                                                                                                                                                                                     <  2:45996‑M1/251‑18 (MQ=255)                 aaaCAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGCTACAGTTCCATGATGTTTCCTTTTACAGTTCCCGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgcgtggaatgcatagtt                                                                                                                                                                                                                                    >  2:306331‑M1/1‑233 (MQ=255)                     aGAATGCCCGCATCACCAATCGCTTTTTGCAGACGCGGCAGAACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGGTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCGTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaa                                                                                                                                                                                                                                <  1:727057‑M1/251‑23 (MQ=255)                      gAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaac                                                                                                                                                                                                                               <  1:415041‑M1/251‑24 (MQ=255)                               catcaCCACTCGCTTTTTGCAGACGCGGCAGCACTTCCCATAGGGACTCCTGGCTGGCAGAGATAATGCTCGGGTGAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttg                                                                                                                                                                                                                       <  2:20630‑M1/250‑32 (MQ=255)                                  caCCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaa                                                                                                                                                                                                                    >  2:70554‑M1/1‑216 (MQ=255)                                     cAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTACCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaccctgtcttgaaaaaat                                                                                                                                                                                                                <  1:288999‑M1/251‑39 (MQ=255)                                          gCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgc                                                                                                                                                                                                           >  1:780230‑M1/1‑208 (MQ=255)                                             ttttGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccg                                                                                                                                                                                                         <  1:838532‑M1/250‑46 (MQ=255)                                                     gcgctgcagcacttccCAGATGCACTCATAGCGGGCCGCGATAATGCTCGGGTTAGTTGGCACCCCGGCAATGGGGAAGATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCAGGATGTTTCCTGTGACAGTTCCAGTTCAGGTTGCAGCTGGCTGGCTATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttcaaaaaatatctccgaatgtaacg                                                                                                                                                                                                <  2:718676‑M1/246‑55 (MQ=255)                                                            agcaCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCTCCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattc                                                                                                                                                                                         <  1:297147‑M1/251‑62 (MQ=255)                                                               aCTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattcctt                                                                                                                                                                                      <  1:450862‑M1/251‑65 (MQ=255)                                                                    ccATATGGACTCCTTTCTGGCAGCGATAATGCTCGGGTTAGTTGTCCACCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTTCAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaagga                                                                                                                                                                                 <  2:586946‑M1/251‑70 (MQ=255)                                                                        atGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatc                                                                                                                                                                             >  2:76346‑M1/1‑178 (MQ=255)                                                                           gACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctg                                                                                                                                                                           >  1:671031‑M1/1‑175 (MQ=255)                                                                           gACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctg                                                                                                                                                                           >  1:301132‑M1/1‑175 (MQ=255)                                                                                   gCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatat                                                                                                                                                                  <  1:268109‑M1/251‑85 (MQ=255)                                                                                            gATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatagccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaattt                                                                                                                                                         >  1:90018‑M1/1‑158 (MQ=255)                                                                                             aTAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttg                                                                                                                                                        <  2:401509‑M1/251‑95 (MQ=255)                                                                                               aaTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGGTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttga                                                                                                                                                       <  2:701746‑M1/250‑96 (MQ=255)                                                                                                aTGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaat                                                                                                                                                     <  1:492355‑M1/251‑98 (MQ=255)                                                                                                        gTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgc                                                                                                                                              >  2:742157‑M1/1‑146 (MQ=255)                                                                                                           agatgtcACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacacttccttaaggatatctgaaggtatattcagaatttgaataaaatgcagac                                                                                                                                          >  2:793206‑M1/4‑143 (MQ=255)                                                                                                           aGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagac                                                                                                                                          >  1:246515‑M1/1‑143 (MQ=255)                                                                                                                 cACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaat                                                                                                                                    >  1:119173‑M1/1‑137 (MQ=255)                                                                                                                 cACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaat                                                                                                                                    <  2:777439‑M1/251‑115 (MQ=255)                                                                                                                      cGGCAATGGGGAAGATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatatt                                                                                                                               <  2:203828‑M1/251‑120 (MQ=255)                                                                                                                       ggCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattg                                                                                                                              >  1:166461‑M1/1‑131 (MQ=255)                                                                                                                             ggggAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacagttccttagggatatctgaaggtatattcagaatttgaataaaatgcagacaggactatattgaaaacg                                                                                                                        >  2:573672‑M1/1‑125 (MQ=255)                                                                                                                               ggAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgag                                                                                                                      <  2:90985‑M1/251‑129 (MQ=255)                                                                                                                                gAATATGCGTGCAAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaagatatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagg                                                                                                                     >  2:817571‑M1/1‑122 (MQ=255)                                                                                                                                  atatGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggt                                                                                                                   >  1:204810‑M1/1‑120 (MQ=255)                                                                                                                                    atGCGGGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaccgagggtg                                                                                                                  <  1:257013‑M1/250‑133 (MQ=255)                                                                                                                                           ccAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaaca                                                                                                          <  2:521516‑M1/251‑141 (MQ=255)                                                                                                                                               aCGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaag                                                                                                      <  2:671820‑M1/251‑145 (MQ=255)                                                                                                                                                    cGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtattt                                                                                                 >  1:798675‑M1/1‑102 (MQ=255)                                                                                                                                                      aCTTCTGCGACGTTAGCGGTGTCCAGATACAGTTACATGATGTTTCCTTTTACAGTTCAAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttc                                                                                                >  1:230037‑M1/1‑100 (MQ=255)                                                                                                                                                           tGCGACGTTGGCGGTGTCCAGATCCAGTTCCATGATGTTTCCTTTTTCAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcataggaagatcctgtcttgaaaaaagatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtctttcgtaaa                                                                                           <  2:250575‑M1/250‑156 (MQ=255)                                                                                                                                                             cGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatatggaaagcgagggtgttagaacagaagtatttcagaaaacc                                                                                        <  2:613203‑M1/251‑159 (MQ=255)                                                                                                                                                               aCGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTGCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGAtgcgtggaatgcatagaaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggctatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccct                                                                                      >  1:708324‑M1/1‑91 (MQ=255)                                                                                                                                                                 gTTAGCGGTGTCCAGATACAGTTCAATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGAAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcg                                                                                    >  1:252204‑M1/1‑89 (MQ=255)                                                                                                                                                                   tatcggtGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgactaaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaacccacgc                                                                                   <  2:96984‑M1/247‑164 (MQ=255)                                                                                                                                                                            ccAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcac                                                                         <  2:163670‑M1/251‑174 (MQ=255)                                                                                                                                                                                     gTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattga                                                                                                                             >  2:846213‑M1/1‑69 (MQ=255)                                                                                                                                                                                     gTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattga                                                                                                                             <  1:846213‑M1/190‑122 (MQ=255)                                                                                                                                                                                          atgatgTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaa                                                           <  2:491171‑M1/251‑188 (MQ=255)                                                                                                                                                                                          atgatgTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcaga                                                            >  1:564837‑M1/1‑64 (MQ=255)                                                                                                                                                                                           tgatgTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcaga                                                            >  1:647846‑M1/1‑63 (MQ=255)                                                                                                                                                                                           ggatgtTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaaggattgcagaaaaccctcgcgcaaaatcactagggtttgcagaag                                                          <  2:302569‑M1/250‑189 (MQ=255)                                                                                                                                                                                                                cagttcATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattac                                     <  2:791662‑M1/251‑210 (MQ=255)                                                                                                                                                                                                                 agttcCTGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacg                                    >  1:802951‑M1/1‑41 (MQ=255)                                                                                                                                                                                                                             agcagGGTGGCGATAGCTACGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatatgcagaatttgaataaagtgcagacataaatatatagaaaacgaggttgttataacagaagaatttcagaaacccctcgcacaaaagcacgagggttggcagaagaggaagagtagctgttattacggatccctgcagc                        <  2:34486‑M1/251‑223 (MQ=255)                                                                                                                                                                                                                             agcagGCTGGCGATAGCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgc                        >  2:846836‑M1/1‑29 (MQ=255)                                                                                                                                                                                                                                agGCTGGCGATAGCGTCGGCAGAGGCtgtgtggcatgcatagtaagaacctgctttgaaaaaatatcgcagaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatatttaaaacgagggggttagaacagaagtatttcagaaaaccctcacgcaaaatcacgtgggtttgcagaagaggaagattagccggtattacgcatacctgccgcaat                     >  2:397594‑M1/1‑26 (MQ=255)                                                                                                                                                                                                                                 ggCTGGCGATAGCGTCGGCAGAGGCtgtggggaatgcatagtaagaacctggcttgaaaaaatatcgccgaatgtaacgacaatacattaaggatatctgaaggtatatccagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaacccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggaattacgcatccctgccgcaatc                    <  2:14043‑M1/251‑227 (MQ=255)                                                                                                                                                                                                                                     ggctatagcGGCGGCAGAGGCtgtgtgaaatgcatagtaagaccctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaaagagggtgttagaacagaagtatttcagaaacccctcgcgcaaaagcccgagggtttgcagaagaggaatattagccggtattccccatccctgccgcaatcccgg                <  2:564837‑M1/247‑231 (MQ=255)                                                                                                                                                                                                                                            gCGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagt         <  1:167373‑M1/251‑238 (MQ=255)                                                                                                                                                                                                                                             cGTCGGCAGAGGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaacccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtg        <  2:739118‑M1/251‑239 (MQ=255)                                                                                                                                                                                                                                                    agagGCtgtgtggaatgcatagtaagaacctgtcttgaaaaaatatcgccgaatgtaacgacaattccttaaggatatctgaaggtatattcagaatttgaataaaatgcagacagaaatatattgaaaacgagggtgttagaacagaagtatttcagaaaaccctcgcgcaaaagcacgagggtttgcagaagaggaagattagccggtattacgcatacctgccgcaatcccggcaatagtgaccatt  >  2:425306‑M1/1‑6 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     GGCTCATGGTCTGAGCAAACAGAATGCCCTCATCACCAATCGCTTTTTGCAGACGCGGCAGCACTTCCCATATGGACTCCTTGCTGGCAGCGATAATGCTCGGGTTAGTTGTCACCCCGGCAATGGGGAATATGCGTGCCAGACGTTCGACTTCTGCGACGTTAGCGGTGTCCAGATACAGTTCCATGATGTTTCCTTTTACAGTTCCAGTTCATGTTGCAGCAGGCTGGCGATAGCGTCGGCAGAGGCGGCGTTAACCAGCGCGTTACGGAATTCTTCGTGCATGATGCGGCGAGCGAGGCGCGAGAAAATGCGCATATGCTGATCGCCCGCAGCGTGTTTGTTCAGGGTTAACATAATGATGAATTGCGCTTCATCATCGCCCCAGCGCACCGGCGCTTGCAGACGCGCCACGCTGATGGTGGATTGCTCAATGTGTTCTGATTTGCTGTGTGGAATGGCAAAACTAAAGCCCAGACCGGTAGAGAAAACG  >  NC_000913/4139522‑4140014

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