Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,542,042	IS1 (+) +8 bp	coding (6‑13/597 nt)	fimE →	regulator for fimA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	41 (0.810)	37/482	0.3	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 4542049	0 (0.000)					coding (13/597 nt)	fimE	regulator for fimA
*	?	NC_000913	= 3584195	NA (NA)	58 (1.150)	53/484	0.1	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	4542042 =	0 (0.000)					coding (6/597 nt)	fimE	regulator for fimA

CTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTG  >  NC_000913/3583952‑3584440                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      cTTTGACGTGGTGATATAGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACAtt                                                                                                                                                                                                                                                 >  1:158185/1‑250 (MQ=17)       cGTGGTGATATAGATGACGGATGGCCGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTa                                                                                                                                                                                                                                          <  1:108155/251‑1 (MQ=25)          ggtgATATAGATGACGGATGGCTGTCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAAGCTGAGGCAGCACCTGGCACGGCGGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑TAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACt                                                                                                                                                                                                                                       <  2:344511/251‑1 (MQ=255)              atatatattagggagggctggcCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAAGCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCCCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTGCCGACATTAAAAGTAACtggtt                                                                                                                                                                                                                                   <  2:305763/236‑5 (MQ=11)               tatagATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACAGACATTAAAAGAAACtttttt                                                                                                                                                                                                                                  >  2:327411/5‑251 (MQ=255)                        cGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAACGTAACTTTTTTGAAAtt                                                                                                                                                                                                                            >  1:5296/1‑248 (MQ=255)                          gaggtatggcGGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCCCGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCGGAGGCAGCCCCGGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga                                                                                                                                                                                                                       <  2:211545/245‑1 (MQ=14)                          gATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAgaga                                                                                                                                                                                                                       <  1:135941/251‑1 (MQ=255)                           aTGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTag                                                                                                                                                                                                                          >  1:72055/1‑247 (MQ=255)                                    cgcagtaggAATCCCGCCTGAAGTGAAAGCTGTACGTAATCATCAAGCGCTATACGCCGCGAATTGAGCGGCATAACCTGAATCCAAGGCAGCCTCTGGCACGGACGGGACGGAAGTCGCTGTTGTTCTCAAAATCG‑TGTGAGCGGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATAAATAAGTGGGAGTCATTACCGACATTAAAAGTAACTTTTTTGACATTAGAGATAATATTGGc                                                                                                                                                                                                             <  2:8502/243‑1 (MQ=14)                                     gCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                                                                                                                                                                                                                   >  2:555458/1‑228 (MQ=40)                                     gCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACttttt                                                                                                                                                                                                                                   <  1:555458/228‑1 (MQ=40)                                     gCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc                                                                                                                                                                                                             >  1:29698/1‑250 (MQ=255)                                        gTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAg                                                                                                                                                                                                         <  1:441742/251‑1 (MQ=255)                                           agatacccGCCTGAAGTGAAAGCTGCACTTAATCAGCAAGCGAGATACGCAGCTAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGGCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAAAAAGTGGGATTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGaaa                                                                                                                                                                                                      <  2:403544/246‑1 (MQ=17)                                                    ccTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTc                                                                                                                                                                                              >  2:89354/1‑250 (MQ=255)                                                       gAAGGGAAAGCTGCACGGAATCAGCAAGCGATATACGCAGCGAATTGGGCTGCATAACCTGAATTGGAGGCAGCTCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCGCAAAATCGTTG‑GTGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCAGTACCGACATTAAAATTATCTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGtt                                                                                                                                                                                          <  2:216049/251‑1 (MQ=17)                                                                    cACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGAACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCATTTAGTGtttt                                                                                                                                                                             >  2:543550/1‑251 (MQ=255)                                                                     aCGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGCCAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTCGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGAtttt                                                                                                                                                                            >  2:61921/1‑250 (MQ=255)                                                                     aCGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTtc                                                                                                                                                                            <  1:572263/251‑1 (MQ=255)                                                                      cGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTtct                                                                                                                                                                           <  1:476766/251‑1 (MQ=255)                                                                              gCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCAttt                                                                                                                                                                   >  2:429117/1‑250 (MQ=255)                                                                                   ggATATTCGCAGCGAATTGATCGGCAGAACCTGATTCTGAGGCTGCACCTGGCACGGCTGGGACGGAAGCTGCTGTCGTTCTCAAAATCTGTA‑GCGCTGCATGATAAAGTCAACTGGCATTATTTGAACATAAATCACTTTCAATAAGTGGGAGGCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGAGAATATTGGCACAGAAAATATATTCAGTGCGTTTAGTGATTTCGCTTCATTTTTaaa                                                                                                                                                              <  2:335787/250‑1 (MQ=255)                                                                                              gCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAACGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGTAATCATgg                                                                                                                                                   >  1:165796/1‑251 (MQ=255)                                                                                                 aaTTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACATAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACa                                                                                                                                                <  1:279261/251‑1 (MQ=255)                                                                                                  agtGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTctct                                                                                                                                                                         <  2:388543/223‑1 (MQ=255)                                                                                                  aTTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTctct                                                                                                                                                                         >  1:388543/1‑225 (MQ=255)                                                                                                            aTAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCata                                                                                                                                      >  1:293953/1‑250 (MQ=255)                                                                                                               aCCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCAtaataa                                                                                                                                  <  1:61921/251‑1 (MQ=255)                                                                                                               aCCTGAATATGAGGCAGCACCTGGCACGGCTGGGCCGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCCTTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTAACGACATTAAAAGTAACTTTTTTGAACTTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCTTTCTTAAACGAAATCATGAACAGTTGTTCatatta                                                                                                                                   >  2:45861/1‑247 (MQ=25)                                                                                                                   gAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACCGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCa                                                                                                                               >  2:506326/1‑250 (MQ=255)                                                                                                                   gAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATc                                                                                                                                <  2:171070/249‑1 (MQ=255)                                                                                                                                                gACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGc                                                                                                 >  1:150307/1‑251 (MQ=255)                                                                                                                                                       gTCGCTGTCGTTCTCAAAAGCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCttt                                                                                          <  1:308321/251‑1 (MQ=255)                                                                                                                                                          gCTGTCGTTTCCAAAATCGGTG‑GAGCTGCATGACAAATTCATCGGGCATTATCTGAACATAAAACACTAACAATAAGGTGGAGTCATTACCGACATTAAAAGTAACTGTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATTATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTa                                                                                       <  2:30796/251‑1 (MQ=38)                                                                                                                                                           cTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCtt                                                                                                                                                                        >  1:63040/1‑169 (MQ=255)                                                                                                                                                                           cGGTG‑GAGCTGCATGACAAAGGCATCGGGCAGTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTATAGATTACCAGCAGaaa                                                                      <  2:450576/251‑1 (MQ=255)                                                                                                                                                                              tG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAAtt                                                                   >  1:326073/1‑251 (MQ=255)                                                                                                                                                                                                                       aaaaCACTAACAATAAGTTGGAATCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATAATGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGTAAAAACtt                           <  2:507926/251‑1 (MQ=255)                                                                                                                                                                                                                             cgaTCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAATCAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACg                     <  2:158185/249‑1 (MQ=255)                                                                                                                                                                                                                                          ggAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCt        >  1:235946/1‑251 (MQ=255)                                                                                                                                                                                                                                           gAGTCATTCCCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTTGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCtt       <  2:156670/251‑1 (MQ=255)                                                                                                                                                                                                                                                aTTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTg  <  1:429117/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                      CTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTG‑GAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTG  >  NC_000913/3583952‑3584440

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