Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	564,107	IS2 (+) +5 bp	coding (370‑374/633 nt)	sfmZ ←	response regulator family protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	381260 =	NA (NA)	35 (1.060)	35/598	0.0	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 564111	0 (0.000)					coding (370/633 nt)	sfmZ	response regulator family protein
*	?	NC_000913	= 382590	NA (NA)	54 (1.640)	48/600	0.0	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	564107 =	0 (0.000)					coding (374/633 nt)	sfmZ	response regulator family protein

GCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATACCAGGAAGGAAGTCTTACTGCTGTCGCCGCCGGAGAACAGG  >  NC_000913/380982‑381549                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  gCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAATCCGCAATAAAACATTGCGGACCAATCGTTGCATCCGCTTCAATAACCGCATTTTTAATCTCCGTGGCAAAACCAATTTTTACGCCATTGCTGATTCTTTTGCCCGGACGAATAAACGCATAATTACTATAATGGCCCCCTGAAtctcta                                                                                                                                                                                                                                                                             >  1:280653/1‑299 (MQ=255)  cGTGACCGCTGACCGATATAACAACCTAATTTATCGCATCCTGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCa                                                                                                                                                                                                                                                                             <  1:262936/300‑1 (MQ=255)          cTGCCCGATATAACACCCTAATTTATCGCATCCGGTAGCGATCATTCCCTCTGGAGTTCGAACAGACACGGACTGTTAATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCaa                                                                                                                                                                                                                                                                    <  2:31002/301‑1 (MQ=255)            gCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAata                                                                                                                                                                                                                                                                  >  1:40841/1‑301 (MQ=255)              ccGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAatat                                                                                                                                                                                                                                                                 >  1:128201/1‑300 (MQ=255)                      actacCTAAATTATCGTATCCGGTCGCGATAATGCCCTTTGGATGTCGAACAGACACGGGCTGCTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCCATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAAt                                                                                                                                                                                                                                                        <  2:148648/298‑1 (MQ=255)                       aaacATAATCTACCGCAGCCGGTAGCGATAGTTCCCTCTGGAGTTCGAACAGACACGTGCTGGTCATCCAGAAGAAGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATa                                                                                                                                                                                                                                                       <  2:20451/300‑1 (MQ=255)                               attATCGCTACCGGATGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAg                                                                                                                                                                                                                                               >  1:179529/2‑301 (MQ=255)                                                cGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATTTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTTGttttt                                                                                                                                                                                                                              >  2:285260/1‑301 (MQ=255)                                                  ataaataacgcaggAGTACGAACAGAGACGGGCTGTTCCATCATACGGTGATTACTAGTACGTGCTAGTGCACCTAAATATGACTCGTTTGCAACTACCAAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttttt                                                                                                                                                                                                                              <  2:139679/287‑1 (MQ=38)                                                  aTAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAg                                                                                                                                                                                                                            <  1:285272/301‑1 (MQ=255)                                                              tgattgtgaaatgacatatggtggtCATCCATACGATGATTACTCATACGGACTTGGGCGCCCAAGTAACACTGATTAGCAACTACCGAGGCCCGAAGAAAGCAATGCGGAGCAACCGTGGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTAATGTGTCAGGACGAATAAACGCCTAATTACTAGACTGTCCCCCTGAATCTCCAGACAACCAGTATCACTTAAATAAGGGATAGTCTTAATACTAGTTTTTAGACTAGTCATTg                                                                                                                                                                                                                 <  2:184630/276‑1 (MQ=21)                                                                  ttCGAACAGACACGGGCTGTTCATACAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGgaga                                                                                                                                                                                                             <  2:128192/300‑1 (MQ=255)                                                                                          ccAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAACCCAATTTTTACGCCATTGCTGATTATTTTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCCACAGCCAACCAATATCACTTAAATAAGTGATAGTCTTACTACTCGATTTTCGCCTCGTCATTTGATAACGTAtgattacttccttaggctc                                                                                                                                                                                     >  2:148799/1‑286 (MQ=255)                                                                                                                tacTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACca                                                                                                                                                              >  1:194374/1‑301 (MQ=255)                                                                                                                tacTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACAAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGATTAGTCATTGGCGAAACGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACca                                                                                                                                                              >  1:171738/1‑301 (MQ=255)                                                                                                                         ccctaagttgtccTGGTTTGACAATACCGAGTCTGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGaaa                                                                                                                                                     <  2:40838/290‑1 (MQ=255)                                                                                                                                                   cctggtcgacaagaagtgaTTGGGGACAAATCCTGGCTTACGCTTCAATAACCGCAATTTTAATCTCGGTGGCGAAACCAACTCGTACGCCATTGCTGATTATTTCGCCAAGACGAATAAACGCATAATTACTAGTCTGGGGCCCTGAATCTCCAGACAACAAATATCACTTAAAGAAGTGATAGTCTTAATACTAGCTTTTAGACTCGTCATTGGAGACCAGATGATTGTTGTCTTAGGGCCGGACAAAAGCGGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGGGCttt                                                                                                                           <  2:275520/283‑1 (MQ=18)                                                                                                                                                   ccGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCAAGGACGAAGAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTATTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACTACACAGGAAAAGATCGCAATTGTTCAGCAGAGCttt                                                                                                                           >  1:74426/1‑301 (MQ=255)                                                                                                                                                          ggCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCCGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGACTCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGACTAACAGATGATTGATGTCTTAGGGCCGTAGAAACGCAGACGGCGTACCACACAGGCAAAGATCGCACTTGTTCAGCAGAGCTTTGAACCgg                                                                                                                    >  1:225272/1‑301 (MQ=255)                                                                                                                                                            cAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAACACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAATAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGCTGTCTTAGGGCCGGAACAACGCAGACGGCGTACCACCCAGGAAACTATCGCAATTGTTctgaagctgcttcacacggtg                                                                                                                  >  2:68598/1‑281 (MQ=25)                                                                                                                                                            cAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACAACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCgggg                                                                                                                  >  1:262642/1‑301 (MQ=255)                                                                                                                                                                                                       cGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGATTTTAGACTAGTCATTGGAGAACAGATGCTTGATGGCTTAGGGCCGGAGAAACGCAGACGGCGTACCAAACAGGAAAAGATCGAAATTGTTCATCAGGACTTTGAACCGGGGATGACGGTCTCCCTCGTTGCTCGCAAACAGGGCGGatcagcaa                                                                       >  2:217785/1‑299 (MQ=18)                                                                                                                                                                                                                    ccGGTGGCAAAACCAATTTTTACGCCATTGTTTGTTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATTCTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTtct                                                          <  2:179519/300‑1 (MQ=255)                                                                                                                                                                                                                       ttGGCAAAACCAATGTTTACACCATATCTCATTATTGTGCCAGTACGAATAAACGCATAATTACCAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGGGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATAGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTg                                                       <  2:194364/300‑1 (MQ=255)                                                                                                                                                                                                                            aaaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGACCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGtt                                                  >  2:39451/1‑300 (MQ=255)                                                                                                                                                                                                                                 cAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCaa                                             >  1:142736/1‑301 (MQ=255)                                                                                                                                                                                                                                     ttttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATAcc                                         <  2:142727/301‑1 (MQ=255)                                                                                                                                                                                                                                      tttACGCCATTGCTCGATATTGTGCCAGGACGAAAAAACGAATAAGTACAAGACTGGCCCCCCGACTCTCCAGACAACCAAAATCACATAAATAAGTGATAGTCTTAATACGTGTTTTTAGACTAGTCAATGGAGAACAGATGATTGATGTCTTACGGCCGGAGAACCGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATACCa                                        <  2:262631/301‑1 (MQ=14)                                                                                                                                                                                                                                                     cgccgctcttccgatctaaTAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACGCAGGAAAAGATCGCAATTGTTCAGCAGGGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATACCAGGAAGGAAGTCTTAc                         <  2:83085/284‑1 (MQ=255)                                                                                                                                                                                                                                                       tattGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATACCAGGAAGGAAGTCTTActg                       <  1:39454/301‑1 (MQ=40)                                                                                                                                                                                                                                                        tttGAGCCAGGAAGAATAAACGCATAATCACTAGGCTGGCCCCCTGAATCTCCAGACAACAAATATCACTTAAATAAGTGATAGTCGTAATACTAGTTTTTAGAATCGTCATTGGAGAACAGAGGATTGAGGTCCTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATACCAGGAAGGAAGTCTTActgc                      <  2:225262/300‑1 (MQ=11)                                                                                                                                                                                                                                                          tgtgCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATACCAGGAAGGAAGTCTTActgctg                    <  1:217795/301‑1 (MQ=39)                                                                                                                                                                                                                                                                    cGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTAAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGTTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGNAATACCAGGAAGGAAGTCTTACTTCTGTCGCCGCCgg          >  2:49986/1‑301 (MQ=25)                                                                                                                                                                                                                                                                      aaTAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTAGTTGCTCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATACCAGGAAGGAAGTCTTAcagatcggaagagcacaa        >  1:83091/1‑284 (MQ=255)                                                                                                                                                                                                                                                                            cGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGAGCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCACGGCAACATGGTGTAGCAGCCAGCCAGTTATTCCTCTGGCGTACGCAATACCCGGAAGGAAGTCGTACTGCTGTCGCCGCCGGAGAACAgg  >  2:144444/1‑301 (MQ=12)                                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                                                                  GCGTGAACGCTGCCCGATATAACAACCTAATTTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAATACCAGGAAGGAAGTCTTACTGCTGTCGCCGCCGGAGAACAGG  >  NC_000913/380982‑381549

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