Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,720,136	IS3 (–) +3 bp :: +TCA	coding (88‑90/213 nt)	cspA →	cold shock protein CspA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	315229 =	NA (NA)	20 (0.380) +TCA	19/480	1.1	100%	noncoding (1/1255 nt)	IS3	repeat region
	?	NC_000913	3720136 =	0 (0.000)					coding (88/213 nt)	cspA	cold shock protein CspA
*	?	NC_000913	= 2171428	NA (NA)	51 (0.970)	50/484	0.1	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	= 3720138	0 (0.000)					coding (90/213 nt)	cspA	cold shock protein CspA

GCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCC  >  NC_000913/2171199‑2171653                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  gctgctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGGGGCGGTCTCAGTCCGGAACAATTGGAAAACAAGAACCTCGCTTAGTTCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCAt                                                                                                                                                                                                               <  1:286059/250‑1 (MQ=255)  ctgctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGCAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAAAAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCATCGTGGATCAttt                                                                                                                                                                                                             >  2:467350/1‑251 (MQ=255)    gctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGAGTCCATATTACGTGGGCAGGATCAGTCATCCCTTCGTGGATCATTTTg                                                                                                                                                                                                           >  1:22405/1‑251 (MQ=255)    gctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACCTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGTCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTGAGGCCTTTGTCCATATTACGTGGGTATGAACAGTCATCCCTTCGTGGATCATTTTg                                                                                                                                                                                                           <  2:456773/251‑1 (MQ=255)      tACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGAc                                                                                                                                                                                                         >  2:150559/1‑251 (MQ=255)       aCGATAATGCCTGCGTGGAAAGCTTCTGTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCGGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACt                                                                                                                                                                                                        <  2:430776/251‑1 (MQ=255)             aTGCCTGCGTGGAAAGCTTCTTTCCTTCGCTGAAAGTGGAATGTCTCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTCATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTAAGTCCGGAACAATTTTATAACAAGAACATCGCATAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTccgggctaattttgaatgccccc                                                                                                                                                                                                  >  2:249027/1‑229 (MQ=11)             aTGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTgccgcc                                                                                                                                                                                                  <  1:442182/251‑1 (MQ=255)                       ggAAAGCTTCTTTCATTCGCTGACAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGGAAAATTTGAAAACAAGAACCTCGCTTAGGACTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTAttt                                                                                                                                                                                        >  2:253231/1‑250 (MQ=255)                       ggAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGCAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGGAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCCGCGCTAttt                                                                                                                                                                                        >  2:60531/1‑250 (MQ=255)                        gaaagcgtccgtcagtCGCTGAAAGTGGAAGGGATCCATGGTAAACCCTTGATGAGCGGGGAAATAATGCGGGGAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGTCGGCACAGTTGGTATGGCGGCGTCAGTCCGGAACAATTTGAAAACAAGATCCTCGCTTAGGCCAGTGTCCATATTCCGTGGCTAGGATCAGTCAGCCCTTCGTGGATCATTTTGACTCCCGCCTGCGCTATTCt                                                                                                                                                                                       <  2:91486/236‑1 (MQ=14)                                  ttCATGCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTGTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCATTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTg                                                                                                                                                                               <  1:253231/249‑1 (MQ=255)                                         gCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATGAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTgg                                                                                                                                                                      <  1:456260/251‑1 (MQ=255)                                            gAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCgg                                                                                                                                                                   >  1:10258/1‑250 (MQ=255)                                            gAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGACTCATTCCGGAACCATTTGAAAACAAGTACCTCGCTTAGGCCTGTGGCCATTTTACGGGGGTAGGATCAGTCATCCCTTCCTGGATCATTTTGACTGCCGCCTGCGCTTTTCTGTTTTTTGGGGCGCTGGCgg                                                                                                                                                                   >  2:283128/1‑250 (MQ=17)                                               aGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCt                                                                                                                                                                <  2:22405/251‑1 (MQ=255)                                                gTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGGAAACAAGAACCTCGCTTAGGCCTGGGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGGGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGGGCTg                                                                                                                                                               >  2:345471/1‑251 (MQ=255)                                                           ccATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCtt                                                                                                                                                    <  2:360951/251‑1 (MQ=255)                                                                gagaACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAAc                                                                                                                                               >  1:153006/1‑251 (MQ=255)                                                                          tATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGTGCTCAGTCCGGAACAATTTGAAAACAAGAACCGCGCTGAGGCCTGTGTCCATATTACGTGTGTAGGATCAGTCGTCCCTTCGTTGATCATTTGGACTGCCGCCTGCGCTATGCTGCTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTcc                                                                                                                                      <  1:434284/250‑1 (MQ=18)                                                                               gCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGtc                                                                                                                                <  1:212460/251‑1 (MQ=255)                                                                                  gggAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTCGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTTAACCCGACTTCCTGGGGCTCg                                                                                                                             >  2:122579/1‑251 (MQ=255)                                                                                         aatGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCAc                                                                                                                      <  2:10258/251‑1 (MQ=255)                                                                                                             aTTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTc                                                                                                  <  2:352293/251‑1 (MQ=255)                                                                                                                       aTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGGGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGc                                                                                        <  1:299796/251‑1 (MQ=255)                                                                                                                         gtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGcc                                                                                       <  1:283128/250‑1 (MQ=255)                                                                                                                                     ggTGGCGGCGGCACAGATGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGAATGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTCTTTTCGCCCCACCGCCATGACAGCGCGt                                                                           >  2:157765/1‑250 (MQ=255)                                                                                                                                        ggcggcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGACTGTGTCAATATTACGTGGGGAGGATCAGTCATCCCTTCGTGGATAATTTTGACTGCCGCCTGCGATATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCCCCGCCATGAAAGCGCGTAAcc                                                                       >  1:317509/1‑251 (MQ=255)                                                                                                                                          cggcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCCCCGCCCTGAAAGCGCGTAACCAc                                                                     >  2:405276/1‑251 (MQ=255)                                                                                                                                           ggcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACg                                                                    <  1:372519/251‑1 (MQ=255)                                                                                                                                                   aGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTAGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGACATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTg                                                            >  1:248751/1‑251 (MQ=255)                                                                                                                                                    gTTGGTGGGGCGGCCTCAGGCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGGGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCGGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCGGGGTCTGGGTTTCACTATTCCCAGACATGATTTTCGCCGCACCGCCATGAAAGCGCGTAACCTCGCCTTTTTGt                                                           <  2:82463/251‑1 (MQ=38)                                                                                                                                                                          ggAACAATTTGAAAACAAGATCCTCGCTTAGTCCTGTGTCCATATTACGGGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAATTCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCAc                                     <  1:330391/251‑1 (MQ=255)                                                                                                                                                                          ggAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCAc                                     <  1:150559/251‑1 (MQ=255)                                                                                                                                                                                                       gggCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGTATCATTTTGACCGCCCCCTGCGCTATTCTGTTTTTGGGCGCGCTGGCGAGCTGAAATCGTTCTTTAACCCCGACTGCCTGGGTCTCGGTTTCCCTATTACCAGACATTATTTTCGCCGCTCCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACt        <  1:467350/250‑1 (MQ=255)                                                                                                                                                                                                       gggCCTGGGTCCATATTGCGGGGGTAGGATCAGTCATCCCTTCGTGGATCATTTGGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTAAACCCCGCCTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCCCACCGCCATGAAAGCGCGTAACCACGCCTGTTTGTTCGAGAAAGCGCAAATCGGCACGGATTTTCGCTTCCGAGGCAGCAAATACt        <  2:304406/250‑1 (MQ=255)                                                                                                                                                                                                       aGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAACCCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATAc         <  1:60531/250‑1 (MQ=255)                                                                                                                                                                                                             gggcccATATTGCGGGGGTAGGATCAGCCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAACCCCGACTTCCTGGTTCTGGGTTCCACTATTACCAGACATTTTTTGCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGcc  <  2:61829/247‑1 (MQ=255)                                                                                                                                                                                                                  cATATTCCGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTCTTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGcc  >  2:230089/1‑246 (MQ=255)                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                  GCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCC  >  NC_000913/2171199‑2171653

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 38 ≤ 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.