Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	4,620,828	43 bp→AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT	intergenic (+2/‑13)	deoB → / → deoD	phosphopentomutase/purine nucleoside phosphorylase

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	4620828	4620870	43	42 [0]	[0] 65	deoB/deoD	phosphopentomutase/purine nucleoside phosphorylase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 4620827	0 (0.000)	42 (0.490) +AATGGGTGCGGCTTGATGCCCTCACCCCGGCCCTCTCCCACAGGGAGAGGGTGAAAACATGTTGAT	22/466	2.9	100%	intergenic (+1/‑56)	deoB/deoD	phosphopentomutase/purine nucleoside phosphorylase
	?	NC_000913	4620871 =	0 (0.000)					intergenic (+45/‑12)	deoB/deoD	phosphopentomutase/purine nucleoside phosphorylase

AAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCG  >  NC_000913/4620722‑4621126                                                                                                          |                                                                                                                                                                                                                                                                                                            cactatcaataagttggagtcattaccaccttcgcggatattggtcagactctggCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtga                                                                                                           <  1:370550‑M1/273‑195 (MQ=255)   ctatcaataagttggagtcattaccaccttcgcggatattggtcagactctGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaag                                                                                                         >  1:193386‑M1/26‑104 (MQ=255)   ctatcaataagttggagtcattaccaccttcgcggatattggtcagactctGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaag                                                                                                         <  2:131220‑M1/275‑197 (MQ=255)     atcaataagttggagtcattaccaccttcgcggatattggtcagactCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtg                                                                                                       >  2:141658‑M1/24‑102 (MQ=255)     atcaataagttggagtcatcaccaccctcgcggatattggtcagactctggcaaaATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgtggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgcaaggcgaccagctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtg                                                                                                       >  2:164240‑M1/28‑102 (MQ=255)             gttggagtcattaccaccttcgcggatattgGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgt                                                                                               <  1:283754‑M1/285‑207 (MQ=255)                   gtcattaccaccttcgcggATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcgg                                                                                         >  1:2099‑M1/10‑88 (MQ=255)                    tcattaccaccttcgcgGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggt                                                                                        <  1:444247‑M1/292‑214 (MQ=255)                     cattaccaccttcgcGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGNTATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagnatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggta                                                                                       >  2:801920‑M1/8‑86 (MQ=255)                     cattaccaccttcgcGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcaganatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggta                                                                                       >  1:597834‑M1/8‑86 (MQ=255)                     cattaccaccttcgcGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagnatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggta                                                                                       >  2:648509‑M1/8‑86 (MQ=255)                     cattaccaccttcgcGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcacccaggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgacgtatattgctgaaactttccttgaagatgcccgtgaagtgaacaaagttcgcggta                                                                                       >  2:538891‑M1/8‑86 (MQ=255)                        taccacnttCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGANATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgc                                                                                    <  2:726287‑M1/296‑218 (MQ=255)                             ccTTCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacaggtagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggc                                                                               <  1:756537‑M1/300‑223 (MQ=255)                             ccTTCGAGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggc                                                                               <  1:860194‑M1/300‑223 (MQ=255)                               ttCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggctt                                                                             >  1:866951‑M1/1‑76 (MQ=255)                               ttCGCGGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatattgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacacgctgcgtgcgaagtatattgctgaaactttccttgaagatgcacgtgaagtgaacaacgttcgcggtatgctgggctt                                                                             >  1:470577‑M1/1‑76 (MQ=255)                                  gcgGATATTGGTCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcac                                                                          <  1:488779‑M1/300‑228 (MQ=255)                                          tggtcagacTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtactt                                                                  >  1:505085‑M1/5‑65 (MQ=255)                                           ggtcagaCTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtactta                                                                 >  1:437271‑M1/4‑64 (MQ=255)                                            gtccgactcTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttnatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttac                                                                >  2:768060‑M1/5‑63 (MQ=255)                                             tcaGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATNTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgnagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttaca                                                               <  2:796569‑M1/299‑239 (MQ=255)                                             tcaGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttaca                                                               <  1:54241‑M1/299‑239 (MQ=255)                                                gACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaag                                                            <  1:248812‑M1/300‑242 (MQ=255)                                                gACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatacccgtgaagcgaacaacgttcgcggtatgttgggcttcaccggtacttacaaag                                                            >  2:276962‑M1/1‑59 (MQ=255)                                                gACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGGTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaag                                                            <  1:782221‑M1/300‑242 (MQ=255)                                                                       cTTCTGATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggt                                     >  1:733528‑M1/1‑36 (MQ=255)                                                                            gATATGGAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacgg                                <  1:63872‑M1/300‑270 (MQ=255)                                                                                  gAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatggg                          >  1:8991‑M1/1‑25 (MQ=255)                                                                                  gAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatggg                          <  2:838215‑M1/300‑276 (MQ=255)                                                                                  gAATATGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctaggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatggg                          >  2:560202‑M1/1‑25 (MQ=255)                                                                                    atatGGCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatgggta                        <  2:89171‑M1/300‑278 (MQ=255)                                                                                        ggCAAAGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggatgaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatgggtatccc                    <  1:207232‑M1/300‑282 (MQ=255)                                                                                            aaGCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatgggtatcccgtcc                <  1:557169‑M1/300‑286 (MQ=255)                                                                                              gCCATGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatgggtatcccgtcctg              >  1:60450‑M1/1‑13 (MQ=255)                                                                                                 aTGTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatgggtatcccgtcctgctc           <  1:573069‑M1/300‑291 (MQ=255)                                                                                                   gTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggaganggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacnaaggccgcaaaatttccgtaatgggtcacggtatgggtatcccgtcctgctcca         >  2:806100‑M1/1‑8 (MQ=255)                                                                                                   gTTCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggatcgaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgaacagctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacnaaggccgaaaaatttccgtaatgggtcacggtatgggtatcccgcccngctcca         >  2:726357‑M1/1‑8 (MQ=255)                                                                                                    ttCTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgctaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcgggatgctgggcttcaccgggactgacaaaggtagcaagatttgcgtaatgggtcgcggggggggggttcgggggggataggt        <  1:232058‑M1/300‑294 (MQ=255)                                                                                                      cTGATaatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatgggtatcccgtcctgctccatct      <  1:836456‑M1/300‑296 (MQ=255)                                                                                                        gATaatgggtgaggcttgaggcccccaccccggccctccaacacagggagagggtgaaaacatgttgataaggataaaacaatgactaccccacacattaatgcagaactgggcgatctcgctgacgtagttttgatgcaaggcgacccgctgcgtgcgaagtatatngctgaaacttcccttgaagatgcacgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttaaaaaggccgcacaatttccgtaatggatcccggtacgggtatcccgtcctgctcactctac    >  2:209454‑M1/1‑3 (MQ=255)                                                                                                          taatgggtgcggcttgatgccctcaccccggccctctcccacagggagagggtgaaaacatgttgataaggataaaacaatggctaccccacacattaatgcagaaatgggcgatttcgctgacgtagttttgatgccaggcgacccgctgcgtgcgaagtatattgctgaaactttccttgaagatgcccgtgaagtgaacaacgttcgcggtatgctgggcttcaccggtacttacaaaggccgcaaaatttccgtaatgggtcacggtatgggtatcccgtcctgctccatctacac  <  1:516138‑M1/300‑300 (MQ=255)                                                                                                          |                                                                                                                                                                                                                                                                                                            AAACCGGGCTCACTGGGTCATCGTGAAACCTTCGCGGATATCGGCCAGACTCTGGCAAAATATTTTGGTACTTCTGATATGGAATATGGCAAAGCCATGTTCTGATGGATTTGGGCGGAGCGTTGACTCCGCCTTTGTTATGTCACAAAAAGGATAAAACAATGGCTACCCCACACATTAATGCAGAAATGGGCGATTTCGCTGACGTAGTTTTGATGCCAGGCGACCCGCTGCGTGCGAAGTATATTGCTGAAACTTTCCTTGAAGATGCCCGTGAAGTGAACAACGTTCGCGGTATGCTGGGCTTCACCGGTACTTACAAAGGCCGCAAAATTTCCGTAATGGGTCACGGTATGGGTATCCCGTCCTGCTCCATCTACACCAAAGAACTGATCACCGATTTCG  >  NC_000913/4620722‑4621126

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