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	42 [0]	[24] 66	[proK]	[proK]

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

AAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCT  >  NC_000913/3708171‑3708650                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                         aaGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGTTCGCGCCGCATCCGACAATGGCCTGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaa                                                                                                                                                                                      <  2:336627‑M1/300‑5 (MQ=255)      ttAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctcc                                                                                                                                                                                 <  1:592405‑M1/300‑10 (MQ=255)          gCAAATATAAAGAAAGACAATTGTTTCTTTCACAAATCCGTTAATGGATAGTTAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgacc                                                                                                                                                                             <  2:380596‑M1/300‑14 (MQ=255)            aaagataaaaaAAGAAAATGGTGTCTGTCAAAAATCCGGTAATGGATAGTGAGATAGGGGGCGCAAATGTGGCAGGGAAATGTGTGGCGTGTAGG‑GCGAATAAGGCGGTCGCGCCGCAGCCGACAATGGCCGGCACATTTGACTGAGGCGACGCTGACGCGTCTTAGCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCATAAGGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgacccc                                                                                                                                                                           <  1:198792‑M1/296‑16 (MQ=255)            aaaTATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgacccc                                                                                                                                                                           <  1:295942‑M1/300‑16 (MQ=255)            aaaTAGAAAAAAAGACAATGGTTCCTGTCACAAACCCTTTAATGGATAGTGAAATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGT‑TGGGTTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTACCAGGCCTACGTCGGGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGAGGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgacccc                                                                                                                                                                           <  1:550372‑M1/300‑16 (MQ=255)                   aaaaaaGACAATGGTTTCTGTCACAAATCCATTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcc                                                                                                                                                                    <  2:199758‑M1/300‑23 (MQ=255)                        agacgatcgtttctgtcACAAATCCGTTAATTGATAGTGAGATAGGGGGCGCAAATTCGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACGTTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTcgaacctccgaccccttcgtcccgaac                                                                                                                                                               <  2:176883‑M1/292‑28 (MQ=255)                                 tttCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgc                                                                                                                                                      <  1:84706‑M1/300‑37 (MQ=255)                                         acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                              <  1:593425‑M1/300‑45 (MQ=255)                                         acaAATCCGTGAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                              <  2:177392‑M1/300‑45 (MQ=255)                                           aaagacgttaaTGGATAGTGAGATATGGGGCGCAGATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAAGAATGCGGTCGCGCCGCATCCGACAATTGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccatgctg                                                                                                                                            <  2:62823‑M1/295‑47 (MQ=255)                                           aaaTCCGTTAATTGATAGTGAGATATGGGGCGCAAATTTGGCAGGTAAATGTGTNTCGTGTAGT‑TCGAATAAGGCGTTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctg                                                                                                                                            <  2:514876‑M1/300‑47 (MQ=255)                                              tCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgc                                                                                                                                         <  1:312846‑M1/300‑50 (MQ=255)                                                   tAATGGATAGTGAGATATGGTTCGCAAATTTGGCAGGTAAATGTGTGTCGTGTAGG‑TCGAATAAGGCAGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGGCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGCCTGTGATAGAGGATTcgaacctccgacccctccgtcccgaacggagtgcgctaccatgctgcgccaacc                                                                                                                                    <  2:521786‑M1/300‑55 (MQ=255)                                                            gTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgc                                                                                                                           >  2:415667‑M1/1‑237 (MQ=255)                                                                gaTATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcgggg                                                                                                                       >  1:180908‑M1/1‑233 (MQ=255)                                                                 acatggGGCACAAATTTGGCAGGGAAATATGTGTCGTGTAGG‑TCGAATAATGCGGTCGCACCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCATGCCTACTTCGAGTCGAACCTCCCTTCCCCCGTAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggc                                                                                                                      <  2:198188‑M1/298‑69 (MQ=255)                                                                          gcAAATTTGGCAGGGAAATGTGTGTCGTGTAGT‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttac                                                                                                             <  1:154070‑M1/300‑78 (MQ=255)                                                                          gcAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttac                                                                                                             <  1:214024‑M1/300‑78 (MQ=255)                                                                            aaaTTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactg                                                                                                           <  1:517458‑M1/300‑80 (MQ=255)                                                                            aaaTTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCACGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactg                                                                                                           <  2:3996‑M1/300‑80 (MQ=255)                                                                               tttGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACTTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgc                                                                                                        <  1:596851‑M1/300‑83 (MQ=255)                                                                                    cAGGGAAATGTGTGTCGTGTAGT‑TCGAATAAGGCGGTCGCGCCGCATCCGACAANGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagata                                                                                                   <  2:516171‑M1/300‑88 (MQ=255)                                                                                    cAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagata                                                                                                   <  1:238329‑M1/300‑88 (MQ=255)                                                                                         aaaTGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgctgggcgcatcttactgcgcagatacgccc                                                                                              <  1:149593‑M1/300‑93 (MQ=255)                                                                                                    gtgtAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcca                                                                                   >  1:59866‑M1/1‑197 (MQ=255)                                                                                                          g‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaata                                                                             >  1:440397‑M1/1‑191 (MQ=255)                                                                                                                                   tCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgag                                                     >  1:557310‑M1/1‑167 (MQ=255)                                                                                                                                                             gATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagcc                           <  1:183277‑M1/300‑160 (MQ=255)                                                                                                                                                               tGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgc                         <  1:237551‑M1/300‑162 (MQ=255)                                                                                                                                                                        gACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttt                >  1:58938‑M1/1‑130 (MQ=255)                                                                                                                                                                            cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcc            <  1:59484‑M1/300‑175 (MQ=255)                                                                                                                                                                            cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcc            <  2:345988‑M1/300‑175 (MQ=255)                                                                                                                                                                            cgTCTTATCAGGCCTAAGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgtggggcgcatcttcctgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcggccggccgcttagttagccgctgactcttttgca            >  2:338211‑M1/1‑126 (MQ=255)                                                                                                                                                                             gTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtccagaacgaagtgcgctaccaggctgcgccaatcaacgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggtcgcttagttagccgctgcctcttttgcct           >  2:134985‑M1/1‑125 (MQ=255)                                                                                                                                                                                 tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgngccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctagtcaatcccttaatagcaaaatgtcttttgatcggcgagaaagtcagcaggcagcttagttagccgctgcctcttttgcctgcga       >  2:158526‑M1/1‑121 (MQ=255)                                                                                                                                                                                      ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgcacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgtatcttactgcgcagatacgccctcgtcaatcccttactcgcaaaatgccttttgatcggcgagaaagtcagcaggccacttagttagccgctgcatcttttgcctacgagatgt  >  1:22068‑M1/1‑116 (MQ=255)                                                                                                                                                                                      ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccctttgatcggcgagaaagtcagcaggccgcttagttggccgctgcctcttttgcctgcgagatgt  >  2:4235‑M1/1‑116 (MQ=255)                                                                                                                                                                                      ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccctaatagcaaaatgccttttgatcggcgagaaggtcagcaggccgcttagttagccgctgcctcttttgcctgcgagatgt  >  2:78051‑M1/1‑116 (MQ=255)                                                                                                                                                                                      ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatctnactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcgagatgt  >  1:507733‑M1/1‑116 (MQ=255)                                                                                                                                                                                      ggCCTAAGTCGCGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGTAAAAAAGCCTGCTCGTTGAGTAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtccagaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttaatgcgcagatacgccctcgtcaatcccttaatagccaaatgtcttttgatcggngagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcnagatgt  >  2:574674‑M1/1‑116 (MQ=255)                                                                                                                                                                                                                                                                                                         |                                                                                                                                                                                         AAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGG‑TCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCT  >  NC_000913/3708171‑3708650

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