Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	572,857	IS2 (+) +5 bp	coding (392‑396/456 nt)	ybcN →	DNA base‑flipping protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 382590	NA (NA)	58 (1.290)	56/490	0.1	96.6%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	572857 =	2 (0.050)					coding (392/456 nt)	ybcN	DNA base‑flipping protein
*	?	NC_000913	= 572861	2 (0.050)	44 (0.980)	42/490	0.2	95.6%	coding (396/456 nt)	ybcN	DNA base‑flipping protein
	?	NC_000913	= 1469240	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

CCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGG  >  NC_000913/382352‑382840                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                           ccGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAg                                                                                                                                                                                                                                                <  2:390643/251‑1 (MQ=33)                      gagagCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTAc                                                                                                                                                                                                                           >  1:172448/1‑251 (MQ=255)                                       gATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACgg                                                                                                                                                                                                          >  1:294434/1‑251 (MQ=255)                                        aTAAAGCGTGACAACATCAGAATCATGAACAAACCTGACGGGCTCACGACCGACACGAACCTTGCAAAGGCGTTCGAGCATTCTAACGAATGGCATACGCATCGAGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGGTGGCTTGTAATGGGTTAAGTGCTAACCGATGTCTGGAAATATAGGGGCAAATCAAATTACCTATCAGGCAGGTTGCGCCAATCACTACCGGACCCTGAACAACGGc                                                                                                                                                                                                         >  2:50539/1‑251 (MQ=17)                                              cGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTa                                                                                                                                                                                                                            <  2:258208/226‑1 (MQ=255)                                              cGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTa                                                                                                                                                                                                                            >  1:258208/1‑226 (MQ=255)                                                                aTGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACATATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTcgc                                                                                                                                                                                 <  2:326868/251‑1 (MQ=255)                                                                  gCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCgcgc                                                                                                                                                                               >  1:209351/1‑251 (MQ=255)                                                                                  ttAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAcc                                                                                                                                                                                                      <  1:401767/212‑1 (MQ=255)                                                                                  ttAACGACAGCAAAGAAACTTGCCGAGGCGTTCGAGCATTATAACGAATGGCATCCGAATAGTGCGCTGGGTTATCGCTCGCCACGGGAATCTCTGCGGCAGCGGGCTTGTAATGGGGTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAcc                                                                                                                                                                                                      >  2:401767/1‑212 (MQ=18)                                                                                    aaCGGCAGCAAAGAACCTTGCCGATGCGATCGATCATTATAACGAATGGCATCCGCATAGTGCTCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTCCCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCtt                                                                                                                                                             <  2:362809/251‑2 (MQ=255)                                                                                    aaCGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTg                                                                                                                                                             <  1:154349/251‑1 (MQ=255)                                                                                                gAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCaa                                                                                                                                                 <  2:270596/251‑1 (MQ=255)                                                                                                      ggCAGAGGCGTGCGAGCATTATAACGAATGGCATCCTCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCCCCATGACAAATACGGGTATTCGCGCCGATAACAACTTCCCCTGCAGTTTCATCAATAAtt                                                                                                                                            <  1:139254/249‑1 (MQ=255)                                                                                                           agGCGTTCGAGTATTATAACGAATTGCAGCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATTCGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAAGTACCTgg                                                                                                                                      <  2:289752/251‑1 (MQ=255)                                                                                                                  cGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGTTCTGCaa                                                                                                                               <  1:375807/251‑1 (MQ=255)                                                                                                                    aGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATc                                                                                                                             >  1:190914/1‑251 (MQ=255)                                                                                                                         tataACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAAGTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGaa                                                                                                                         <  1:241300/250‑1 (MQ=255)                                                                                                                               gAATGGCATCCGCATAGTGCGCTGGGTTATCTCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATAttttt                                                                                                                  <  2:109426/251‑1 (MQ=255)                                                                                                                               gAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGGCTGCAATCTGAATAttttt                                                                                                                  <  1:376704/251‑1 (MQ=255)                                                                                                                                 aTGGCATCCGCATAGGGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCGGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGTTATTCTCGCCGATAACAACTTCACCTGCAGTTTCATCAAGAATTACCTGGTCTGCAAGCTGAATATTTTtct                                                                                                                <  2:411878/251‑1 (MQ=255)                                                                                                                                      aTCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCttt                                                                                                            >  2:366885/1‑250 (MQ=255)                                                                                                                                           cATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTa                                                                                                                                           <  1:19785/214‑1 (MQ=255)                                                                                                                                           cATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTa                                                                                                                                           >  2:19785/1‑214 (MQ=255)                                                                                                                                             tAGTGCGCTGGGTTCTCGATCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGCCCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTa                                                                                                     >  2:189498/1‑250 (MQ=255)                                                                                                                                                 gcgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCt                                                                                                 <  1:1857/250‑1 (MQ=255)                                                                                                                                                 gcgcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCAACAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCt                                                                                                 <  1:315117/250‑1 (MQ=255)                                                                                                                                                   gcTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTCATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTc                                                                                               >  2:179085/1‑250 (MQ=255)                                                                                                                                                      gggTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTc                                                                                               >  1:178215/1‑247 (MQ=255)                                                                                                                                                        gTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAg                                                                                         <  1:366885/251‑1 (MQ=255)                                                                                                                                                            tctctcgacccaggagactCTGCAGCAGCGGGCTTGTAATGGGTTTCGTGCTAACAGCTGTCTGGAACTATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGCCTGCCATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGtttt                                                                                     >  2:69984/19‑251 (MQ=255)                                                                                                                                                             cgctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGtttt                                                                                     <  2:190914/250‑1 (MQ=255)                                                                                                                                                              gctcgcCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTACGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTGTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGAAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTACTCTTCAGATAGTTTTTc                                                                                   >  2:370700/1‑251 (MQ=255)                                                                                                                                                                       gggAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATACTTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCtgt                                                                          >  2:367208/1‑251 (MQ=255)                                                                                                                                                                                     cAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGttt                                                                                      >  2:116778/1‑225 (MQ=255)                                                                                                                                                                                     cAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGttt                                                                                      <  1:116778/225‑1 (MQ=255)                                                                                                                                                                                          ggCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATaa                                                       <  2:172448/251‑1 (MQ=255)                                                                                                                                                                                              tGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCgg                                                   >  1:221093/1‑251 (MQ=255)                                                                                                                                                                                                          aGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTcc                                        >  2:345274/1‑250 (MQ=255)                                                                                                                                                                                                                   aGATGTCTGGAAATATAGGGGCCAATCCAATTACCTATCACGCAGTTTGCGCCAATCACTACCGGCACCTGAATAACGGCACCATGTCAAATACGGGCCTTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTCACTGGTCTGCAATCTGAATATTTTTCTCTTTAAAATTAAGCTTCAGATAGTTTTTCAGTTCTTGTGTCAGAAGGTTGGAAATAAGCGGTTGACGCTTCAACCAGTGCgg                              >  2:120721/1‑251 (MQ=255)                                                                                                                                                                                                                                 ataGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtgggt                >  2:397620/1‑251 (MQ=255)                                                                                                                                                                                                                                      ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTc           <  2:343491/251‑1 (MQ=255)                                                                                                                                                                                                                                      ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTc           >  2:443267/1‑251 (MQ=255)                                                                                                                                                                                                                                           aTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGTCAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCGGTAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTGTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCt      <  1:261642/251‑1 (MQ=255)                                                                                                                                                                                                                                              cAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAg   <  2:178215/251‑1 (MQ=255)                                                                                                                                                                                                                                               aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAgg  <  2:294434/251‑1 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                           CCGGAGAGTAACGGAATAGCAGAGAGCTTCGTGAAAACGATAAAGCGTGACTACATCAGTATCATGCCCAAACCAGACGGGTTAACGGCAGCAAAGAACCTTGCAGAGGCGTTCGAGCATTATAACGAATGGCATCCGCATAGTGCGCTGGGTTATCGCTCGCCACGGGAATATCTGCGGCAGCGGGCTTGTAATGGGTTAAGTGATAACAGATGTCTGGAAATATAGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGG  >  NC_000913/382352‑382840

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