Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,586,576	IS5 (–) +4 bp	intergenic (‑90/+241)	ydeP ← / ← ydeQ	putative oxidoreductase YdeP/putative fimbrial adhesin protein YdeQ

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1396044 =	NA (NA)	34 (0.710)	33/486	0.5	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	1586576 =	0 (0.000)					intergenic (‑90/+244)	ydeP/ydeQ	putative oxidoreductase YdeP/putative fimbrial adhesin protein YdeQ
*	?	NC_000913	= 1397238	NA (NA)	28 (0.590)	27/484	0.7	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	= 1586579	0 (0.000)					intergenic (‑93/+241)	ydeP/ydeQ	putative oxidoreductase YdeP/putative fimbrial adhesin protein YdeQ

TGCTGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTCGGTTGTGACTCGCGCGACCATCGGCGGCGTTATAGAA  >  NC_000913/1396998‑1397481                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    tGCTGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGGTTCGGGCGGACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAAt                                                                                                                                                                                                                                            <  1:88841/250‑1 (MQ=14)                       ggCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCa                                                                                                                                                                                                                      <  1:313964/250‑1 (MQ=255)                       ggCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGAACAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAGATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAg                                                                                                                                                                                                                     <  2:230176/251‑1 (MQ=255)                             gTTACTCACGCTGGCCAACCTTATTTGGGCGGACCAAATGATC‑CGTCAGTACGAGAGATCTCACTAAAATCTGGGAATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGGTTATCCAAGTCATTTACGGGAGAAAAAATCGTCTTAAACATGAAGAAATGAAATGTCTCAGTCAGTAGAGAAGAATTTCCCCGCTTATTAGCCGCTTCCCCAACTAAGAAAGGCGCTGATTGTAAATCCaacgaaga                                                                                                                                                                                                               <  2:37027/251‑8 (MQ=2)                                  tCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCCTTTACGGGAGAAAAACTCGGCTCAAACATGAAGAAATGCAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATAAATGCGTTGATTGTAAATCCAGCTagacgggg                                                                                                                                                                                                           >  2:498666/1‑243 (MQ=255)                                  tCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCt                                                                                                                                                                                                                   <  1:144012/242‑1 (MQ=255)                                  tCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCt                                                                                                                                                                                                                   >  2:144012/1‑242 (MQ=255)                                         ggCCAACCTGTTTCGGGCGGACCAAACGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGtttt                                                                                                                                                                                                    <  1:368230/250‑1 (MQ=255)                                                  gttacGGGCGGACCAAA‑GTTAACGGCAGTGGGAGAGATCTCACTAAAAACTGGTGATAACGCCTTAAATGGCGAAGAAACGGGCTAAAGAGGCGGATTCAAGGCATTTACGGGAGAAAAAATCGGCGCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATGCGCACCTTTCCTAACTTATCAATGCGTGGATTGTAAATCCAGCTAACAGGTGAGGGTTTCAGAGCAGAc                                                                                                                                                                                          <  2:473120/247‑1 (MQ=14)                                                      cGGGCGGACCAAATGATA‑CGTCAGTGGGATAGATCTCACTAAAAACTGGGGATAACGCCGTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGTGAGAAAAAAGCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCGAACTAATCAATGCTTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAAc                                                                                                                                                                                       <  2:440445/250‑1 (MQ=255)                                                      cGGGCGGACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGGTTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATACAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCACCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGAGAGGTTTTCagataatacac                                                                                                                                                                                        >  2:149830/1‑241 (MQ=17)                                                           ggACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCTCCTTCCCTAACTAATCAATGCGTTg                                                                                                                                                                                                                                  >  2:185043/1‑202 (MQ=38)                                                           ggACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTg                                                                                                                                                                                                                                  <  1:185043/202‑1 (MQ=38)                                                            gACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGaa                                                                                                                                                                                 >  1:227533/1‑250 (MQ=255)                                                                                    gagaTCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCATCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAg                                                                                                                                                         <  1:257907/251‑1 (MQ=255)                                                                                      gaTCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGt                                                                                                                                                       >  1:246076/1‑251 (MQ=255)                                                                                                       gggATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTTAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCTTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGGTACGTTTATGTAACAAGAAATTTGCCTg                                                                                                                                      <  2:106522/251‑1 (MQ=255)                                                                                                           aaacgacctAAATGGCGAAGAAACGGGCTAAATAGGCTGATGCAAGGCTGTTTCGGGCGCAAAAAGCGGCTCAAACATGAAGAGAGGAAATTACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCGTACCTATCTAAGCAATGCGTTGATTGTAAAACCATCTAAGAGTTGAGGTATTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTGCGTTTATGTAACAAGAAATTTGCCTGcaca                                                                                                                                  <  2:53910/243‑1 (MQ=18)                                                                                                                   tAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGCCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCGTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTGACGTTTAGGTAACAAGATATTTGTCTGCACAAGGAtt                                                                                                                            <  2:312169/249‑1 (MQ=255)                                                                                                                        ggCGAAGAAACGGTCTAAATAGGCTGATTTAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTTCACAAGGATTAcatca                                                                                                                      <  2:500385/250‑1 (MQ=255)                                                                                                                           gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGAt                                                                                                                  >  2:377775/1‑251 (MQ=255)                                                                                                                           gaagaaACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGAATGAGTCAGCCGAGCAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTACATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGAt                                                                                                                  >  1:456799/1‑251 (MQ=255)                                                                                                                                aaCGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATgg                                                                                                             >  1:125749/1‑251 (MQ=255)                                                                                                                                 aCGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGc                                                                                                            <  2:246076/251‑1 (MQ=255)                                                                                                                                  cGGTCTAGATAGGCTGATTCAAGGCATTTACGGGCGACACAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGACGGATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTTCATCATGATTATGGCg                                                                                                           >  2:206797/1‑251 (MQ=37)                                                                                                                                                      aaGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGCGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAACTGTCATGGAATTGTAACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAgg                                                                                       >  1:270752/1‑251 (MQ=255)                                                                                                                                                              tACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGaaa                                                                                <  2:500289/250‑1 (MQ=255)                                                                                                                                                                                   aaaCATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCgg                                                          <  1:482237/251‑1 (MQ=255)                                                                                                                                                                                     aCATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGtt                                                        <  2:402261/251‑1 (MQ=255)                                                                                                                                                                                         gaagaaATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTGGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAATGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTc                                                     <  2:456799/250‑1 (MQ=255)                                                                                                                                                                                               atgaaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGCAGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGAAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAGACTGAAGTCAGCGATAGGGAAGAAATTTCTCTTTGATTTGTTGGCGGTTTTCattattt                                              >  1:228544/1‑250 (MQ=255)                                                                                                                                                                                               atgaaatgaCTGAGTCAGCCGAGAAGAAGTTCCCCGCTTATTCGCACCGTCCCTAACTAAACAATGCGTTGAGTGTAAAACCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTg                                              <  2:135225/251‑1 (MQ=255)                                                                                                                                                                                                  aaatgaCTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAAAGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGAGAAACTGAAGACAGCGAAAGGGAAGAAATTTATCTTTGGTTTGTTGGCGGTTTTCattaattcgg                                           >  1:380150/1‑245 (MQ=255)                                                                                                                                                                                                        cTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTCgg                                     >  1:40865/1‑251 (MQ=255)                                                                                                                                                                                                                                  gCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGCGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTAGGTTGTGACTCGCGCGACCATCgggggc           >  1:474139/1‑251 (MQ=255)                                                                                                                                                                                                                                   cTTATTCGCACCTTCACTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTCGGTTGTGACTCGCGCGACCATCGgcggcg          >  1:177549/1‑251 (MQ=255)                                                                                                                                                                                                                                    ttATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTCGGTTGTGACTCGCGCGCCCATCGgcggcgg         >  1:172007/1‑250 (MQ=255)                                                                                                                                                                                                                                           aaCCTTCCCTAACTAATCAATGCGATGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGCGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGGTTTCATTAGTGAGGCGTCGGTTGTGACTCGCGCGACCATCGGCGGCGTTATAGaa  <  2:78622/250‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                    TGCTGAAAAACGATAACCAACTGGCGATGTTATTCACGCTGGCCAACCTGTTTCGGGCGGACCAAATGATA‑CGTCAGTGGGAGAGATCTCACTAAAAACTGGGGATAACGCCTTAAATGGCGAAGAAACGGTCTAAATAGGCTGATTCAAGGCATTTACGGGAGAAAAAATCGGCTCAAACATGAAGAAATGAAATGACTGAGTCAGCCGAGAAGAATTTCCCCGCTTATTCGCACCTTCCCTAACTAATCAATGCGTTGATTGTAAATCCAGCTAAGAGGTGAGGTTTTCAGAGCAGACAACGGTGAAATGTCATGGTATTGTTACGTTTAGGTAACAAGAAATTTGTCTGCACAAGGATTACATCATGATTATGGCGAAACTGAAGTCAGCGAAAGGGAAGAAATTTCTCTTTGGTTTGTTGGCGGTTTTCATTATTGCGGCGTCGGTTGTGACTCGCGCGACCATCGGCGGCGTTATAGAA  >  NC_000913/1396998‑1397481

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