Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,541,425	IS1 (–) +8 bp	coding (469‑476/603 nt)	fimB →	Type 1 fimbriae regulatory protein FimB

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	50 (1.050)	44/486	0.2	96.0%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	4541425 =	2 (0.040)					coding (469/603 nt)	fimB	Type 1 fimbriae regulatory protein FimB
*	?	NC_000913	= 3584195	NA (NA)	57 (1.200)	49/486	0.1	96.5%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	= 4541432	2 (0.040)					coding (476/603 nt)	fimB	Type 1 fimbriae regulatory protein FimB

CCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTG  >  NC_000913/3583950‑3584440                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      ccggttggagtggtgAAAAGGATGACGGAGGGCGGGCCGCGGTATGAAACCCTCCTGAAGGGAAAGCGGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACTGCTGGGACGGAAGTCGCTGTCGTTGTCTAAAGCGGAGGAGCTGCATGACAATTTCAGCGGGCATTATCTGAACATACAACACTATCAATAATTTTGAGTCATTACCGACAt                                                                                                                                                                                                                                                  <  2:190898/242‑1 (MQ=11)     ttGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTa                                                                                                                                                                                                                                                <  1:447405/249‑1 (MQ=32)       gACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCGGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAg                                                                                                                                                                                                                                            <  2:229168/251‑1 (MQ=21)         cGTGGAGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGACTCTGAAGCCGCACCTGGCACGGCTGGGACGGAAATCGCTGTCGTTCTCAAAATCGGTGGCGCTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATACATAAGTTGGAGTCATTACCGACATTAAAAga                                                                                                                                                                                                                                           >  2:352163/1‑249 (MQ=11)                  atGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACtttttt                                                                                                                                                                                                                                  <  1:271181/250‑1 (MQ=39)                    ggaggACGGAGGGCTGGCCGCTGAATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaa                                                                                                                                                                                                                               <  2:1759/247‑1 (MQ=255)                     gATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGaaa                                                                                                                                                                                                                              <  1:144227/251‑1 (MQ=40)                       ggaaggATGGCGGGCGGCGGTATGAACCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCGGGGACGGAAGTCGCTGTCGTTCTAAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACAGAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                                                                                                                                                                                                            <  2:111584/247‑1 (MQ=14)                                   ggcCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAatat                                                                                                                                                                                                                 <  1:496451/250‑1 (MQ=255)                                      cGCTGTAAGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGc                                                                                                                                                                                                             <  2:122716/251‑1 (MQ=255)                                          gTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGACATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGcac                                                                                                                                                                                                           <  1:352163/249‑1 (MQ=255)                                          gTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGAAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGCCATCGGGCATTATCTGAACATAAAACAATATCAATAAGTTGGAGCCATTACCGACATTAAAAGTAAATTTTTTGAAATTAGAGATAATATTGacac                                                                                                                                                                                                           >  2:461692/1‑249 (MQ=255)                                              gatgcccGTCTGAAGGGATAGCTGCACGTAATCATCAATCGATATACGCAGCGAATTGAGCGGCATAACCTTAATCTGAGTCAGCACCTGGCACGGCGGTTAAGGATGGCGCTGGCGTCCTCAAAATTGGTGGAGCTGTATTACAAAGTCATCGGGCATTAGCTGAAAATAAAACCCTATCAATAATTTGGAGTCATGACCGACATTAAAAGTAACTTTTTTGAAAGTAGAGCTAATATAGGCCCATaaaa                                                                                                                                                                                                     <  2:488978/247‑1 (MQ=9)                                                      ccTGAAGGGAAAGCGGCCCGTAATCAGCAAGCGATATACGCAGCGAAGTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTAACAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCa                                                                                                                                                                                             <  2:362923/251‑1 (MQ=255)                                                           agtgAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCg                                                                                                                                                                                        <  1:188503/248‑1 (MQ=255)                                                                           aaTCAGCAAGCGATATACGCAGCGTATTGGGCGGCATAACCTGAATCTGAGGCAGCTCCTGGCACGGCTGGGACGGTAGTCTCTTTCGTTCTCAAAATCGGTTGAGATGCATTATAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGTAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCtt                                                                                                                                                                        <  2:382312/251‑1 (MQ=21)                                                                                 cgaGCGATATACGCAGCGAATTGAGCGTCATAACCTGAATCTGAGGCAGCACCGTGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGTTGGAGCTTAATGACAAAGTCATCGGGCATTATCTGAACATAAATCACTAGCAATAAGTTGGATTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATAGTCAGTTCGTTTAGTGATTGCTCGTCATTCt                                                                                                                                                                  <  2:290675/249‑1 (MQ=18)                                                                                 cAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGTCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCt                                                                                                                                                                  >  2:141378/1‑251 (MQ=255)                                                                                       atatACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACg                                                                                                                                                            >  1:358477/1‑251 (MQ=255)                                                                                            cGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAAGCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATc                                                                                                                                                       <  1:38841/251‑1 (MQ=255)                                                                                             gCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCTGTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTGCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCa                                                                                                                                                      <  2:77938/251‑1 (MQ=255)                                                                                                  gAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGAc                                                                                                                                                 >  1:495659/1‑251 (MQ=255)                                                                                                   aaTTGAGCGGCATAACCTGAATCTGAGTCAGCACCTGGCACGGCTGGGACGGAAGTCTCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACa                                                                                                                                                <  1:330930/251‑1 (MQ=255)                                                                                                              aTAACCTGAATCTGAGGCATCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGATCTCAAAAACGGTGGAGCTGCATTACAAAGTCATCGTGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCATAGAAAATATATTCAGTTCGTTTAGTGATTGCTCTGCATTCTTAAACGAAATCATGTACCGTTGTTCataa                                                                                                                                     <  2:42261/251‑1 (MQ=21)                                                                                                              aTAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCataa                                                                                                                                     <  1:428607/251‑1 (MQ=255)                                                                                                                     gAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCa                                                                                                                               >  1:371135/1‑250 (MQ=255)                                                                                                                            aGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGTTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACGTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGTAATCATGGACAGTTGTTCATAATAATCAAAAGGaaa                                                                                                                       <  1:260724/251‑1 (MQ=255)                                                                                                                             ggCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGCCAGTTGTTCATAATAATCAAAAGGAAAc                                                                                                                      >  2:292431/1‑251 (MQ=255)                                                                                                                                 gcaCCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACaaaa                                                                                                                  <  1:331650/251‑1 (MQ=255)                                                                                                                                       ggCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCCTGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCAtt                                                                                                            >  2:325418/1‑251 (MQ=255)                                                                                                                                       ggCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAAtt                                                                                                                                                                                                                            <  1:498648/139‑1 (MQ=255)                                                                                                                                       ggCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGAAACTTTTTTGAAAtt                                                                                                                                                                                                                            >  2:498648/1‑139 (MQ=255)                                                                                                                                        gCACGGCTGGGACGGAAGGCGCTGTCGGTCTCAAAATCGGTGGAGCTGCAGGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAATTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTg                                                                                                           <  1:15186/251‑1 (MQ=255)                                                                                                                                         cACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTAATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGt                                                                                                          <  1:365859/251‑1 (MQ=255)                                                                                                                                            ggCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTaa                                                                                                        >  1:48573/1‑250 (MQ=255)                                                                                                                                                         gTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAATGCATCGGGCATTATCTGAAAATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTGTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCttt                                                                                          <  2:455537/251‑1 (MQ=255)                                                                                                                                                                                                    accGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAg                                               <  2:424401/249‑1 (MQ=255)                                                                                                                                                                                                        ggCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCa                                           <  1:65605/251‑1 (MQ=255)                                                                                                                                                                                                                                             aGTCATTACCGACATTAAAAGTAACGTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTa      <  1:488903/251‑1 (MQ=255)                                                                                                                                                                                                                                                     ccGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTg  >  2:296871/1‑247 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                      CCCTTTGACGTGGTGATATGGATGACGGATGGCTGGCCGCTGTATGAATCCCGCCTGAAGGGAAAGCTGCACGTAATCAGCAAGCGATATACGCAGCGAATTGAGCGGCATAACCTGAATCTGAGGCAGCACCTGGCACGGCTGGGACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGACATTAAAAGTAACTTTTTTGAAATTAGAGATAATATTGGCACAGAAAATATATTCAGTTCGTTTAGTGATTTCTCTTCATTCTTAAACGAAATCATGGACAGTTGTTCATAATAATCAAAAGGAAACAAAAATCATTGTAAACCATGCCATCTTTTTAGAGATTACCAGCAGAAAATTCCTTGTTCGAAATTTAGCAGTGCAAGCTATGGAAAAACTTACAACGGAACTACACTCCTTAAGTG  >  NC_000913/3583950‑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.