Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,720,136	IS3 (–) +3 bp	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	2170171 =	NA (NA)	32 (0.700)	31/482	0.4	100%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	3720136 =	0 (0.000)					coding (88/213 nt)	cspA	cold shock protein CspA
*	?	NC_000913	= 2171428	NA (NA)	50 (1.090)	44/482	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

ATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACAA  >  NC_000913/2171211‑2171666                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 aTGCCTGCCTGGAATGCTTCTTTCATTCGCTGAAAGTGGGATTTATCCATGGAGAACACTGTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAAGTATATCGACTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGTCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGGTTAGGCCTGTGTCCATATTACGTGGTTAGGATCAGTCATCCGTGCGTGG‑ATCATTTTAACT‑gccgcc                                                                                                                                                                                                               <  2:321014/251‑1 (MQ=14)   gCCTGCGTGGAAAGCTTCTTGCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCt                                                                                                                                                                                                              <  2:367660/250‑1 (MQ=255)              aaGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTTATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGt                                                                                                                                                                                                  <  1:56070/251‑1 (MQ=255)                gCTTCTTTCATTCGATGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGttt                                                                                                                                                                                                >  1:150077/1‑251 (MQ=255)                               tGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCg                                                                                                                                                                                 >  1:226279/1‑251 (MQ=255)                                 aaaGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCg                                                                                                                                                                                 >  1:187599/1‑249 (MQ=255)                                    gTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCCGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTg                                                                                                                                                                            >  1:373879/1‑251 (MQ=255)                                     tGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATCATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGCGGCGGCGGCACAGTTGGTGTGGCGGCCTCCGGCCGGAACAATTTGAAAACAAGAACCCCGCTTAAGCCTGTGTCCATATTACGAGGGTAGGATCAGTCATCCCTTCGTGG‑ATAATTTTGACT‑GCCGCCTGCGCTATTCTGCTTTTTGGCGCGTTGGCGAGCTg                                                                                                                                                                            >  1:302413/1‑250 (MQ=17)                                      ggAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAAAGGTGATTAATTATATCGAATGTGATTACAATAGGAGGCGGCGGCACAGTTGGTGTGGCGGACTCAGTCCGGAACAATTTGGAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCATGCGCTATTCTGTTTTTTGAAGCGCTGGCGAGCTGa                                                                                                                                                                           >  1:183239/1‑250 (MQ=17)                                       gAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGCGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGCGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTTGCGCGCTGGCGTGCTGaa                                                                                                                                                                          >  2:400347/1‑250 (MQ=255)                                                       aaCACTTTATCAGCCGGTAAATAATGCGGTCAGCGGTGTTTAATTATATCGAATGTTATTACAATCGGTGGCGGCGGCACAGTTTGTGTGGCGGCCTCAGTCCGTAACAATTTGAAAACAAGAACCGCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCAGCCCTTCGTGG‑AGCATTTTGACT‑GCCGCCTGCGCTAGGCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCg                                                                                                                                                         <  2:224364/251‑1 (MQ=18)                                                               aTCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGCTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTg                                                                                                                                                 >  1:338636/1‑251 (MQ=255)                                                                 cAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTgg                                                                                                                                                <  1:441655/250‑1 (MQ=255)                                                                     cGGGACATAATGCGGGCAACGGTGTTTAATTATCTCGAATGTGATTACCATCGGTGGCGGCGGCACAGTTGGAGTGGCGGCCTCAGTACGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGACCATATTCCGTGGGTCGGATCAGTTATCCCTTCCTGG‑GTCCTTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAAACCGACTTCCTGGGtctc                                                                                                                                           >  2:281457/1‑251 (MQ=21)                                                                     cGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTAGTCTGTTGTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGtctc                                                                                                                                           <  1:442383/251‑1 (MQ=255)                                                                              atGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCCTATTACGTGGGTAGGATCAGTCATCCCTTCGGGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGGTGGCGAGCGGAAATCGCTCTTTTAACCCGGCTTCCTGGGTCTCGGttttcct                                                                                                                                  >  2:371914/1‑247 (MQ=21)                                                                              atGCGGGCAACGGTGTTTAAGTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGGGGGTATGATCAGTCATCCCTTCCTTG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCCCTCTGGCCAGCTGAAAACGCTCTTTAACCCCGTCTTCCCGGGTCTCGTTTTCAcc                                                                                                                                  >  2:450785/1‑250 (MQ=18)                                                                               tGCGGGCAACGGTGTTTAATTATATCGAATGTGCTTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGGGTCCATATTACGTGGGTAGGATCAGTCATCCCCTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGTCGCGCTGGCGAGCTGAAATCGCTCTTTACACCCGCCTTCCCGGGGCTCGGTTTctct                                                                                                                                  >  2:100807/1‑247 (MQ=21)                                                                                            gtTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACACGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGGGCTATTCTGTTTTTTTGCGCGCTGGCGGGCTGAAATCGCCTTTTAAACCCGAc                                                                                                                                                       >  2:200038/1‑216 (MQ=21)                                                                                            gtTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGAc                                                                                                                                                       <  1:200038/216‑1 (MQ=255)                                                                                                    atatCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTAGAGGCATCAGTCCGGAACAATTTGAAAACAAGACCCTCGCTTAGGCCTGTGTCCAACTTAAGTGGGTAAGATCAGTAATACCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAAACCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTcgcc                                                                                                            >  1:435069/1‑251 (MQ=37)                                                                                                      atCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCttt                                                                                                                                                                <  1:461911/197‑1 (MQ=255)                                                                                                           aTGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGc                                                                                                     >  1:232768/1‑251 (MQ=255)                                                                                                                                       acttGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGTCTGTGTCCAGCTTACGTGGGTAGGATCAGTCACCCCTACTTGG‑TTAATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTTGCGCGCTGGCGAGCTGACATCGCTCTTTAAACCCGACTTCCTGGGTCACGGTTTCCCTATTACCAGACATTATTTTCGCCGCCCCGCCATGTAACTGCGTCACCACGCCtttttt                                                                         >  2:100532/3‑250 (MQ=25)                                                                                                                                                    gcCTCAGTCCGGAACCATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCCTGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAgcgc                                                            >  2:25331/1‑251 (MQ=255)                                                                                                                                                    gcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGAATAAGTAATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAgcgc                                                            >  1:82270/1‑251 (MQ=255)                                                                                                                                                    gcCTCAGGCCGGGACGCTTTGAAAACAAGAACCTAGCTTAGGCCTGTGTCCACATTACGTGGGTAGGCTCAGTCATCCCTTCGTGG‑CTCATTTTGACT‑GCCGACTGCGCTATTCTGTTTTTTGGCGAGCTGGCGAGCTGAAATCGCTCTTTAAAACCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCCCCGCCATGCAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAgcgc                                                            >  2:299644/1‑251 (MQ=255)                                                                                                                                                      cTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCa                                                           >  2:92513/1‑250 (MQ=255)                                                                                                                                                                    aTTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTg                                            >  1:227588/1‑251 (MQ=255)                                                                                                                                                                       ggAAAACAAGAACCCCGCTTGGGCCTGGGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTCCCTGGGTCTCGGTTCCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTc                                          <  2:232768/249‑1 (MQ=255)                                                                                                                                                                           aaCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑CTCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGCGAAAGCGCAAATCGGCCCGGATTGTCGCTTc                                     >  2:46350/1‑251 (MQ=255)                                                                                                                                                                            aCAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTcc                                    <  1:256191/251‑1 (MQ=255)                                                                                                                                                                             aaaGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGGTAT‑ATTTTGCCT‑GCCGCCTGCGCTATTCTGTTTTTTGGCCCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGCCTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTGGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCg                                   <  1:46350/250‑1 (MQ=255)                                                                                                                                                                                 aaCCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTa                                                                                                                                 >  1:437388/1‑153 (MQ=255)                                                                                                                                                                                 aaCCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCGGGGTCTCGGTTTCACTa                                                                                                                                 <  2:437388/153‑1 (MQ=255)                                                                                                                                                                                           gggCCTGTGCCCATATTACGTGGGTAGGAGAAGTATTCTCTTCGTGG‑‑TTAGTTTGGCTTGCCGCCTGCGCTATTTTGTTTTTTGGCGCGCTGGCGAGCTAAATTCGCTCTTTAACCCCGACTTCCTGGGTCTCGGTTTCACAATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTGTTTGTGCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATAc                      <  2:227588/249‑1 (MQ=255)                                                                                                                                                                                              ccgttgtCCAGATTCCGTGGGTAGGATCAGCCATCCCTCCGTGA‑ATCATTTTGACT‑GCGTCCTGCGGTATCCGGTGTTTTGGCGCGCGGGGGGGTGGAAATCGCGCTTAAACCCCGACTTCCGGGTTCGGGGTTCCTCTATTCCAAGGCATTTTTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTTTTCGAGAAAGCGCAAATCGGCACGGATTGTCTCTACCGAGGCATCAAGTACTccc                  <  2:49808/247‑1 (MQ=255)                                                                                                                                                                                                       atatTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTTGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAACTACTCCCGCCAGATcc         >  2:140208/1‑251 (MQ=255)                                                                                                                                                                                                              gTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGTTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTAAAACCCGACTTCCTGGGTCTCGGTTTCACTATTCCCAGACATTATTTCCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACaa  <  1:29305/251‑1 (MQ=255)                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                 ATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTCATCCCTTCGTGG‑ATCATTTTGACT‑GCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACAA  >  NC_000913/2171211‑2171666

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