Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,582,442	Δ986 bp	IS1‑mediated	[yhhZ]–[yrhA]	[yhhZ], [yrhA]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3582442	3583983–3583427	986–1542	44 [1]	[11] 12	[yhhZ]–[insB6]	[yhhZ],yrhA,insA6,[insB6]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3582441	1 (0.020)	41 (0.630)	39/478	0.6	97.6%	coding (579/1179 nt)	yhhZ	putative endonuclease YhhZ
	?	NC_000913	3583428 =	NA (NA)					noncoding (1/768 nt)	IS1	repeat region

TTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTGAAGGATTCAGAGATACTGCGGAAAAACTGGGTGGAGATGCAATAAAAGGCTATGATCAAATACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGACATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCT  >  NC_000913/3582193‑3582690                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          ttATCGGCTATTCATCAACAGATCATTGACAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTgg                                                                                                                                                                                                                                                         >  1:496557‑M1/1‑249 (MQ=255)  tATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggg                                                                                                                                                                                                                                                        >  1:410355‑M1/1‑248 (MQ=255)                                              aaaacAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgag                                                                                                                                                                                                            <  2:650555‑M1/250‑48 (MQ=255)                                                 aCAATAACTATTAGTTATGAATTTATCCCCTGTCAACATCTTAGCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtg                                                                                                                                                                                                         <  2:410355‑M1/251‑51 (MQ=255)                                                             aGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggctt                                                                                                                                                                                             >  1:265553‑M1/1‑189 (MQ=255)                                                              gTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttc                                                                                                                                                                                            >  2:702138‑M1/1‑188 (MQ=255)                                                                   gAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCACATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgttt                                                                                                                                                                                       >  1:54462‑M1/1‑183 (MQ=255)                                                                                gTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtcc                                                                                                                                                                          >  1:702185‑M1/1‑170 (MQ=255)                                                                                           aTCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctttcagctgtccctcctgttc                                                                                                                                                                 >  1:286937‑M1/1‑159 (MQ=255)                                                                                                    aCCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgac                                                                                                                                                       <  1:255129‑M1/250‑101 (MQ=255)                                                                                                        aGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACAtt                                                                                                                                                                                                                                                           <  2:79197/146‑1 (MQ=255)                                                                                                        aGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACAtt                                                                                                                                                                                                                                                           >  1:79197/1‑146 (MQ=255)                                                                                                                 agttggCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaa                                                                                                                                         <  1:247148‑M1/249‑115 (MQ=255)                                                                                                                            ccTGAAAATTATAACCGTTAGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGATACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtggatgatggggtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcc                                                                                                                              >  2:415393‑M1/1‑126 (MQ=255)                                                                                                                             cTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcac                                                                                                                             >  1:41861‑M1/1‑125 (MQ=255)                                                                                                                                           cGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagc                                                                                                               >  1:533352‑M1/1‑111 (MQ=255)                                                                                                                                               tgtttTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgagtaacggcaaaagcaccgccggacatcagcgct                                                                                                            >  1:587399‑M1/1‑107 (MQ=255)                                                                                                                                                 tttttACCAAATTCAAAAAACCACACAAATAATCGTTTCAAAACTTTAAACAGCAAAGCTATTAGCAGGCTACTTGCTGCTGGTGGCGTATCCAAAGGGCACATTtgtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagatactgacggggtgttgcgtaacggcaaaaaccccgccggacatcagcgctatc                                                                                                         >  2:551134‑M1/1‑105 (MQ=255)                                                                                                                                                                     ccAAACAAATAATCGTTTCAAAACGTTCAACAGCAACGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtca                                                                                      >  2:262659‑M1/1‑85 (MQ=255)                                                                                                                                                                       aaacaaaTAATCTTGTAAAAACGTTAAACAGAAAAGCTATTGGCAGGCTTCTTGCTGCTGTTGGCGTATACAATGGGAACATGggtgatgctgccaacttactgatttagtgtatgattgtgttgttgaggtgctccagtggcttctgttgctaccagctgtcccccctgtgcagctactgacggggtggtgcgtaacgccaaaagcaccgccggacatcagcgctatctcttctctcactgccgtaaaac                                                                                   <  2:55820‑M1/251‑169 (MQ=255)                                                                                                                                                                         acaaaTAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTGAAGGATTCAGTGATACTGCGGAAAAACTGGGTGGAGATGCAATAAAAGGCTATGATCAAATACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATCCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACt                                                                                 >  1:219736/1‑251 (MQ=255)                                                                                                                                                                              taatCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcac                                                                            >  1:88513‑M1/1‑76 (MQ=255)                                                                                                                                                                                       cAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttc                                                                   >  1:439917‑M1/1‑67 (MQ=255)                                                                                                                                                                                        aaaaCGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactggtttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggccaccgcagttc                                                                   >  2:63078‑M1/1‑66 (MQ=255)                                                                                                                                                                                         aaaCGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcac                                                                 >  1:73970‑M1/1‑65 (MQ=255)                                                                                                                                                                                          aaCGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttcctgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgccggggtggggcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcact                                                                >  2:249719‑M1/1‑64 (MQ=255)                                                                                                                                                                                                   cAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttaacttacaccgct                                                       >  1:304122‑M1/1‑55 (MQ=255)                                                                                                                                                                                                      cAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttct                                                    >  1:83180‑M1/1‑52 (MQ=255)                                                                                                                                                                                                      cAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacgggggggtgcgtaacggcaaaagcaccgccgaacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccacttct                                                    >  2:555935‑M1/1‑52 (MQ=255)                                                                                                                                                                                                            tATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaaccc                                              >  1:131423‑M1/1‑46 (MQ=255)                                                                                                                                                                                                                ggCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggagcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcacctgcagttcacttacaccgcttctcaccccggt                                           >  1:594244‑M1/1‑42 (MQ=255)                                                                                                                                                                                                                      cTACTTGCTGCTGGTGGCGTATACAATGGGACCATTggtgatgctgccaacttcctgctttcgtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcacc                                    >  1:446027‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                      cTACTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcacc                                    >  1:392390‑M1/1‑36 (MQ=255)                                                                                                                                                                                                                        gcgtgcggctggtggcgtATACAATGGGAACATTggtgatgctgccaacttactgatttagtggatgatggtgtttttgaggtgctccagtggcttgtgtttctatcagctgcccctcctgttcagctactgacggggtggggcgtaacggcaaaagcaccgccggacatcagcgctatctcggctctcactgccgtaaaacatggcaactgcagtgcacttacaccgcttctcatcccggtacgcaccag                                  <  2:439917‑M1/244‑218 (MQ=255)                                                                                                                                                                                                                         cTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggac                                                                                                                     <  1:483360‑M1/167‑135 (MQ=255)                                                                                                                                                                                                                         cTTGCTGCTGGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggac                                                                                                                     >  2:483360‑M1/1‑33 (MQ=255)                                                                                                                                                                                                                            gctgctGGTGGCGTATACAATGTGATCATGggggatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaa                              <  1:551134‑M1/251‑222 (MQ=255)                                                                                                                                                                                                                             ctgctgGTGGCGTATACAATGGGAACATTggagatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgccggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaac                                                >  1:400595‑M1/1‑29 (MQ=255)                                                                                                                                                                                                                             cggctggtGGCGTATACAAGGGGACCATTgtggatgcggccaacttactggttttgggtatgatggtgtttttgaggggccccagtggcttcttgttctttcagctgtccctcctggtcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcatcgctatctctgctctcactgccgtaaaacatggcaactgcagttcactgacaccgcttctcaac                                                <  2:400595‑M1/230‑204 (MQ=255)                                                                                                                                                                                                                              tgctgGTGGCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaccatggcaactgcagttcacttccaccgcttctcaacccggtacgcaccagaaaatc                            >  2:221983‑M1/1‑28 (MQ=255)                                                                                                                                                                                                                                     ggCGTATACAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcattgat                      <  1:702138‑M1/250‑230 (MQ=255)                                                                                                                                                                                                                                           taCAATGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctcccgtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaacagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatgtcaactgcagttcacttacaccgcttctcaacccggtccgccacagaaaatcattgatatggcc                >  2:595280‑M1/1‑15 (MQ=255)                                                                                                                                                                                                                                               aTGGGAACATTggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcattgatatggccatga            >  1:79479‑M1/1‑11 (MQ=255)                                                                                                                                                                                                                                                         tggtgatgctgccaacttactgaattagtgtatgatggtgtttttgaggtgctccagtggcttctgtttatatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctatcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcattaatatggccatacatggcgttgg  >  2:643349‑M1/1‑1 (MQ=255)                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                          TTATCGGCTATTCATCACCAGATCATTGAAAACCAACTGGATACAGAAACAATAACTATTAGTTATGAATTTATCCTCTGTCAACATCTTATCGCAAATACCGAGTTCAGCTATTTGGCACTCCCTGAAAATTATAACCGTTTGTTTTTACCAAATTCAAAAAACCAAACAAATAATCGTTTCAAAACGTTAAACAGCAAAGCTATTGGCAGGCTACTTGCTGCTGGTGGCGTATACAATGGGAACATTGAAGGATTCAGAGATACTGCGGAAAAACTGGGTGGAGATGCAATAAAAGGCTATGATCAAATACTAAATGAAAAAACAGCGGGCATAGCGATAGCAACAGCATCTATTCTTTTAACAAAGCGTTCTAATGTTGATACATATACAGAAATAAATAGTTACTTAGGCAAACTTAGAGGTCAACAAAAACTTCTTGATGGTATAGACATAATAGAAATAATATACATTAAGAGACCTTCAAAAGACTTAGCT  >  NC_000913/3582193‑3582690

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