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)	48 (0.740) +TCA	44/482	0.5	98.0%	noncoding (1/1255 nt)	IS3	repeat region
	?	NC_000913	3720136 =	1 (0.020)					coding (88/213 nt)	cspA	cold shock protein CspA
*	?	NC_000913	= 2171428	NA (NA)	53 (0.810)	49/486	0.4	100%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	= 3720138	0 (0.000)					coding (90/213 nt)	cspA	cold shock protein CspA

GAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACA  >  NC_000913/2171185‑2171665                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               gAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATcc                                                                                                                                                                                                                                         <  2:362170/251‑1 (MQ=32)        aaaGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAAAAAGAACCTCGCTTAGGCCTGGGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTg                                                                                                                                                                                                                                  >  2:5855/1‑251 (MQ=25)               gctgctACGATAATGCCTGCGTGGAAAGCTTCTTTGAGTCGCTGAAAGTGGAATGTATCCATGGAGAACTCTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATGACAATCGGTGGCGGCGGCACAGTTGGTGTGGCTGCCTCAGTCCGGAACTATGTGAAAACAAGAACCTCTCTTAGGCCTGTTTCCATATTACGAG‑GGTAGGATCTGTCATCCCTTCGTGGATCAtt                                                                                                                                                                                                                           <  2:311264/251‑1 (MQ=11)               gctgctACGATAATGCCTGCGGGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACGTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGTTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCATTCATCCCTTCGTGGATCAt                                                                                                                                                                                                                            <  2:615119/250‑1 (MQ=255)                  gctACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCAttt                                                                                                                                                                                                                          <  1:5855/249‑1 (MQ=40)                         tAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGGCGGAAAATAATGCGGGCAACGGTGTTTAATTATATCGAGTGTAATTACAATCGGTTGCGGCGGCACAGTTGGTGTGGCGTCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGt                                                                                                                                                                                                                                              <  1:423537/222‑1 (MQ=11)                                ggCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTgcgc                                                                                                                                                                                                          <  2:16462/250‑1 (MQ=255)                                     ggAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTc                                                                                                                                                                                                     >  1:524802/1‑251 (MQ=255)                                      gAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCt                                                                                                                                                                                                    <  1:487738/251‑1 (MQ=255)                                         aGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGtt                                                                                                                                                                                                 <  1:29944/251‑1 (MQ=255)                                                          gAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGGTCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCgg                                                                                                                                                                                >  2:490442/1‑250 (MQ=255)                                                               gggAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑Ga                                                                                                                                                                          <  2:596422/250‑1 (MQ=255)                                                                ggAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑Ga                                                                                                                                                                          <  2:210835/250‑1 (MQ=255)                                                                   aTGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCg                                                                                                                                                                                 >  1:263752/1‑241 (MQ=255)                                                                   aTGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCg                                                                                                                                                                                 <  2:263752/241‑1 (MQ=255)                                                                          cATGGATAACACTTTGTCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATACCGAATGTGATTCCAATCGGTGGCGGCTGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAATAACCTCTCTTGTGCCTGTGTCCATATTAC‑TGTGGTAGGATCTGTCATCCCTTCGTTGATCATTTTGACTGCCGCCTGCGCTATTCTATTTTTTGGCGCTCTATCGAGCT‑GAAATCGCTCttt                                                                                                                                                               <  2:427284/251‑1 (MQ=17)                                                                              gagaACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAAAAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAAc                                                                                                                                                           >  2:339046/1‑251 (MQ=255)                                                                                aatcacTGTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAAccc                                                                                                                                                         <  1:162794/248‑1 (MQ=255)                                                                                 aaCACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCg                                                                                                                                                        <  1:243485/251‑1 (MQ=255)                                                                                    acTTTATCAGACGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACt                                                                                                                                                     >  2:612434/1‑251 (MQ=255)                                                                                        tATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGAGTGGCGGCCTCAGACCGGTATAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTTCGGG‑GGTAGGATCAGTCATCCCTTCGTGGTTCATTTTGACTGCCGACTGCGCTTTTCTGGTTTTTGGTGCGCCGGAG‑GATAGAAATCGCTCTTTTAACCCGCCTTCCt                                                                                                                                                 >  2:392627/1‑251 (MQ=17)                                                                                           cAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGGGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGAGCATTTTGACTGCCGCCTGCGCTAGTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCCCTCGTTAAACCCGACTTCCTgg                                                                                                                                               <  2:95636/250‑1 (MQ=255)                                                                                                    aataatGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTTAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCCTATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCCTTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGttt                                                                                                                                     >  2:458470/1‑251 (MQ=255)                                                                                                    aataatGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTATTTTGGCGCGCTGGCGAGCT‑GAAAACGCTCTTTAAACCCGTCTTCCTGGGTCTCGGttt                                                                                                                                     <  2:495024/251‑1 (MQ=255)                                                                                                       aatGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAATCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTccc                                                                                                                                  >  2:233923/1‑249 (MQ=255)                                                                                                        atGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCGCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTGGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGc                                                                                                                                                                                  <  1:117048/203‑1 (MQ=255)                                                                                                        atGCGGGCAACGGGGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGACCCCGTCCGGAACCATTTGAAAACAAGCACCTCGCCTAGGCCTGTGTCCATATTACGTG‑GGTAGGCTCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGc                                                                                                                                                                                  >  2:117048/1‑203 (MQ=18)                                                                                                            gggCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTAtt                                                                                                                              >  1:571100/1‑250 (MQ=255)                                                                                                             ggCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCGGGGTCTCGGTTTCACTATTAc                                                                                                                            <  1:541849/251‑1 (MQ=255)                                                                                                                  cGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCGCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGAc                                                                                                                       <  2:181560/251‑1 (MQ=255)                                                                                                                      gtTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGTCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTTCGCGATTCTTTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTGACTATTACCAGACatt                                                                                                                    <  1:482821/250‑1 (MQ=255)                                                                                                                        ttAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACattat                                                                                                                  <  2:526473/250‑1 (MQ=255)                                                                                                                         tAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAAcc                                                                                                                                                          >  1:218416/1‑209 (MQ=255)                                                                                                                         tAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAAcc                                                                                                                                                          <  2:218416/209‑1 (MQ=255)                                                                                                                          aaTTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGGAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGGGCGCTGGCGAGCT‑TAAATCGCTCTTTAAACCCGGCTTCCGGGGTCTCGGTTTCACTATTTCCAGACATTAttt                                                                                                                <  2:691630/250‑1 (MQ=255)                                                                                                                           aTTATATCGAATGTTAGTACAAGCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCTGAACAATTTGAAAACAATAACCTCGCTTAGGCCTGCGTCCATATTACGTG‑GGTAGTATCAGTCAGCCCGTCGTGGATCATTTTTACTGCCGCCTGCGCTAGTCTGTGTTTTGTCGCGCTGTCGAGCG‑GAAATCGctct                                                                                                                                                                 <  2:625709/200‑1 (MQ=17)                                                                                                                           aTTATATCGAATGTGATTACACTCGGTGGCGGCGGCAGAGTTGGTGTGGCGGCCTCAGTCCGGAAAAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCCTATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTTACTGACGCATGCGCTCTTCTGTTTTTTGGCGCGCTGGCGGGCT‑GAAATCGctct                                                                                                                                                                 >  1:625709/1‑200 (MQ=18)                                                                                                                                       gtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCa                                                                                                  <  1:670446/251‑1 (MQ=255)                                                                                                                                       gtgATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGGCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGcc                                                                                                   <  2:430126/250‑1 (MQ=255)                                                                                                                                            tACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGaa                                                                                              >  1:476211/1‑250 (MQ=255)                                                                                                                                               aaTCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCCGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGGTCATTTTGACTGCCGCCTGCGCTCTTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAAACCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAgcg                                                                                          >  1:650861/1‑251 (MQ=255)                                                                                                                                                    gTGGCGGCGGCACAGTTGGGGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGGAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGCCTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTa                                                                                      >  1:87728/1‑250 (MQ=255)                                                                                                                                                          gcggcACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGc                                                                               <  1:288001/251‑1 (MQ=255)                                                                                                                                                                                                       cAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCg                                  >  1:153169/1‑251 (MQ=255)                                                                                                                                                                                                                cGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGgcagca                          >  1:481825/1‑250 (MQ=255)                                                                                                                                                                                                                     aGGCCTGTGTCCATATTCCGGG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAACCCCGACTTCCGGGGTCTCGGTTTCACTATTCCCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATAc                     <  1:523356/250‑1 (MQ=255)                                                                                                                                                                                                                             gtCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCAGAGGCAGCAACTACTCCCGcc              >  2:1751/1‑249 (MQ=255)                                                                                                                                                                                                                               ccATTTTACGTG‑GGTAGTACCAGTCCTCCTTTCTTGAATCATTTTGCCGGCCGCCTGCGCTATTCTGTTTTTGGGCCCGCGGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTGCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGAt          <  2:524802/251‑1 (MQ=255)                                                                                                                                                                                                                                   atTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTg      <  1:458470/251‑1 (MQ=255)                                                                                                                                                                                                                                      acggg‑ggTAGGATCAGTCACCCCTTCGTGGATCATTTGGCCGGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAACCCCGACTTCCGGGGTCTCGGTTTCACTATAACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCg                <  2:540272/234‑1 (MQ=255)                                                                                                                                                                                                                                      aCGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCg                >  1:540272/1‑238 (MQ=255)                                                                                                                                                                                                                                         tG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCATGAACa  >  1:594025/1‑249 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               GAGCGCAAAAGGTTGCTGCTACGATAATGCCTGCGTGGAAAGCTTCTTTCATTCGCTGAAAGTGGAATGTATCCATGGAGAACACTTTATCAGCCGGGAAATAATGCGGGCAACGGTGTTTAATTATATCGAATGTGATTACAATCGGTGGCGGCGGCACAGTTGGTGTGGCGGCCTCAGTCCGGAACAATTTGAAAACAAGAACCTCGCTTAGGCCTGTGTCCATATTACGTG‑GGTAGGATCAGTCATCCCTTCGTGGATCATTTTGACTGCCGCCTGCGCTATTCTGTTTTTTGGCGCGCTGGCGAGCT‑GAAATCGCTCTTTAAACCCGACTTCCTGGGTCTCGGTTTCACTATTACCAGACATTATTTTCGCCGCACCGCCATGAAAGCGCGTAACCACGCCTTTTTGTTCGAGAAAGCGCAAATCGGCACGGATTGTCGCTTCCGAGGCAGCAAATACTCCCGCCAGATCCTGAACA  >  NC_000913/2171185‑2171665

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