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)	27 (0.590) +TCA	26/482	0.7	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)	45 (0.980)	44/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

ATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACAAGCACGGT  >  NC_000913/2171183‑2171673                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        aTGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACCATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑At                                                                                                                                                                                                                                                   >  1:145132/1‑251 (MQ=14) aTGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑At                                                                                                                                                                                                                                                   <  1:253000/251‑1 (MQ=31)         aaaaGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTGAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGt                                                                                                                                                                                                                                           <  2:411532/251‑1 (MQ=25)             ggTTGCTGCTGCGATAATGCCTGCGTGGAAAGCTTCTTTCATGCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATGACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGAt                                                                                                                                                                                                                                       <  2:244039/251‑1 (MQ=255)                tgctgcTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTAGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCCGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCAt                                                                                                                                                                                                                                    >  1:22877/1‑251 (MQ=255)                 gctgctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCAtt                                                                                                                                                                                                                                   <  1:316262/251‑1 (MQ=39)                   tgctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCAtttt                                                                                                                                                                                                                                 >  1:192280/1‑251 (MQ=40)                     ctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATAATTTTg                                                                                                                                                                                                                                >  1:159740/1‑250 (MQ=39)                                   gCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTgcgc                                                                                                                                                                                                                  <  2:269162/250‑1 (MQ=255)                                    cGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAAGGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCCGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAAGTTGAAAACAAGAACCGCGCTTAGGCCTGGGGCCATATTACGTGGGTAGGATCAGTC‑GTCCCTTCGTGGATCATTTTGACCGCCGCCTGCGCTa                                                                                                                                                                                                                <  2:208742/251‑1 (MQ=14)                                     gTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTTTTTAATTAGATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGGGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGCCTGCCGCCTGCGCTAt                                                                                                                                                                                                               <  2:138717/251‑1 (MQ=255)                                            gCTTCTTTCATTCGCTGAAAGTGGCATGTATCCATGGCGCACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCCGACCGGAACAATTTGACAAAAAGAAACTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGAACAGTC‑AGCCCTTCGTGGGTACTTTTGACTGCCGCCTGCGCTCATCTGttt                                                                                                                                                                                                        >  2:64855/1‑251 (MQ=14)                                                 tttCATTCGCTGAAAGTGAAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACCATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATTTTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTgg                                                                                                                                                                                                   >  2:168282/1‑251 (MQ=255)                                                            gAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGc                                                                                                                                                                                          <  1:131549/249‑1 (MQ=255)                                                                         atccagtgagagcccttgaaccgCCGGGAAAGAATGCGGGCAACGGTGTTTAAGTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTGCGTGGATCATTTTGACTTCCGCCTGCGCTATTCTGTTTTTTGGCGCGGTTGCGAGCTGAAAGCGct                                                                                                                                                                            <  2:295153/228‑1 (MQ=17)                                                                         aTCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGctc                                                                                                                                                                           >  1:24348/1‑251 (MQ=255)                                                                                    acacTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCg                                                                                                                                                                 >  1:88854/1‑250 (MQ=255)                                                                                         ttATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTc                                                                                                                                                            <  2:159740/250‑1 (MQ=255)                                                                                           agccggcGGGAAATAATTCGGGCAACGGTGTTTAATTATATCGAATTTGAGTACAATCGGTGGCGGCGGCACAGGTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCAATTGGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCGGAAATCGCTCTTTAAACCCGACTACCt                                                                                                                                                          <  2:337124/244‑1 (MQ=21)                                                                                              aGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGCATGTGATTAAACTAGGTGGAGGCGGCACAGTTGGTGTGGCGGCCTCCGTCAGGAACAATTTGAAAACAAGAACACAGCCTACGCCTGAGTCAATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTCTTTTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGt                                                                                                                                                      >  1:53063/1‑251 (MQ=21)                                                                                                      aataatGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGttt                                                                                                                                             <  2:88854/251‑1 (MQ=255)                                                                                                      aataatGCGGGCAACGGTGCTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGttt                                                                                                                                             >  1:130855/1‑251 (MQ=255)                                                                                                                     ggTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACa                                                                                                                              <  1:415713/251‑1 (MQ=255)                                                                                                                                    cGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCAc                                                                                                               <  1:409318/251‑1 (MQ=255)                                                                                                                                    cGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCAc                                                                                                               <  1:397482/251‑1 (MQ=255)                                                                                                                                               aCAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTACGGAACACTTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGaaa                                                                                                     >  1:216077/1‑250 (MQ=255)                                                                                                                                                    cGGTGGCGGCGGCACCGTTGGTGCGGAGGCATCAGTCCGGAACAATTTGAAAACAAGAAACTCGATTAGGACTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGGTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGt                                                                                               >  2:50707/1‑251 (MQ=255)                                                                                                                                                      gTGGCGGCGGCACAGTTGGTGTGGCGGCCACAGACCGGAACAATTTGAAAATAAGAACCTCCCTTAGGCCTGTGTCCCTATTCCGTGGGTAGGATAAGAC‑CTCCCCTAGCGGGTCCTTTTGCATGCCGCCTGCGCTATTCTGTTTTTTTGCGCGCTGGCGAGCTGAAATCGGTCTTTAAACCCGACTTCCTGGGTC‑TCGGGTTCACTATTACCAGACATTATTTTTGCCGCACCCCCATGAAAGCGCGTa                                                                                              >  2:316891/1‑250 (MQ=21)                                                                                                                                                                    gTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGGGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTGGACTGCCGCCTGCGCTATTCTGTTGTTTGGCGCGCGGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGTTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGt                                                                               <  1:228046/251‑1 (MQ=255)                                                                                                                                                                                gcCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAgcgc                                                                   >  1:81476/1‑251 (MQ=255)                                                                                                                                                                                    cAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAAt                                                               >  1:257296/1‑251 (MQ=255)                                                                                                                                                                                        ccGGAACACTTCGAAAACAAGAACCTCGCTTCCGCCAGTCACCATATTCCGCGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAACCCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGc                                                           >  2:455455/1‑251 (MQ=255)                                                                                                                                                                                                   tGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCGGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTGCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTc                                                 <  1:168282/250‑1 (MQ=255)                                                                                                                                                                                                       aaCAAGACCCCCGCTTAGGCCGGTGTCCATATTACGGGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAACCGCTCTTTAACCCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTc                                            <  2:130741/251‑1 (MQ=255)                                                                                                                                                                                                        aCAAGACCCGCGTTTAGGCCTGTGTCCATATTACGGGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTAAAACCCGACTTCCGGGGTC‑CCGGTTTCACTATTACCAGACATTATTTGCGCCGCTCCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAAGCGGCACGGATTGTCGCTTcc                                           <  2:216077/251‑1 (MQ=255)                                                                                                                                                                                                        aCAAGAACCTCGATTAGGCCTGTGTCCATATTACGTGGGTAGCATCCGTC‑ATCCCTTCGTGAATCATTTTGACAGACGCCTGAGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAATTCGCTCTTTAAACCCGACTTCCTGGGGC‑TCGGATTCACTATTACAAGACATTATTTTCGCCGAACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGGAAGCGCAAATCGCCCCGTATTGTCGCCTcc                                           >  1:67890/1‑251 (MQ=255)                                                                                                                                                                                                           aaacccTCGCTGTGGCCTGTTTCGATATTACGGGGGTAGGATCAGGC‑TTCCCTTCGTGGATCATTTTGACTGCCGCCTGGGCTATTCTGTTTTTTGGCGCGCGGGCGAGCTGAAATCGTTCTTTAAACCCGACTTCC‑GGGTGGTCGGTTGCACTATGCCCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTTTCGCTTCCGAg                                        <  2:269284/247‑1 (MQ=255)                                                                                                                                                                                                           aGAACCTCGCTTAGGCTGGTGTCCATATTACGTGGGTAGGATCAGTC‑ACCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCTCGCTGGCGAGCTGAAATCTCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAg                                        <  2:22877/251‑1 (MQ=255)                                                                                                                                                                                                                                tCCATATTACGTGGGTAGGATC‑GTCCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCGGGCGAGCTGAAATCGCTCTTTAACCCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTCCGCCGCCCCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAg                    <  2:364469/250‑1 (MQ=255)                                                                                                                                                                                                                                     atTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTg              >  1:18550/1‑251 (MQ=255)                                                                                                                                                                                                                                                 ggATCAGTC‑ATCCCTTCGTGGATCATTTTTACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTCAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAACCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACAAGCACGGt  >  1:96824/1‑251 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                        ATGAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTGGGTAGGATCAGTC‑ATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCTGAAATCGCTCTTTAAACCCGACTTCCTGGGTC‑TCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACAAGCACGGT  >  NC_000913/2171183‑2171673

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.