Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,907,451	IS2 (+) +5 bp	coding (137‑141/144 nt)	yobF ←	DUF2527 domain‑containing protein YobF

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1651548	NA (NA)	37 (0.920)	36/480	0.2	100%	noncoding (706/706 nt)	IS2	repeat region
	?	NC_000913	1907451 =	0 (0.000)					coding (141/144 nt)	yobF	DUF2527 domain‑containing protein YobF
*	?	NC_000913	= 1907455	0 (0.000)	23 (0.570)	22/480	0.8	100%	coding (137/144 nt)	yobF	DUF2527 domain‑containing protein YobF
	?	NC_000913	= 2070271	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

TACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAAC  >  NC_000913/2070030‑2070510                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                tACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGa                                                                                                                                                                                                                                          >  2:244279/1‑249 (MQ=32) tACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATGCAGGGGGCCAGTCTAGTGGCGa                                                                                                                                                                                                                                          <  2:189722/249‑1 (MQ=17)               aaCTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGACGCATCATCCCGTGtt                                                                                                                                                                                                                          >  2:253569/1‑251 (MQ=25)                 ctggctggCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGtttt                                                                                                                                                                                                                        >  1:19794/1‑251 (MQ=40)                              acacCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTg                                                                                                                                                                                                             <  2:266607/249‑1 (MQ=255)                                acCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCCAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATgg                                                                                                                                                                                                         >  1:399001/1‑251 (MQ=255)                                     gTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCag                                                                                                                                                                                                     <  1:205459/250‑1 (MQ=255)                                                                      attCAAAGCTCTGGTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCGCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTATGGGCGAAGCCTCCTCCCGTGTTTTACTTCTCATTGATGGTCAGGAGGTAACTTCTCAGCGCTCCGGAGATAATTATg                                                                                                                                                                   >  2:226215/2‑251 (MQ=17)                                                                         cAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGgtg                                                                                                                                                                <  1:297377/251‑1 (MQ=255)                                                                                  gctgAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGt                                                                                                                                                       >  1:82415/1‑251 (MQ=255)                                                                                       aCAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTa                                                                                                                                                                     >  2:56708/1‑232 (MQ=255)                                                                                       aCAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTa                                                                                                                                                                     <  1:56708/232‑1 (MQ=255)                                                                                        cAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAg                                                                                                                                                 <  1:67115/251‑1 (MQ=255)                                                                                                 tCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTg                                                                                                                                        >  1:95574/1‑251 (MQ=255)                                                                                                           tgtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCg                                                                                                                              >  1:239572/1‑251 (MQ=255)                                                                                                           tgtgGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCg                                                                                                                              >  2:136749/1‑251 (MQ=255)                                                                                                                                        cTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGTGTTAATTCTGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAATGGTCCATATTCCg                                                                                                 >  2:395547/1‑251 (MQ=255)                                                                                                                                            gACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACt                                                                                             <  2:167503/251‑1 (MQ=255)                                                                                                                                              cATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGt                                                                                           <  2:305202/251‑1 (MQ=255)                                                                                                                                                               ccAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAAt                                                                           <  1:136749/250‑1 (MQ=255)                                                                                                                                                                    gACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCgg                                                                     >  1:290914/1‑251 (MQ=255)                                                                                                                                                                        aGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTa                                                                   <  2:399001/249‑1 (MQ=255)                                                                                                                                                                         gTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGCGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTg                                                                >  1:261808/1‑251 (MQ=255)                                                                                                                                                                            gacgtgtgctctccgatcttATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTa                                                             <  1:187278/232‑1 (MQ=255)                                                                                                                                                                             aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGt                                                               <  1:229818/248‑1 (MQ=255)                                                                                                                                                                                            gACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGGTAACTTCCTCACCaaaa                                              >  2:235228/1‑250 (MQ=255)                                                                                                                                                                                               tATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGCTGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTa                                                             >  2:187278/1‑232 (MQ=255)                                                                                                                                                                                               tATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAggg                                          >  2:133331/1‑251 (MQ=255)                                                                                                                                                                                                aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAggg                                          <  2:19794/250‑1 (MQ=255)                                                                                                                                                                                                    cTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATCGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAAt                                     >  1:196576/1‑251 (MQ=255)                                                                                                                                                                                                            aGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGAc                             <  2:251464/251‑1 (MQ=255)                                                                                                                                                                                                             gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTAtt                                                                 >  1:354030/1‑214 (MQ=255)                                                                                                                                                                                                             gTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTAtt                                                                 <  2:354030/214‑1 (MQ=255)                                                                                                                                                                                                                     ggTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATTAGTCTGCGCTGTAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCtata                     <  1:162841/250‑1 (MQ=255)                                                                                                                                                                                                                             ggAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTaa             <  1:244279/250‑1 (MQ=255)                                                                                                                                                                                                                                gaTTCAGGGGGCCAGTCTAGTGGCGAAGCATCATACAGTGTTTTAATTCTCATTGATGGTAAGGAGGTAAATTATCAGCGCGCCGGAGATACTTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTAACAGGCAATTGGCGGAATTGTTAACTTCATCACCAAAAAGGGTGAATATCCGGACTCCCTATATCAATTAAATTg         >  1:258921/1‑251 (MQ=255)                                                                                                                                                                                                                                    cAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATGGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATa      <  2:347718/250‑1 (MQ=255)                                                                                                                                                                                                                                       ggggCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGGGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGGTTTTTTAACTTCATCACCAAAAAGGGTGAATTTCCGGACTCCCTATATCACTTTAATTTATACAAc  >  2:222340/1‑251 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                TACGCCAGAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAAC  >  NC_000913/2070030‑2070510

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