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.620)	30/486	0.6	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)	46 (0.890)	45/486	0.2	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	= 3720138	0 (0.000)					coding (90/213 nt)	cspA	cold shock protein CspA

ACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGT  >  NC_000913/2169927‑2170410                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                aCCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCt                                                                                                                                                                                                                                            <  1:420061/250‑1 (MQ=255)   cATCAGCACCGAAAAACATTAAATCGGCATTAATATCCTATAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTAc                                                                                                                                                                                                                                         >  1:441268/1‑251 (MQ=255)           ccGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGAACAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAAt                                                                                                                                                                                                                                 >  1:33893/1‑251 (MQ=255)              aCAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAACATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAgg                                                                                                                                                                                                                              >  1:50300/1‑250 (MQ=255)                   cATTAAATCGGCATTAATATACTGTAATGAACGCTCGGCACTGGAACCATGAATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGA‑CTTTATTTTCGGAAAGGGCAAACGCTGACGGGAGACTGTTGGTCATCACCGTGATATTTTTGGCGGTCATTCATCCTTCGGCGATAAGCATTGTCGCACTTCCGATGTTGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGac                                                                                                                                                                                                                          >  1:317320/1‑250 (MQ=255)                        aaTCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCCGGGAAATGA‑CTTTTTTTTCGGAAAGGGCAAACGCTGCCGGGGGGGTGGTGGTAATCACCGTGATATTTTTGGCGGTCCTTGATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGGGAATAACGGTCTGGTCTTTCCCCGCAATAGTgggaccgcgg                                                                                                                                                                                                                    >  2:215243/1‑243 (MQ=255)                         aTCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGGAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCCCCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCAATGTCGTACTTCCGCTGTCGAGAATAACGGACTGATCTTACCCAGCAATAGTGGACAAGCggg                                                                                                                                                                                                                   >  1:516386/1‑250 (MQ=255)                                 tAATATCCTGTAATGAACGATCGGCAATGGAACCATGCATCAAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGCCCAGGGTCATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCAGGGAGATTGTTGGTAATCACCGTGGTATTTTTGGCGGTGATTAATTATTCGGCGGTAAGCAGTGTCGTACTTCCTCTGTCGAGAATAACGGTGTGATCTTACCCAGCAATAGTGGACACGCGGCTCAGTg                                                                                                                                                                                                             >  1:509218/1‑249 (MQ=255)                                     atCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGacac                                                                                                                                                                                                                        <  2:292532/234‑1 (MQ=255)                                     atCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGacac                                                                                                                                                                                                                        >  1:292532/1‑234 (MQ=255)                                            aaTGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGCGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACtctc                                                                                                                                                                                                  >  1:259242/1‑249 (MQ=255)                                               acctcattgacggcatggccATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTTTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGTCTAAGTGAGTAAACTCtcagtc                                                                                                                                                                                              <  1:459175/232‑1 (MQ=255)                                                   gCTCGGCAATGGAACCATGCATCGAGCGCGTTGTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAgg                                                                                                                                                                                         <  2:188415/251‑1 (MQ=255)                                                           acggAACCATGCATCGAGCGCGTTTTATGGCGTACCGTGCCACCACAGACGACCAGAGTAAGGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAAGCTCCGTGATATTTTTGGCGGTCATTAATCCTTCAGCGATAAGCATTGTCGTGCTTCCGCTGTCGAGAATATCGGGCTGATCTTACCCAGCAATAGTGGACACGCGGCGAAGTGAGTAAACTCGCAGTCAGAGGTTACTcac                                                                                                                                                                                 <  2:490775/249‑1 (MQ=255)                                                                aCCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGAc                                                                                                                                                                            <  2:170306/251‑1 (MQ=255)                                                                       aTCGAGCGCGTTTGATCGCGGGCCGTGCCACCTCAGACGACCAGAGTA‑TGTTCGTTAGTTTCGGAAAGGGCAAGCGTTGTCGGGAGACTGTGTGTAATCTCCGTGATATTTTGGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGCCTGATCTTACCCAGCAATTGTGGACACTTGGCTAAGTGAGTAATCTCGCAGTGAGATGTGACTCACTTTACTAAAACa                                                                                                                                                                     <  2:481204/251‑1 (MQ=255)                                                                          gAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGCGAAATGT‑ATTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACAGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATATTACCCAGCAATAATGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGt                                                                                                                                                                   >  1:481864/1‑250 (MQ=255)                                                                              gcgTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                                                                                                                                               >  1:352546/1‑234 (MQ=255)                                                                              gcgTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACAt                                                                                                                                                                               <  2:352546/234‑1 (MQ=255)                                                                                    tATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCGTTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGt                                                                                                                                                         <  1:238530/250‑1 (MQ=255)                                                                                                             gacCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGGCGTACTTCCGCTGTCGAGAAGAACGGTCTGATCTTTCCCAGCAATAGTGGACACGCGGCTAAGTGATTGAACGCTCAGTCAGAGGTGACTCAGAGGACAATAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGcc                                                                                                                               <  1:352354/251‑1 (MQ=37)                                                                                                             gacCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATAGTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCgac                                                                                                                               <  1:336557/251‑3 (MQ=255)                                                                                                                            cTTTATTTTCGGAACGGGCACACGCTGCAGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAACAAACCCCGTAAACAGCATTCGCCTGACTTTCGCAGTg                                                                                                                 >  2:63009/1‑251 (MQ=255)                                                                                                                                      ggAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGGGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGaa                                                                                                       <  2:314425/251‑1 (MQ=255)                                                                                                                                        aaaGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTGAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGc                                                                                                     >  1:379051/1‑251 (MQ=255)                                                                                                                                                          ggAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCt                                                                                                    <  1:538459/234‑1 (MQ=255)                                                                                                                                                          ggAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCt                                                                                                    >  2:538459/1‑234 (MQ=255)                                                                                                                                                                                gATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACt                                                             >  2:109650/1‑251 (MQ=255)                                                                                                                                                                                      tttGGCGGTCCTTAATCCTTCAGCGATAAGCATTGTCGTACTTCCGCTGCCGACAAAAACGGTCTGATCTTACCCAGCAATCGTGGCTACGCGGCTAAGTGAGGAACATCTCAGTCAGAGGAGACTCACATTTCAAAAACAGTATCAACCAGTAAAAAACCCCGTACACAGCACTCGCCTGAATTTCTCAGTGTAGCACTGAAGCTTGCTGCCCGCATCGGTGTTTCTGCCGCacgccgcgcaatcagcct                                                       >  2:486530/1‑234 (MQ=12)                                                                                                                                                                                      tttGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTAGCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGcc                                                        <  2:50300/250‑1 (MQ=255)                                                                                                                                                                                               cATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATc                                              >  1:169838/1‑251 (MQ=255)                                                                                                                                                                                                aTTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATc                                              >  2:253136/1‑250 (MQ=255)                                                                                                                                                                                                      ccTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAActc                                       <  1:28436/251‑1 (MQ=255)                                                                                                                                                                                                             cGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACt                                <  2:394822/251‑1 (MQ=255)                                                                                                                                                                                                               aTAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTgg                              >  1:408792/1‑251 (MQ=255)                                                                                                                                                                                                                       tGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTa                       >  2:280398/1‑250 (MQ=255)                                                                                                                                                                                                                        gTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGCCACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTaa                      >  2:273191/1‑250 (MQ=255)                                                                                                                                                                                                                                  cGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGCCACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACACAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGTTGTTACTGCCGCAGCCCGTGAACTCAGCCGGTATGAATCACAACTCTACAACTGGCACAGTCAACAGCAAAATc           >  2:68884/1‑251 (MQ=255)                                                                                                                                                                                                                                         gagaATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGAc    >  1:165702/1‑251 (MQ=34)                                                                                                                                                                                                                                           gaATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTAAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCCTTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGt  >  1:75952/1‑251 (MQ=18)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                ACCATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGT‑CTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAAACAGCAAAATCAGCAGACGT  >  NC_000913/2169927‑2170410

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