Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,468,298	IS3 (+) +3 bp	coding (85‑87/921 nt)	yfdH →	bactoprenol glucosyl transferase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2170171 =	NA (NA)	41 (0.790)	38/492	0.4	95.3%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	= 2468300	2 (0.040)					coding (87/921 nt)	yfdH	bactoprenol glucosyl transferase
*	?	NC_000913	= 2171428	NA (NA)	42 (0.810)	37/490	0.4	95.4%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	2468298 =	2 (0.040)					coding (85/921 nt)	yfdH	bactoprenol glucosyl transferase

ACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGTCTTCTG  >  NC_000913/2169927‑2170416                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       cccaacAGCACCGACAATCAATAAATCGGCATTAATCTCCTGTAAAGACCGCTCGGCAATGGAACCAGTCATCGATCGCGTTTTATGGCGGTCCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTGTCGTAAAGGTCAAACGCGGCCGGGAGACTGTTGTGAATCACCTTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAATCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCtt                                                                                                                                                                                                                                                  <  2:500172/246‑1 (MQ=255)            cGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTTGCGGTCATTAATCCTTCGGCGCTAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAAt                                                                                                                                                                                                                                        >  1:495813/1‑250 (MQ=255)                   cATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCTGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCTTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGaca                                                                                                                                                                                                                                <  1:218082/251‑1 (MQ=255)                     ttAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGCCCAGGGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGGGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGactcg                                                                                                                                                                                                                              >  1:371140/1‑248 (MQ=255)                                       ccTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAAc                                                                                                                                                                                                            <  1:403435/250‑1 (MQ=255)                                       catGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGAAGACCAGAGGAATGTCTTTATTTTCGGCAAGGGGAAACGCTGGCGGGAGACTGTTGGTAATCAACGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGAGTTGTCGTACTTCCGCTGTGGAGAATAACGGTATGCTTTTAACCAGCAATGGAGGACCCGAGGC‑TAAGTAAGTaaatt                                                                                                                                                                                                           >  1:240627/3‑249 (MQ=255)                                          gTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAATGGACACGCGGC‑TAAGTGAGTAAACtctc                                                                                                                                                                                                        >  1:266667/1‑251 (MQ=255)                                               gAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCACACGCTGCCGGGAGACTGTTGGTAATCAACGTGATATTTTTGGCGGGCATTAATCCTTCGGCGATAAGCATTTTCGTACTTCCGCTGTCGGAAAAAACGGTCTTATTTTACCCCGCCATAGGGGACACGCGGC‑GAAGAGAATaaaactccatgcc                                                                                                                                                                                                   >  2:435889/1‑241 (MQ=255)                                                aaCGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGACACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTTTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCtcagtc                                                                                                                                                                                                    >  1:115817/1‑249 (MQ=255)                                                  cGCTCGGCAAGGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCagag                                                                                                                                                                                                <  2:189053/251‑1 (MQ=255)                                                       ggCAATGGAACCATGCAACGGGCGCGTTTTATGGCGGACCGTGCCTCCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCGTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGAc                                                                                                                                                                                           <  2:247834/251‑1 (MQ=255)                                                              gAACCATGAAACGATCGCGTTTTATGGCGGGCCGTGCCCCCACATACTACTATAGTAATGGCTTTATTTTCGGAAAGGGCAAACGCTGCCGGAAATCTGTCAGCAATCACCGTGAGATTTTAGGCTGTCATTAACCCTTCGGCGATATCCATGGTCGTACTTCCGCTGTCTAGAATAACGGTCTGATCCTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACACTCAGTCAGAGGTGATTCACATg                                                                                                                                                                                    <  2:442500/251‑1 (MQ=19)                                                                aCCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACCAACGAACAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGGGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCCTTTTCGTACTTCCGCTGTCGAGAATAACGGTCTTATCTTACCCCGCAATTGTGGACACGCGGCATAAGTG‑GAAATATCTCAATCAAAGGTGACTCACAtaac                                                                                                                                                                                  >  2:368444/1‑248 (MQ=255)                                                                   aTGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTTGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGAATTGTCGTACTTCCGCTGTCGAGAATACTGGTCTTATCTTACCCAGCAATAGTGGACACGAGGC‑TAAATTAGTAAACCCTCAGTAAGAGGTGACTCCAATGACaa                                                                                                                                                                                >  1:132248/1‑250 (MQ=255)                                                                       aTCGAGAGCGTTTTATGGCGGACCGGGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGAATAGCATTGTCGTACTTTCGCTGTCGAGAATACAGGACTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAAGGAGTAAACTCTCCATCAGAGGTTGCTCAAAAGACAAAAAc                                                                                                                                                                            >  2:97008/1‑250 (MQ=255)                                                                                       ggCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTaaa                                                                                                                                                            <  1:151587/250‑1 (MQ=255)                                                                                                              acCAGAGTTATGTCTTTATTTGCGGAAAGGTCAAACGCTGCCGGGAGACGGTTGGTAATCGCCGTGATATTTTTGGCGGTTAGTAATCCGTCGGCGATAAGCATTGTCGTACTTCCGCTGGCGAGAATAACGGTCTGATCTGACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGGGACTCACATGACAAAAACAGTATCAACCAGTAAAAGACCCCGTAAACAGCATTCGccc                                                                                                                                    <  2:15256/251‑2 (MQ=21)                                                                                                                   agTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAAtt                                                                                                                               <  1:435889/251‑1 (MQ=255)                                                                                                                   agTAATGTCTTTAGTTTCGGAAAGGGCCAACGCTGCCGGGAGTCTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATGGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAGCCCCGTAAACAGCATTCGCCTGAAtg                                                                                                                               <  2:371140/251‑2 (MQ=255)                                                                                                                                 tttCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCCGTCATTAATCCTTCGGCGATAAGCATTGACGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAAGTGACTCACATGACAAAACCAGTATCAACCAGTAAAGAACCCCGTAAACAGCCTTCGCCTGAATTTCACAGTGAAGccc                                                                                                                 >  1:347214/1‑251 (MQ=25)                                                                                                                                       aaaGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGc                                                                                                           >  1:193068/1‑251 (MQ=255)                                                                                                                                        aaGGGCAAACGCTGCCGGGAGACTGTTGGTAATCCCCGTGATATTTTTGGCGGTCATTAATCCTGCGGCGATAAGCATTGGCGGACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCGCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCt                                                                                                          <  2:58519/251‑1 (MQ=255)                                                                                                                                                       cGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCCGTGAAGCCCGGAAGCTTGCTGAACGCATCgg                                                                                           >  2:192185/1‑251 (MQ=255)                                                                                                                                                          gagaCTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGgtgt                                                                                        <  1:479862/251‑1 (MQ=255)                                                                                                                                                            gaCTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGtt                                                                                       >  1:64659/1‑250 (MQ=255)                                                                                                                                                               ggTTGGTAATCACCTTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTATTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGATGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTg                                                                                   <  2:379955/250‑1 (MQ=255)                                                                                                                                                                gTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGAAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTCCTGc                                                                                  >  2:318778/1‑251 (MQ=255)                                                                                                                                                                         aaCCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATGGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCATAGGTGACTCACATGACAAAAACAGTATCACCCAGTAAAAATCCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGccc                                                                          <  2:326261/249‑1 (MQ=255)                                                                                                                                                                             gTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGAATCGGTGTTACTGCCGCAGCCCGTGa                                                                      >  1:157689/1‑250 (MQ=255)                                                                                                                                                                                      ttGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCt                                                             <  1:503581/250‑1 (MQ=255)                                                                                                                                                                                        ggCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAACCCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGACCTCAGCCTGt                                                           >  1:240032/1‑250 (MQ=255)                                                                                                                                                                                                      cTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCCTCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAActcc                                            >  1:69201/1‑250 (MQ=255)                                                                                                                                                                                                              aTAAGCATTGTCGTTCTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTgg                                    <  1:420712/251‑1 (MQ=255)                                                                                                                                                                                                                aaGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACTCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGgcg                                  <  2:142202/251‑1 (MQ=255)                                                                                                                                                                                                                                    tGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATcagc              >  2:219115/1‑251 (MQ=35)                                                                                                                                                                                                                                        gagaATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGAc          <  2:254541/251‑1 (MQ=34)                                                                                                                                                                                                                                         agaATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACg         <  2:113499/251‑1 (MQ=33)                                                                                                                                                                                                                                                 ggTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACCACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGTCTTCTg  >  2:131638/1‑250 (MQ=11)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       ACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGC‑TAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGTCTTCTG  >  NC_000913/2169927‑2170416

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.