Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,093,481	IS2 (+) +5 bp	coding (606‑610/1359 nt)	dgcT →	putative diguanylate cyclase DgcT

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	381260 =	NA (NA)	47 (0.580)	43/486	0.8	93.9%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	= 1093485	3 (0.040)					coding (610/1359 nt)	dgcT	putative diguanylate cyclase DgcT
*	?	NC_000913	1093481 =	3 (0.040)	60 (0.740)	52/490	0.5	95.1%	coding (606/1359 nt)	dgcT	putative diguanylate cyclase DgcT
	?	NC_000913	= 1651548	NA (NA)					noncoding (706/706 nt)	IS2	repeat region

TTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAA  >  NC_000913/381013‑381506                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       ttATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACt                                                                                                                                                                                                                                                        >  1:532945/1‑248 (MQ=255)     cGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTgg                                                                                                                                                                                                                                                 <  2:819420/251‑1 (MQ=255)     cGCATCCGGTAGCGATAAATACCTCTGTAGTTCGATCAGACAAGGGTTGGTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGATTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTg                                                                                                                                                                                                                                                  <  2:309751/250‑1 (MQ=255)       cATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGcc                                                                                                                                                                                                                                               >  1:446844/1‑251 (MQ=255)               aGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAAt                                                                                                                                                                                                                                       <  1:327083/251‑1 (MQ=255)                gCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAAtc                                                                                                                                                                                                                                      <  1:188012/251‑1 (MQ=255)                     aaTTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTcc                                                                                                                                                                                                                                   >  1:740686/1‑249 (MQ=255)                        tCCCTCTGGAGTTCGAACAGACACGGGCTGTTCAGCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTGTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCGCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGAc                                                                                                                                                                                                                               <  2:78757/250‑1 (MQ=255)                               ggatatCGAACAGACGCGGGCTATTCATCTAGACGACGATTACTGGTACGTACTTGTGCGACTAAATATGCCTGGCGTGCAATTACCGAGTCGGCAATAAAACATTTCTGACCAATCGTTGCTTCCTCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTCTTGTGCCATGACGAATAAATGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCaa                                                                                                                                                                                                                         <  2:280470/244‑1 (MQ=255)                                      gAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACtt                                                                                                                                                                                                                 >  1:453723/1‑250 (MQ=255)                                      gAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTAAGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGCCAACCAATATCACtt                                                                                                                                                                                                                 >  2:113107/1‑250 (MQ=255)                                                        cATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAAGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCCGTCAACCAATATCACTTAAATAAGTGATAGTCTTa                                                                                                                                                                                               >  2:312675/1‑250 (MQ=255)                                                          tCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAt                                                                                                                                                                                             >  1:2377/1‑250 (MQ=255)                                                              gACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAg                                                                                                                                                                                        <  1:312675/251‑1 (MQ=255)                                                                 gatgatTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGttt                                                                                                                                                                                     >  1:677078/1‑251 (MQ=255)                                                                   tgatTACAGGTACGCACCTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTAGGCCATAAAACATTGCGAACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCCACAACAATTTTTACGTCATTGATGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAATTAAATAAGTGATAGTATTAATACTAGttttt                                                                                                                                                                                   >  2:119450/1‑251 (MQ=37)                                                                           ggTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGt                                                                                                                                                                           <  1:198220/251‑1 (MQ=255)                                                                               cgtacTTGTGCGCCTAAATATGCCTGTTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTTCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCAtt                                                                                                                                                                       <  2:37586/251‑1 (MQ=255)                                                                                                 tatGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTTCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGCtgattgat                                                                                                                                                     >  1:361468/1‑251 (MQ=255)                                                                                                      cTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCtt                                                                                                                                                >  1:849193/1‑251 (MQ=255)                                                                                                             gCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGAATTTTTAATCTCGGGGGCAAAACCAATTTTTCCGCCATTGCTGATTATTTTGACCTGACGGATAAACGCATAATTACTAGACTAGACCCCTGAATCTCCAGAAACCCAATATCACTTAAATAAATGATAGTCTTAATACTCGTTTTTAGACAAGTCATTGGCGAACAGATGATTGATGACTTAgagc                                                                                                                                           >  1:578191/1‑246 (MQ=255)                                                                                                                aCTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTTCTTCCGCTTCAATAACCGCATTTTTAATCTCGGGTGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGAGCCGgag                                                                                                                                      <  2:426334/251‑1 (MQ=255)                                                                                                                  caCCGATTCTGCAATAAAACATTGCAGACCAAGCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCAATGCTGATTATAGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGaa                                                                                                                                    <  2:205728/250‑1 (MQ=255)                                                                                                                        gTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCTCTTCAATAACCGCATTTTTAATCTCGATGGCAAAACCAATTTTTACGCCATTGCTGATTATTATGCCAGTACGAATAAACGCATAATTACTAGACTGGACCCCTGATTCTACAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAATCTAGTCATTGGAGACAAGATGATTGATGTCTTAGAGCCGGAGaacctcgg                                                                                                                              >  2:76467/1‑245 (MQ=25)                                                                                                                             cAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGc                                                                                                                         <  2:337897/251‑1 (MQ=255)                                                                                                                             cAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGc                                                                                                                         >  1:337897/1‑251 (MQ=255)                                                                                                                                           ggACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGaa                                                                                                           <  1:14053/251‑1 (MQ=255)                                                                                                                                                     ttGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAACCCAATTTTTACGACATTGCTGATTATTGTGACAGGCCGACTAAACGCATAATTACTAGACTGGCCACCTGACTCTCAAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACCGATGCTTGATGTCTTAGGGCCGGAGAAACGCAGACGGAGTGCCACACAGGAAAAGATCGCa                                                                                                  >  2:820441/1‑250 (MQ=18)                                                                                                                                                        cTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCACTTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTg                                                                                              >  1:726067/1‑251 (MQ=255)                                                                                                                                                                     aaCCGCTTTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTACAGACAAACAATATCACTTAACTAAGTGATAGTCTTAATACTAGATTTTAGACTAGTCATTGGAGTACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAAATCGCCATTGTTCCGCAGAGCtt                                                                                 >  2:192220/1‑251 (MQ=18)                                                                                                                                                                            tttttAATCTCGGTGGAAAAACCAATTTTTAAGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATACTTACTAGCCTGGCCCCCTGAATATCCAGACAACCAATATCACTTAAATAAGGGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGACCAACTGATTGATGTCTTAGGGACGGAAAAACGCCGACGTCGTACCACACAGGCAAAGATCGACATTGTACAGCAGAGCTTTGAACCg                                                                          >  2:535641/1‑251 (MQ=18)                                                                                                                                                                                tAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGATTTGAACCgggg                                                                       >  2:106172/1‑250 (MQ=255)                                                                                                                                                                                    ctcGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTAAGCAGAGCTTTGACCCGGGGATGAc                                                                  >  2:653783/1‑251 (MQ=255)                                                                                                                                                                                      cGGTTGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCATACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGTAGAGCTTTGAACCGGGGATGa                                                                   <  1:529206/248‑1 (MQ=255)                                                                                                                                                                                        gTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACAAATATCACTTAACTAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCCGACGGCGTACCACACAGGAAAAGATAGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGtc                                                              >  2:824780/1‑251 (MQ=255)                                                                                                                                                                                           gCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACca                                                                                                                   >  2:267967/1‑195 (MQ=255)                                                                                                                                                                                           gCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACca                                                                                                                   <  1:267967/195‑1 (MQ=255)                                                                                                                                                                                              aaaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTc                                                        >  1:636357/1‑251 (MQ=255)                                                                                                                                                                                               aaCCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCg                                                       <  1:409158/251‑1 (MQ=255)                                                                                                                                                                                                   aaTTTTTACGACATTGCTGATTATTGTGCCAGGACGAATCAACGCATAATTACTAGACTGGCCCCCTGAATCTCCTGACAACCAATATCAATTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAACAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGCCGGACTCCCTCGTTGc                                                   >  2:554099/1‑251 (MQ=17)                                                                                                                                                                                                     tttttACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCaga                                                                                     <  1:572988/215‑1 (MQ=255)                                                                                                                                                                                                          aCGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCaa                                            <  1:5267/251‑1 (MQ=255)                                                                                                                                                                                                              cATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGTCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATg                                        <  2:2377/251‑1 (MQ=255)                                                                                                                                                                                                               aTTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGaa                                                                                                                                    <  1:821852/158‑1 (MQ=255)                                                                                                                                                                                                               aTTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGaa                                                                                                                                    >  2:821852/1‑158 (MQ=255)                                                                                                                                                                                                                      attattGTGCCAGGGCGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTGCCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGCTGCCGGTCTCCCTCGTTGCCCGGCAACATGGTGTagc                                 >  2:132080/1‑250 (MQ=255)                                                                                                                                                                                                                              gCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAAAAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCagccagcc                         >  2:713572/1‑250 (MQ=255)                                                                                                                                                                                                                               ccAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCagccagcc                         >  1:517738/1‑249 (MQ=40)                                                                                                                                                                                                                                cAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGCCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCagccagcc                         >  2:660917/1‑248 (MQ=255)                                                                                                                                                                                                                                cAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACAAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAAc                                                                            <  2:405768/197‑1 (MQ=255)                                                                                                                                                                                                                                cAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACAAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAAc                                                                            >  1:405768/1‑197 (MQ=255)                                                                                                                                                                                                                                   gACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTAt                   <  1:529456/251‑1 (MQ=38)                                                                                                                                                                                                                                     cGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTAttt                 >  1:261461/1‑251 (MQ=37)                                                                                                                                                                                                                                               cATAATTACTAGACTGGCCTCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTTGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGGTGTGAACCGGGGATGACGGTCTCCTTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTa       <  1:81625/251‑1 (MQ=255)                                                                                                                                                                                                                                               cATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTa       <  1:42635/251‑1 (MQ=32)                                                                                                                                                                                                                                               cATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTAATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGt        >  1:813449/1‑250 (MQ=17)                                                                                                                                                                                                                                                  aaTTACTAGTCCTGCCCCCTGAATCTTCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCATACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCATCCAGCCAGTTATTTCTCTGTCGTAAg     <  2:585459/250‑1 (MQ=11)                                                                                                                                                                                                                                                  aaTTACTAGACTGGCCCCCTGAATCTCCAGACAACAAATATCACTTAAATACGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACCGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGAGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCAGGACACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGc    >  2:804883/1‑251 (MQ=11)                                                                                                                                                                                                                                                   aTTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGTGCACAGATGATTGATGTCTTAGTTCCGGAGAAACGCAGACGGCGTACCACACAGGAAACGAACGCAATGGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGAAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGtcaga    >  2:852062/1‑246 (MQ=11)                                                                                                                                                                                                                                                   aTTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGc    >  1:219820/1‑250 (MQ=31)                                                                                                                                                                                                                                                     tACTTGACTGGACTCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCGTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCaa  <  2:446844/250‑1 (MQ=2)                                                                                                                                                                                                                                                     tACTAGACGGGCCCCCTGAATCTCCAGACACCCAATATCACTTAAATAAGTGATAGTCTTAATACTAGGTTTTACACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGATTTTAACCGGGGATACCGGTCTCCATCTTTGCCCGGCACCATGGTGAAGCAGCCAGCCAGTTATTTCTACGGCGTacgaac  >  2:683521/1‑245 (MQ=11)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                       TTATCGCATCCGGTAGCGATAATTCCCTCTGGAGTTCGAACAGACACGGGCTGTTCATCCAGACGATGATTACTGGTACGTACTTGTGCGCCTAAATATGCCTGGTTTGCAACTACCGAGTCGGCAATAAAACATTGCGGACCAATCGTTGCTTCCGCTTCAATAACCGCATTTTTAATCTCGGTGGCAAAACCAATTTTTACGCCATTGCTGATTATTGTGCCAGGACGAATAAACGCATAATTACTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGATGTCTTAGGGCCGGAGAAACGCAGACGGCGTACCACACAGGAAAAGATCGCAATTGTTCAGCAGAGCTTTGAACCGGGGATGACGGTCTCCCTCGTTGCCCGGCAACATGGTGTAGCAGCCAGCCAGTTATTTCTCTGGCGTAAGCAA  >  NC_000913/381013‑381506

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