Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
MC JC	NC_000913	3,708,467	Δ168 bp		[proK]	[proK]

Missing coverage evidence...
			seq id	start	end	size	←reads	reads→	gene	description
*	*	÷	NC_000913	3708467	3708634	168	36 [0]	[28] 64	[proK]	[proK]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3708466	0 (0.000)	36 (0.540)	20/374	1.9	100%	intergenic (‑761/+150)	dppA/proK	dipeptide ABC transporter periplasmic binding protein/tRNA‑Pro
	?	NC_000913	3708635 =	0 (0.000)					noncoding (58/77 nt)	proK	tRNA‑Pro

TTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCT  >  NC_000913/3708175‑3708650                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                         tttAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctc                                                                                                                                                                                  <  1:764958‑M1/300‑9 (MQ=255)                     gACAATGTTTTCTGTCACAAATCCGTTAATAGATAGTGAGATATTGGGCACAAATTTGGCAGTGAAATGTGTGTCGTGTAGTTCGAATAAGGCGTTCGCGCCGCATCCGACAATGGCCTGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTcgaacctccgaccccttcgtcccgaac                                                                                                                                                               <  2:564181‑M1/299‑28 (MQ=255)                       cAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGGTATGGGGCGCAAATTTGGCAGGGAAATGTGTGTAGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCATGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccacttcgtccctaacgaa                                                                                                                                                            >  2:154000‑M1/1‑270 (MQ=255)                           ggTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgc                                                                                                                                                        >  1:13663‑M1/1‑266 (MQ=255)                                     acaAATCCGTTAATGGATAGTTAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGTTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACTTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCACCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                              <  2:255184‑M1/300‑45 (MQ=255)                                     acaAATCCGTTAATGGATAGTTAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                              <  2:505976‑M1/300‑45 (MQ=255)                                     acaAATCCGTTAATGGATAGTTAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                              <  1:680835‑M1/300‑45 (MQ=255)                                     acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                              <  1:787202‑M1/300‑45 (MQ=255)                                     acaAATCCGTTAATGGATAGTGAGATATGGGGCACAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGNTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                              <  1:181988‑M1/300‑45 (MQ=255)                                     acaAATCCGTTAAGGGAGAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                              <  2:571549‑M1/300‑45 (MQ=255)                                            cGTTAATGGATAGTGAGATATGGGGTGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgcca                                                                                                                                       >  2:176922‑M1/1‑249 (MQ=255)                                            cGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgcca                                                                                                                                       >  1:74145‑M1/1‑249 (MQ=255)                                                        gttagatatGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgc                                                                                                                           <  1:572613‑M1/297‑64 (MQ=255)                                                        gTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTTTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgc                                                                                                                           <  1:116415‑M1/300‑64 (MQ=255)                                                        gTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgc                                                                                                                           <  1:116644‑M1/300‑64 (MQ=255)                                                          gagaTATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgc                                                                                                                           >  1:816966‑M1/1‑235 (MQ=255)                                                                  gggCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaaccaccgaatgcggggcgcatc                                                                                                                 <  1:350744‑M1/300‑74 (MQ=255)                                                                           tttGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAATGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgc                                                                                                        <  1:73549‑M1/300‑83 (MQ=255)                                                                           tttGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCTACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccaaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgc                                                                                                        <  2:45280‑M1/300‑83 (MQ=255)                                                                           tttGGCAGGGAAATGTGTGTCGGGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgc                                                                                                        <  1:75114‑M1/300‑83 (MQ=255)                                                                                 aGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagata                                                                                                   <  1:227830‑M1/299‑88 (MQ=255)                                                                                 aGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgacccattcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagata                                                                                                   <  1:811332‑M1/299‑88 (MQ=255)                                                                                          tgtgtCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtc                                                                                         <  1:159954‑M1/300‑98 (MQ=255)                                                                                               cGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcc                                                                                    >  1:828109‑M1/1‑198 (MQ=255)                                                                                                                           gcATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttgtgatc                                                           <  2:168108‑M1/297‑128 (MQ=255)                                                                                                                             aTCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcga                                                      >  2:829545‑M1/1‑168 (MQ=255)                                                                                                                                                       tGATGCTACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagtta                              <  2:420291‑M1/298‑157 (MQ=255)                                                                                                                                                       tGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagtta                              <  2:861688‑M1/298‑157 (MQ=255)                                                                                                                                                          tGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTACCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaataccttttgatcggcgagaaagtcagcaggccgcttagttagccgc                         <  2:189634‑M1/300‑162 (MQ=255)                                                                                                                                                              aCGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcc                     >  2:719556‑M1/1‑135 (MQ=255)                                                                                                                                                                            tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgc         >  2:810234‑M1/1‑121 (MQ=255)                                                                                                                                                                                 ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccnatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaaccccttaatagcaaaatgccttttgatcggcgagaaagtaagcaggccgcttagttagccgctgcctattttgcctgcgggatgt  >  2:291078‑M1/1‑116 (MQ=255)                                                                                                                                                                                 ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcgggatgt  >  1:36429‑M1/1‑116 (MQ=255)                                                                                                                                                                                 ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggcagcttagttagccgctgcctcttttgcctgcgggatgt  >  1:312608‑M1/1‑116 (MQ=255)                                                                                                                                                                                 ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtccagaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcgggatgt  >  1:246615‑M1/1‑116 (MQ=255)                                                                                                                                                                                 ggCCTAAGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGTAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccaaccccttcgtaccgaacgaagtgcgctaacaggctgcgcaaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcgatgagaaagtccgtaggccgcttcgttagccgctgcctcttttgcctacgggatgt  >  2:52978‑M1/1‑116 (MQ=255)                                                                                                                                                                                                                                                                                                    |                                                                                                                                                                                         TTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCT  >  NC_000913/3708175‑3708650

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: ‑