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	41 [0]	[24] 65	[proK]	[proK]

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 3708466	0 (0.000)	41 (0.790)	28/374	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

AAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTC  >  NC_000913/3708168‑3708652                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                           aaaaaGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCNAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTc                                                                                                                                                                                           <  1:124856‑M1/300‑2 (MQ=255)      gATTGAAGCAAATATAACACAAGGCAATGGTTTCTTTCGCAAGTCCGATAATGGATAGTTAGATGTGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAATGCAGTCGCGCCGCATCCGACAACGGCCGGCACATTTGCCTGATGCGACGCTTACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGTAAGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGGTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTcgaacc                                                                                                                                                                                      <  2:317927‑M1/300‑7 (MQ=255)       aTTTAAGCACATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacc                                                                                                                                                                                      <  2:317700‑M1/299‑7 (MQ=255)                      aaaaaaGACAATGTTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCACAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcc                                                                                                                                                                      <  1:6624‑M1/300‑23 (MQ=255)                      aaaaaaGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcc                                                                                                                                                                      <  2:582340‑M1/300‑23 (MQ=255)                           aGACAATGTTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaac                                                                                                                                                                 <  1:499310‑M1/300‑28 (MQ=255)                           aGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaac                                                                                                                                                                 <  1:499261‑M1/300‑28 (MQ=255)                                            acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGTAAATGTGTGTCGTGTAGTTCGAATAAGGCGTTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACTTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                <  2:157639‑M1/300‑45 (MQ=255)                                            acaAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                <  1:388312‑M1/300‑45 (MQ=255)                                            acaAATCCGTTAATGGATAGGGAGATATGGTGCACAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTNTTATCAGGCCTACTTCGAGTCGAACCTCCCTGCCCCCGTAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                <  2:199771‑M1/300‑45 (MQ=255)                                              aaaTCCGTTAATGGATATTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggc                                                                                                                                                <  2:152737‑M1/298‑45 (MQ=255)                                                 tCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCNCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgc                                                                                                                                           >  2:502754‑M1/1‑251 (MQ=255)                                                         tGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcacc                                                                                                                                   >  1:401720‑M1/1‑243 (MQ=255)                                                                 gagaTATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgc                                                                                                                             >  1:153659‑M1/1‑235 (MQ=255)                                                                 gagaTATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtccagaacgaagtgcgctaccaggctgcgccaatcaccgaatgc                                                                                                                             >  1:448756‑M1/1‑235 (MQ=255)                                                                   gaTATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCNTGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGNGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcgggg                                                                                                                         <  2:495964‑M1/300‑68 (MQ=255)                                                                   gaTAGGGGGCGCAAGTTTGGCAGGGAAATGTGTGTCGTGTAAGTCGAAGAAGGCGTTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACACTGACGCGTCTTATCAGGTCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgactcctccgtcccgaacg                                                                                                                                                                <  2:46132‑M1/261‑29 (MQ=255)                                                                    atatGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggc                                                                                                                        <  2:652515‑M1/300‑69 (MQ=255)                                                                                  tttGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGNTGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgc                                                                                                          <  1:124403‑M1/300‑83 (MQ=255)                                                                                  tttGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgc                                                                                                          <  1:498948‑M1/300‑83 (MQ=255)                                                                                           gAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgcc                                                                                                 >  1:625407‑M1/1‑209 (MQ=255)                                                                                            aactgtgtgtctTGTAGTTCGCATAATGCGTTCGCGCCACATCCGACAATGGTCGGCACATTTGCCCGATGCGACGCTGACGCGTCTCATCAGGCCTACGTCGAGGCAAACCTCCCTGCCCCCGTAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTGGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaactaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccc                                                                                                <  2:340920‑M1/297‑93 (MQ=255)                                                                                                  gtgtCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtca                                                                                          >  1:577797‑M1/1‑202 (MQ=255)                                                                                                     tCGTGTAGGGTGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATTCGACGCTGATGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGTTGATAGATGATTcgaacctccgaccccttcgtcccgaactaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgc                                                                                                          <  1:171475‑M1/281‑83 (MQ=255)                                                                                                      cGTGTAGATCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccacttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaat                                                                                        >  1:160542‑M1/1‑198 (MQ=255)                                                                                                       gtgtAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgactgcggggcgcatcttaatgcgcagataagccctcgtcaatccc                                                                                     >  2:176675‑M1/1‑197 (MQ=255)                                                                                                       gtgtAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatc                                                                                       >  2:178241‑M1/1‑197 (MQ=255)                                                                                                            ggTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaat                                                                                >  1:153542‑M1/1‑192 (MQ=255)                                                                                                                                                  ggCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcag                                            <  1:82985‑M1/298‑145 (MQ=255)                                                                                                                                                           gCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagtt                                 >  1:31282‑M1/1‑145 (MQ=255)                                                                                                                                                           gCCTGATGCGAAGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtccagaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggtcgcttagtt                                 >  1:35493‑M1/1‑145 (MQ=255)                                                                                                                                                               gATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagcc                             <  1:86603‑M1/300‑160 (MQ=255)                                                                                                                                                                          gACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttt                  <  1:329361‑M1/300‑171 (MQ=255)                                                                                                                                                                              cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcc              >  1:313437‑M1/1‑126 (MQ=255)                                                                                                                                                                              cgTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTG‑AGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagntacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcct             >  2:503715‑M1/1‑125 (MQ=255)                                                                                                                                                                                  ttATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAACCTGCTCGTTGAGCAGGTTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacggagtgtgctaacaggctgcgcccaccaacgaatgcggggcgcatcttactgcgcagatacgccctcgtcaaccccctaatagcaaaatgactgttgatcggtgagacagtcagcaggccgcttagttcgccgctgcttcttttacatgc           >  2:14509‑M1/1‑122 (MQ=255)                                                                                                                                                                                   tATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgctcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgc           >  1:314265‑M1/1‑121 (MQ=255)                                                                                                                                                                                        ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcgggatgt    >  1:86529‑M1/1‑116 (MQ=255)                                                                                                                                                                                        ggCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgcattagtcaatcccttaatagcaaaatgcctcttgatcggcgcgaaagtcagcaggccgcttagttagccgctgcctcttttgactgcggggcgt    >  2:269346‑M1/1‑116 (MQ=255)                                                                                                                                                                                         gCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcgggatgt    >  1:517678‑M1/1‑115 (MQ=255)                                                                                                                                                                                             aCGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTAGGTGATAGAGGATTcgaacctccgaccccttcgtcccgaacgaagtgcgctaccaggctgcgccaatcaccgaatgcggggcgcatcttactgcgcagatacgccctcgtcaatcccttaatagcaaaatgccttttgatcggcgagaaagtcagcaggccgcttagttagccgctgcctcttttgcctgcgggatgtga  >  2:177387‑M1/1‑111 (MQ=255)                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                           AAAAAGATTTAAGCAAATATAAAAAAAGACAATGGTTTCTGTCACAAATCCGTTAATGGATAGTGAGATATGGGGCGCAAATTTGGCAGGGAAATGTGTGTCGTGTAGGTCGAATAAGGCGGTCGCGCCGCATCCGACAATGGCCGGCACATTTGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACGTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTACTCGCCACTTCGTGGCTCGCCCTGCGGGCCGTTGCTGGCGCAACGTTCTCTCGCTTCGCTCGAGTCGAACCTCCCTGCCCCCGGAGGTTCTCATCCTCATGAGCTGCAGATGCAAAAAAGCCTGCTCGTTGAGCAGGCTTTTCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTC  >  NC_000913/3708168‑3708652

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