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	47 [0]	[24] 71	[proK]	[proK]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3708466	0 (0.000)	47 (0.830)	31/382	0.5	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

AGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCG  >  NC_000913/3708172‑3708653                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                            aGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaac                                                                                                                                                                                        <  1:35766‑M1/300‑6 (MQ=255) aGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTACAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaac                                                                                                                                                                                        <  1:35121‑M1/300‑6 (MQ=255)        aGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgac                                                                                                                                                                                 <  1:270102‑M1/300‑13 (MQ=255)         gCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgacc                                                                                                                                                                                <  1:267901‑M1/300‑14 (MQ=255)         gCAAATATAAAAAAAGACAACGTTTTCTGTCACAAATCCGTTAAGGGATAGTGAGATATGGGGCACAAATTTGGCAGGGAAATGTGTGTCGTGTAGTTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGTTCGGTGATAGAGGATTcgaacctccgacc                                                                                                                                                                                <  2:298002‑M1/300‑14 (MQ=255)          cAAATATAAAAAAAGACAATGTTTTCTGTCACAAATCCGTGAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTcgaacctccgaccc                                                                                                                                                                               <  1:528665‑M1/300‑15 (MQ=255)             atatAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgacccctt                                                                                                                                                                            >  1:65031‑M1/1‑283 (MQ=255)                  aaaaaaGACACCGTTTTCTGCCACAAATCCGTTAATGGATAGTGAGATATTGGGCGCAAATTTGGCAGGGAAATGTGTTTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcc                                                                                                                                                                       <  2:24254‑M1/300‑23 (MQ=255)                  aaaaaaGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcc                                                                                                                                                                       <  1:445757‑M1/300‑23 (MQ=255)                       aGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaac                                                                                                                                                                  >  1:349814‑M1/1‑273 (MQ=255)                          cAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGACTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCTGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaa                                                                                                                                                               >  2:589002‑M1/1‑270 (MQ=255)                                   cTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGAATcgaacctccgaccccttcgtccagaaaaaagtgcgctac                                                                                                                                                      >  2:600762‑M1/1‑261 (MQ=255)                                        acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGTTCGAATAATGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                 <  2:261995‑M1/300‑45 (MQ=255)                                        acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGTGCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                 <  2:217869‑M1/300‑45 (MQ=255)                                        acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAATGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                 <  1:529155‑M1/300‑45 (MQ=255)                                        acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                 <  1:461495‑M1/300‑45 (MQ=255)                                        acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATNTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                 <  1:472407‑M1/300‑45 (MQ=255)                                        acaAATCCGTTAATGGAGAGTGAGATATGTGGCGCAAATTTGGCATGGAAATGTGTGTCGTGTAGTTCGAATAAGGCGGTCACGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGTAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                 <  2:56747‑M1/300‑45 (MQ=255)                                        acaAATCCGTTAAGGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCTGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCTTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                 <  1:604748‑M1/300‑45 (MQ=255)                                                    aTGGATAGTGAGATATGTTGCGCAAATTTGGCATGTAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACTTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTGGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcac                                                                                                                                     <  2:2629‑M1/300‑57 (MQ=255)                                                      ggATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaacccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccg                                                                                                                                   >  2:327440‑M1/1‑242 (MQ=255)                                                         tAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAANGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaat                                                                                                                                >  1:487746‑M1/1‑239 (MQ=255)                                                          aGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatg                                                                                                                               <  1:667864‑M1/300‑63 (MQ=255)                                                                         gcAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttac                                                                                                                <  2:393049‑M1/300‑78 (MQ=255)                                                                                tGGCAGGTAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCTTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcag                                                                                                         <  1:308238‑M1/300‑85 (MQ=255)                                                                                   cAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagata                                                                                                      <  1:650523‑M1/300‑88 (MQ=255)                                                                                     gggAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcnctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacg                                                                                                    >  2:524533‑M1/1‑211 (MQ=255)                                                                                        aaaTGTGTGTCGTGTAGTTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGTAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccc                                                                                                 <  2:218633‑M1/300‑93 (MQ=255)                                                                                                gtCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaat                                                                                         <  1:249639‑M1/300‑101 (MQ=255)                                                                                                  cGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACTTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcc                                                                                       <  1:63151‑M1/300‑103 (MQ=255)                                                                                                                                       aaTGGCCGTCACATTTGCCTGATGCGACGCCGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGTAGTTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCATGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgacctcttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtc                                                  <  2:145598‑M1/300‑140 (MQ=255)                                                                                                                                       aaTGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgacccnttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaanagcaaaatgccttttgatcggcgagaaagtc                                                  <  2:488488‑M1/300‑140 (MQ=255)                                                                                                                                         ctgcgagatgtgtataagagacagacgctgacgcgtcttatcaggccTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcag                                                >  1:38186‑M1/24‑159 (MQ=255)                                                                                                                                          ggCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaacgaatgcggggcgcatcttactgcgcagatacgccctagtcaatccattaatcgcaaagtgccttttgatcggcgagaaagtcagc                                               >  2:285359‑M1/1‑158 (MQ=255)                                                                                                                                                           gATGCGNCGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCNCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggtgcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagcc                              <  2:579275‑M1/300‑160 (MQ=255)                                                                                                                                                            aTGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccg                             >  1:683115‑M1/1‑140 (MQ=255)                                                                                                                                                                 aCGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtccatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcc                        >  2:681638‑M1/1‑135 (MQ=255)                                                                                                                                                                 aCGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcc                        >  1:521533‑M1/1‑135 (MQ=255)                                                                                                                                                                          cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGACTTTGGTCGGTGATAGAGGATTcgaacttccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcacagatacgccctcgtcaatcccttaatagcaaaatgtcttttgatcggcgagaaagtcagcaggccgattcgttagccgatgcctgttttgcc               >  2:62131‑M1/1‑126 (MQ=255)                                                                                                                                                                          cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaccgaagtgcgctaccaggctgcgccaatcaccgaatgcgnggcgcatcttactgcgcagatacgccctcgtcaatccattaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcc               >  1:500730‑M1/1‑126 (MQ=255)                                                                                                                                                                          cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcc               <  2:695052‑M1/300‑175 (MQ=255)                                                                                                                                                                               tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatccattaatagcaaaatgccttttgatcggcgagaaagtcagcaggcagcttagttagccgctgcctcttttgcctgcgg          >  1:587976‑M1/1‑121 (MQ=255)                                                                                                                                                                                    ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTTGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagaaaaatgtcttctgatcggcgagaaagtaagcaggccgcttagttagcagctgcctctttttcctgcgggatgc     >  2:319861‑M1/1‑116 (MQ=255)                                                                                                                                                                                    ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcgggatgt     >  1:375137‑M1/1‑116 (MQ=255)                                                                                                                                                                                    ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtccagaacgaagtgcgctaacaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcgggatgt     >  2:201858‑M1/1‑116 (MQ=255)                                                                                                                                                                                    ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaacccttcgtcccgaacgaagtgcgctaccaggctgcgccantcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatccattaatagcaaaatgcattttgatcggcgagaaagtcagcaggccgcttagctagccgctgcctattttgcatgcgggatgt     >  2:151154‑M1/1‑116 (MQ=255)                                                                                                                                                                                       cTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcgggatgtgac  <  2:407784‑M1/300‑188 (MQ=255)                                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                            AGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCG  >  NC_000913/3708172‑3708653

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