Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,649,932	IS2 (+) +5 bp	coding (139‑143/192 nt)	ydfD →	lysis protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1467910 =	NA (NA)	29 (0.500)	26/492	0.8	100%	noncoding (1331/1331 nt)	IS2	repeat region
	?	NC_000913	1649932 =	0 (0.000)					coding (139/192 nt)	ydfD	lysis protein
*	?	NC_000913	= 1649936	0 (0.000)	32 (0.560)	30/492	0.6	100%	coding (143/192 nt)	ydfD	lysis protein
	?	NC_000913	= 2070271	NA (NA)					noncoding (1/1331 nt)	IS2	repeat region

GAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAACTTTTTAGAGT  >  NC_000913/2070037‑2070520                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          gagaAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCAtcctc                                                                                                                                                                                                                                           <  1:136593/251‑1 (MQ=35)   gaAATTACTGGCTGTCAGCTACACCATGTTGCCGGGCAACGCGGGAGACCGTCATCCCCGGTTCAAAGCTCTTCTGAACAAGTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACAACAATCATCTGTTCTCCAATTACTAGTCTAAAATCTAGGATTAAGACTATCGCTTATTTAAGTGATATTGTTTGTCTGGAGATTCAGGGGGCCATTCTAGTGGCGAAGCATCCTccc                                                                                                                                                                                                                                         <  2:275007/251‑1 (MQ=11)                 gctgctACACCATGTTGACGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATtctc                                                                                                                                                                                                                           >  2:250066/1‑251 (MQ=255)                   tgctACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCAt                                                                                                                                                                                                                         >  1:18337/1‑251 (MQ=255)                        cacCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAGGCCATCAATCATCTGTTCTCCAATGACTCGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGTGATTCCGGGGGCCAGTCTAGTGGCGAACCATCCTCCCGTGTTTTAATTCTCCTTGATg                                                                                                                                                                                                                    >  1:316300/1‑251 (MQ=14)                         acCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATATGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATgg                                                                                                                                                                                                                   >  1:579479/1‑251 (MQ=255)                           cATGTTGCCGGGCAACGAGGGAGATCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGGTTTTATTCTCATTGATGGTc                                                                                                                                                                                                                 <  2:501828/251‑1 (MQ=255)                                 gCCGGGCAACGAGGGAGTCCGTCATCCCCGTCTCAAAGCTCTGCTGAACAATTCCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGGATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCaggag                                                                                                                                                                                                            <  2:460651/250‑1 (MQ=255)                                         aCGAGGGAGACCTTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCTATCTGTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGTCCCTATGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACGAGTATTAAGTCTATCACTTATTTAAGTTATAGTTTTTGTCTGTTGGTTCAGGGGTCCAGTCTAGTGGCGAATCAGCCTCGCTTGTGTTAATTCTGGTGGATGGTCagcaataaactt                                                                                                                                                                                                     <  2:535958/249‑11 (MQ=2)                                             gggAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAgc                                                                                                                                                                                               >  1:358501/1‑251 (MQ=255)                                                agaCCGGCATCCGCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTGCCTGTTTGGTACGCCGTCTGCGTTTCGCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTAGTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGGGGTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAgcgcg                                                                                                                                                                                            <  2:358501/251‑1 (MQ=17)                                                        aTCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGAt                                                                                                                                                                                    <  1:338340/251‑1 (MQ=255)                                                          ccccGGTTCAAAGCTCTGCTGAACAATTGCGATTTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCTCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGTGGCCAGTCTAGTGGCTAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATa                                                                                                                                                                                   <  2:367764/250‑1 (MQ=255)                                                           cccGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAAt                                                                                                                                                                                 >  1:128196/1‑251 (MQ=255)                                                                  cAAAGCTCTGCTGAACAATTGCTATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTTGTTGTCTGGAGATTCAGGGGGCCAGTCTACTGTCGAAGCAACCTCCCGTGTGTTAAtt                                                                                                                                                                                                                              <  2:552697/199‑1 (MQ=11)                                                                  cAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAAtt                                                                                                                                                                                                                              >  1:552697/1‑199 (MQ=255)                                                                    aaGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGtgtg                                                                                                                                                                        >  1:361398/1‑251 (MQ=255)                                                                     aGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTgg                                                                                                                                                                       <  1:322766/251‑1 (MQ=255)                                                                         ctgctgAACAATCGCGATATTTTCCTGAGTGGCACGCCGGCTGCGTTTCTCCGGCCCTAACACAACAATCATCTGTTCTCCAATGACTAGTATCAAAACTAGTATTAACACTATCACTTATTTACGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGCACCATCCTACCGTGTTTTAATTCTCATTGATGGTCAGGCGGTAACTTATCCGCGCGCCGGCGACAATTATGGGGTGGGACt                                                                                                                                                                   >  2:65894/1‑251 (MQ=17)                                                                                                   gtgtgGTACGAAGTCTGCGTTTCTCCGGACCTAAGACATCAATCCTCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGAATATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTCGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGCCTGTTGATCGATGAGTCTGCGCTGGAGc                                                                                                                                         >  2:517326/1‑251 (MQ=25)                                                                                                              cGTCAGCGTTTCTCCGAACCTAAGACATCACTCATCTGTTCTCCAATGCCTAGTCTAAAAACTAGTATTAAGACTATCACTTCTTTAAGTGATATTGGATGTCTGGAGATACAGGGGGACAGTATAGTGGCGAAGCATCATCCCGTGTTTTAATTCTCATTGATGGTGAGGAGGTAACATATCAGAGCACCGGAGATAATTATGGTGTGAGACTGTTGATAGATGAGTCTGAGCGGGTGAGTGTTGAgtt                                                                                                                               >  2:229207/1‑248 (MQ=21)                                                                                                                                     gACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACt                                                                                                       <  1:346012/251‑1 (MQ=255)                                                                                                                                                   gTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAgg                                                                                         >  1:330646/1‑251 (MQ=255)                                                                                                                                                     tctcCAATGACTAGTCTAAAAACTAGTATTAAGACTCTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGCGGTAACTTATCAGCGCGCCGGAGATAATTATGTTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGCGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGc                                                                                        >  2:268470/1‑250 (MQ=255)                                                                                                                                                           aTGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGc                                                                                 >  1:186384/1‑251 (MQ=255)                                                                                                                                                                      aaaaaCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGAGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCAGAGGCAATTGGCGGTATTGTTa                                                                       >  2:493990/1‑250 (MQ=255)                                                                                                                                                                          aCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCg                                                                                >  2:164431/1‑237 (MQ=255)                                                                                                                                                                          aCTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATAGTCCGTACTGTACGGTTCACAGGCAATTGGCg                                                                                <  1:164431/237‑1 (MQ=255)                                                                                                                                                                             aGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGGTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGGTCACAGGCAATTTGCGGTATTGTTAACGtcatc                                                               <  2:361398/251‑1 (MQ=255)                                                                                                                                                                             aGTATTAAGACTAGCACTTATTTAAGTGATATTTGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGTAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGAGAGATGAGTCTGCGCTGGAGCGTGTTGAGTTAGTGAAAGTTCCATATGCCGTGCTGTACTGGTCACAGGCAAGTGGCGGTATTGTTAACTtca                                                                 <  2:579479/249‑1 (MQ=38)                                                                                                                                                                                     gACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACGTCCTCACCaaaa                                                        >  2:163597/1‑250 (MQ=255)                                                                                                                                                                                      aCTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTgg                                                                                                                                            <  2:224693/165‑1 (MQ=255)                                                                                                                                                                                      aCTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTgg                                                                                                                                            >  1:224693/1‑165 (MQ=255)                                                                                                                                                                                         aTCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAgg                                                     <  2:160417/249‑1 (MQ=255)                                                                                                                                                                                            aCTTAATTAAGTGCTATTCGTTGTCTGGAGATTCAGGGGGCCCGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATACTTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAggggga                                                 >  2:555155/1‑247 (MQ=255)                                                                                                                                                                                                   tAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCg                                          <  1:143218/250‑1 (MQ=255)                                                                                                                                                                                                         atatTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGTCTccc                                   >  1:13905/1‑251 (MQ=255)                                                                                                                                                                                                               gTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGCTGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATc                             >  1:269615/1‑251 (MQ=255)                                                                                                                                                                                                                    cTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGAATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGAATTGTTAACTTCATCACCCAAAAGGGGGAATCTCCGGACTCCCTATATCACtt                         >  2:451289/1‑250 (MQ=255)                                                                                                                                                                                                                          aTTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTg                   >  1:571885/1‑250 (MQ=255)                                                                                                                                                                                                                              aGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTCTTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATAc               >  1:43775/1‑250 (MQ=255)                                                                                                                                                                                                                               gggggCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACa              >  1:85358/1‑250 (MQ=255)                                                                                                                                                                                                                                           aGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAACTTTTTAGagt  <  1:250066/250‑1 (MQ=255)                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                          GAGAAATAACTGGCTGGCTGCTACACCATGTTGCCGGGCAACGAGGGAGACCGTCATCCCCGGTTCAAAGCTCTGCTGAACAATTGCGATCTTTTCCTGTGTGGTACGCCGTCTGCGTTTCTCCGGCCCTAAGACATCAATCATCTGTTCTCCAATGACTAGTCTAAAAACTAGTATTAAGACTATCACTTATTTAAGTGATATTGGTTGTCTGGAGATTCAGGGGGCCAGTCTAGTGGCGAAGCATCCTCCCGTGTTTTAATTCTCATTGATGGTCAGGAGGTAACTTATCAGCGCGCCGGAGATAATTATGGTGTGGGACTGTTGATAGATGAGTCTGCGCTGGAGCGTGTTGAGGTAGTGAAAGGTCCATATTCCGTACTGTACGGTTCACAGGCAATTGGCGGTATTGTTAACTTCATCACCAAAAAGGGTGAATCTCCGGACTCCCTATATCACTTAAATTGATACAACTTTTTAGAGT  >  NC_000913/2070037‑2070520

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