Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,974,440	IS30 (+) +2 bp	coding (833‑834/882 nt)	rffH →	glucose‑1‑phosphate thymidylyltransferase 2

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1469296 =	NA (NA)	73 (1.010)	68/492	0.2	100%	noncoding (1/1221 nt)	IS30	repeat region
	?	NC_000913	= 3974441	0 (0.000)					coding (834/882 nt)	rffH	glucose‑1‑phosphate thymidylyltransferase 2
*	?	NC_000913	3974440 =	0 (0.000)	99 (1.370)	86/492	0.1	100%	coding (833/882 nt)	rffH	glucose‑1‑phosphate thymidylyltransferase 2
	?	NC_000913	4507459 =	NA (NA)					noncoding (1221/1221 nt)	IS30	repeat region

CCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCACTATGTTAAGGGATACTGGCGGCATAAAACCCCATGAGCGTAAGCGGGCTGTAGCTCACCTGACACTG  >  NC_000913/1469056‑1469540                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     ccGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCa                                                                                                                                                                                                                                             >  1:530892/1‑250 (MQ=255)   gTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAAtt                                                                                                                                                                                                                                          <  2:43342/251‑1 (MQ=255)     cATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTg                                                                                                                                                                                                                                         <  1:589391/250‑1 (MQ=255)      aTCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGCCTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGt                                                                                                                                                                                                                                       >  1:510531/1‑251 (MQ=255)           cGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACg                                                                                                                                                                                                                                  <  1:182925/251‑1 (MQ=255)                   aGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGtt                                                                                                                                                                                                                          <  1:358555/251‑1 (MQ=255)                            gAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGaaaa                                                                                                                                                                                                                 <  1:632302/251‑1 (MQ=255)                              aCAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAAGGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTAGGTGAAAACa                                                                                                                                                                                                               <  1:565350/251‑1 (MQ=255)                               cAATTGCGATCGTTACCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTGAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCATCGCAACAGTTATGTGAAAACAt                                                                                                                                                                                                              <  1:265827/251‑1 (MQ=255)                                     cGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGtt                                                                                                                                                                                                        >  1:494468/1‑251 (MQ=255)                                         cctttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCACCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCg                                                                                                                                                                                                     <  2:169624/248‑1 (MQ=255)                                           tttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCg                                                                                                                                                                                                     <  2:275491/248‑1 (MQ=255)                                            ttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAg                                                                                                                                                                                                  <  1:454525/250‑1 (MQ=255)                                            ttCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAg                                                                                                                                                                                                  <  2:494468/250‑1 (MQ=255)                                               cTGTGTGGAACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGCCTAGTCTAAAAACTAGTCTTAAGACTATCACTTATTTAAGTGATATTGGTTGTATGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCACCAGTTATGTGAAAACATGGGGTTGCGGAGGttt                                                                                                                                                                                              >  2:383831/1‑251 (MQ=255)                                                                   gTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGaaaa                                                                                                                                                                                                                 <  2:157852/212‑1 (MQ=255)                                                                   gTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGaaaa                                                                                                                                                                                                                 >  1:157852/1‑212 (MQ=255)                                                                               cTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGa                                                                                                                                                               >  1:513490/1‑250 (MQ=255)                                                                                        caatcaTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCctct                                                                                                                                                     <  2:306717/251‑1 (MQ=255)                                                                                                   tctcCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAAc                                                                                                                                           <  2:525226/250‑1 (MQ=255)                                                                                                      ccAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGCCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATg                                                                                                                                       >  2:639178/1‑251 (MQ=255)                                                                                                        aaTGACTAGTCTAAAAACTAGTATTAAGACTATCACTGATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGga                                                                                                                                     <  2:530892/251‑1 (MQ=255)                                                                                                                    aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAAc                                                                                                                           >  1:340495/1‑249 (MQ=255)                                                                                                                       aaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCt                                                                                                                      >  1:571619/1‑251 (MQ=255)                                                                                                                               ttAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGa                                                                                                               >  2:93742/1‑250 (MQ=255)                                                                                                                                    ccTATCATTTATTTAAGTGATATTGGTTGTCGGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGTTGCAAAAATTGTTAATAAAGTTACTTGTGGATCCAATTGTTCAACGCAACGTTTAGGTGAAAACATGGGGTTGCG‑AGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTTTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCg                                                                                                        <  1:78628/250‑1 (MQ=21)                                                                                                                                     cTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCg                                                                                                        >  1:45629/1‑251 (MQ=255)                                                                                                                                       aTCACTTATTTAAGTGGTATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGa                                                                                                       >  1:200759/1‑250 (MQ=255)                                                                                                                                       aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCg                                                                                                        <  1:93742/249‑1 (MQ=255)                                                                                                                                                 tacgTGATATTGGTTGTCTGGAGATTCAGGTGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGTTGCGAAAATTGTTAATAAAGATACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGATACGCACCTTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAAAACGGAACAGGCTTCAGTGAAATCGCGAATATCCTGGGt                                                                                            >  2:415788/4‑251 (MQ=21)                                                                                                                                                       atatTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCaaa                                                                                       >  1:561132/1‑250 (MQ=255)                                                                                                                                                                               gggCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCACTATGTTa                                                              <  1:132136/251‑1 (MQ=255)                                                                                                                                                                                  ccAGTCTAAACCCTAGAACTGACTTATGCTGAGTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGACTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGTGACGTACAGTTACAGCAGAGGAAAAAGCTTCGGTGTATGAACCATGGAAGAACGGAACAGGCTTTATTGAAATAGCGAATATTATGGGTTCAAAAGCCGGAACGTTCGTCGCTAGGTTAAgg                                                           <  2:679650/251‑1 (MQ=17)                                                                                                                                                                                  ccAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCACTa                                                                   >  2:408017/1‑243 (MQ=255)                                                                                                                                                                                  ccAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCACTa                                                                   <  1:408017/243‑1 (MQ=255)                                                                                                                                                                                                        cttaTGCTGATTCCGGTGCGAAAATTGTTAATAAAGGAACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGCGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGCAAAAGCATCTGTTTTTGAACTATGGACGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCACTATGTTAAGGGATACTGGCGGCATAAAAcccc                                     >  1:47637/1‑251 (MQ=255)                                                                                                                                                                                                           aTGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCACTATGTTAAGGGATACTGGCGGCATAAAACCCCATg                                  >  2:536437/1‑251 (MQ=255)                                                                                                                                                                                                                       gTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAAGATCCTGGGTTCAAAACCCGGAACGAGCTTCACTATGTTAAGGGATACTGGCGGCATAAAACCCCATGAGCGTAAGCggg                      <  2:45629/251‑1 (MQ=255)                                                                                                                                                                                                                                     aataaAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCACTATGTTAAGGGATACTGGCGGCATAAAACCCCATGAGCGTAAGCGGGCTGTAGCTCACCTg        >  1:374789/1‑251 (MQ=34)                                                                                                                                                                                                                                           ggTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCACTATGTTAAGGGATACTGGCGGCATAAAACCCCATGAGCGTAAGCGGGCTGTAGCTCACCTGACACTg  >  1:340936/1‑251 (MQ=32)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                     CCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAAACCCCAGAACTTACTTATGCTGATTCCGGTGCGAAAATTGTTAATAAAGGTACTTGTAGATTCAATTGGTCAACGCAACAGTTATGTGAAAACATGGGGTTGCGGAGGTTTTTTGAATGAGACGAACATTTACAGCAGAGGAAAAAGCCTCTGTTTTTGAACTATGGAAGAACGGAACAGGCTTCAGTGAAATAGCGAATATCCTGGGTTCAAAACCCGGAACGATCTTCACTATGTTAAGGGATACTGGCGGCATAAAACCCCATGAGCGTAAGCGGGCTGTAGCTCACCTGACACTG  >  NC_000913/1469056‑1469540

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