Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	1,333,262	Δ3 bp	coding (2215‑2217/2598 nt)	topA →	DNA topoisomerase 1

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1333261	2 (0.050)	20 (0.510)	20/482	0.9	90.9%	coding (2214/2598 nt)	topA	DNA topoisomerase 1
	?	NC_000913	1333265 =	2 (0.050)					coding (2218/2598 nt)	topA	DNA topoisomerase 1

AAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTAAGTTGCATGTCTGTGGTAATAACCCAACCTGCGACGGTTACGAGATCGAAGAGGGCGAATTCCGCATTAAAGGTTATGACGGCCCGATCGTTGAGTGTGAAAAATGTGGCTCTGAAATGCACCTGAAAATGGGGCGTTTCGGTAAATACATGGCCTGCACCAACGAAGAGTGTAAAAACACACGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCCTGGCTGCCAACACTTTCCCGAAATCGCGTGAAACGCGTGCGCCACTGGTGGAAGAGCTTTATCGCTTCCGCGACCGTCTGCCGGAAAAACTGCGTTATCTGGCCGATGCGCCAC  >  NC_000913/1333022‑1333502                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               gcaaaacgtcgttgcccgaaatgcggcacgtcgatggacagctatctcatcgatccgaaacgtaagttgcatgtctgtggtaataacccaacctgcgacggttacgagatcgaagagggcgaattccgcattaaaggttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGa                                                                                                                                                                                                                                          >  1:399875‑M2/244‑250 (MQ=255)                                     gacagctatctcatcgatccgaaacgtaagttgcatgtctgtggtaataacccaacctgcgacggttacgagatcgaagagggcgaattccgcattaaaggttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAagaaga                                                                                                                                                                                                    <  1:327442‑M2/44‑1 (MQ=255)                                           tatctcatcgatccgaaacgtaagttgcatgtctgtggtaataacccaacctgcgacggttacgagatcgaagagggcgaattccgcattaaaggttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGACAGAAGATCcggg                                                                                                                                                                                              >  1:198714‑M2/202‑250 (MQ=255)                                                     atcataaacgtatgttgcattacagtggtaataacctaatcttcgacggttacgagaacgaagatggcgaatcccgcattaaaggttatgacgtcccgatcgttgagtgtgaaaaatgtggctctgaaattcacctgaaaatggggcgttgcgttaaatacatggccttcaccaccgaagagtgtaaaaacCGTAAGATTTTACGTACCGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCt                                                                                                                                                                                    <  2:327922‑M2/60‑1 (MQ=255)                                                       ccgaaacgtaagttgcatgtctgtggtaataacccaacctgcgacggttacgagatcgaagagggcgaattccgcattaaaggttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgacgagtgtaaaaacCGTAAGATTTTACGTACCGGCGAAGTGGCACCAC‑CGAAAGACGATACGGTGCCATTACCTg                                                                                                                                                                                   >  2:396814‑M2/190‑250 (MQ=255)                                                           aaagtaagtggcatgtcggtggtaataacccaacctgcgacggttacgagatcgaagagggcgaattccgcattaaaagttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTg                                                                                                                                                                              <  2:398762‑M2/66‑1 (MQ=255)                                                             cgtaagttgcatgtctgtggtaataacccaacctgcgacggttacgagatcgaagagggcgaattccgcattaaaggttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGcc                                                                                                                                                                            <  2:266546‑M2/68‑1 (MQ=255)                                                                             gtggtaataacccaacctgcgacggttacgagatcgaagagggcgaattccgcattaaaggttatgacggcccgatcgtggagggtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGAt                                                                                                                                                            <  2:399875‑M2/84‑1 (MQ=255)                                                                                          aacctgcgacggttacgagatcgaagaggtcgacttccgcattaaaggttatgacggcccggtcgttgagagagaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGCAGTGGCCCCAC‑CGAAAGAAGATCCAGTGCCATTACCCGCGCCACCGTGCGAAAAATCAGATGTCTatttctggc                                                                                                                                               >  2:298198‑M2/155‑247 (MQ=255)                                                                                              ggcgaaggtgacgagatcgaagagggcgaattccgcattaaaggttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCg                                                                                                                                           <  2:198714‑M2/101‑1 (MQ=255)                                                                                                       tacgagatcgaagagggcgaattccgcattaaaggttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGtg                                                                                                                                   <  1:251707‑M2/109‑1 (MQ=255)                                                                                                                      ggcgaattccgcattaaaggttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaactacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTCACGGCGAAGTTGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGAGCGACAAATCAGATGCTTATTTCGGGCTGCGTGCCgatgctggcggggggatact                                                                                                                   >  2:388334‑M2/127‑242 (MQ=255)                                                                                                                           attccgcattaaaggttatgacggcccgatcgttgagtgtgaaaaatgtggctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGAGTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAAGCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCCTGGCTg                                                                                                              <  1:382377‑M2/130‑1 (MQ=255)                                                                                                                                                 cctgtCGTTGAGTTTTAAAAAAGGGGGTCTGAAAGGCACCTGAAAAGGGGGCATTTCGGGAAAAACATGGCCTGCACCAACGAAGAGGGTAAAAAC‑‑‑CGTAAGATTTTACGTAAGGGCGAAGGGGCACCAC‑CGAAAGAAGATCCGGTGCCATTGCCTGAGGTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCATGGCTGCCAACTCGTTCCCGAAATCGCGTGa                                                                                     <  2:313974/247‑1 (MQ=255)                                                                                                                                                                tgtgaaaaatgtggctctgaaatgcacctgaaaatggggagtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCCTGGCTGCCAACACTTTCCCGAAATCGCGTGAAACGCGTGCGcc                                                                         >  1:159250‑M2/85‑251 (MQ=255)                                                                                                                                                                             gctctgaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCCTGGCTGCCAACACTTTCCCGAAATCGCGTGAAACGCGTGCGCCACTGGTGGAagag                                                            <  1:323235‑M2/180‑1 (MQ=255)                                                                                                                                                                                  gaaatgcacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCCTGGCTGCCAACACTTTCCCGAAATCGCGTGAAACGCGTGCGCCACTGGTGGAAGAGCttt                                                        <  1:396814‑M2/184‑1 (MQ=255)                                                                                                                                                                                        cacctgaaaatggggcgtttcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTACCGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCGGCCGGTGTGTTCCTGGCTGCCAACACTTTCCCGAAATCGCGTGAAACGCGTGCGCCACTGGTGGAAGAGCTTTATCGCtt                                                 <  1:207405‑M2/191‑1 (MQ=255)                                                                                                                                                                                                 atggggcgtttcggtaaacacatggcctgcaccaacgaagaagttaaaaacCGTAACATCTTACTTTACGGCGAAGTGGCCCCAC‑CAAAAGAAGATCCTGTGCCATTACCTGAGCGGCCGTGCGAAAAATCAGATGCTTATTGCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCCTGGCTGCCAACACTTTCCCGAAATCGCGTGAAACGCGTGCGCCACTGGTGGAAGATCTTTATCTCTTCCGCGACCg                                        <  2:156521‑M2/200‑1 (MQ=255)                                                                                                                                                                                                           tcggtaaatacatggcctgcaccaacgaagagtgtaaaaacCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCCTGGCTGCCAACACTTTCCCGAAATCGCGTGAAACGCGTGCGCTACTGGTGGAAGAGCTTTATCGCTTCCGCGACCGTCTGCCGGa                               <  1:15929‑M2/209‑1 (MQ=255)                                                                                                                                                                                                                 aCATGGCATGCAACAAAGAAGAGTGTAAAAACA‑‑‑GTAAGATTTTAAGTAACGGCGAAGTGGCAC‑ACACGAAAGAAGATCCGGTGCCATTACCGGAGCTGCCGTGCGAGAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGAGTTACTGGCTGCCAACTCTTTCACGAAATCGCGTGACACGCGGGCGCCACTGGTGGAAGAGCTTTCTCGCCTCCGCGACTGTCGGCCGGAAAAAcagagg                     >  2:82220/1‑246 (MQ=255)                                                                                                                                                                                                                                        gagtgtaaaaccCGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCCTGGCTGCCAACACTTTCCCGAAATCGCGTGAAACGCGTGCGCCACTGGTGGAAGAGCTTTATCGCTTCCGCGACCGTCTGCCGGAAAAACTGCGTTATCTGGCCGATGCGCcac  >  2:50691‑M2/13‑250 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                               AAACGTCGTTGCCCGAAATGCGGCACGGCGATGGACAGCTATCTCATCGATCCGAAACGTAAGTTGCATGTCTGTGGTAATAACCCAACCTGCGACGGTTACGAGATCGAAGAGGGCGAATTCCGCATTAAAGGTTATGACGGCCCGATCGTTGAGTGTGAAAAATGTGGCTCTGAAATGCACCTGAAAATGGGGCGTTTCGGTAAATACATGGCCTGCACCAACGAAGAGTGTAAAAACACACGTAAGATTTTACGTAACGGCGAAGTGGCACCAC‑CGAAAGAAGATCCGGTGCCATTACCTGAGCTGCCGTGCGAAAAATCAGATGCTTATTTCGTGCTGCGTGACGGTGCTGCCGGTGTGTTCCTGGCTGCCAACACTTTCCCGAAATCGCGTGAAACGCGTGCGCCACTGGTGGAAGAGCTTTATCGCTTCCGCGACCGTCTGCCGGAAAAACTGCGTTATCTGGCCGATGCGCCAC  >  NC_000913/1333022‑1333502

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