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)	43 (0.740)	42/490	0.4	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)	77 (1.310)	72/492	0.1	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	4542042 =	0 (0.000)					coding (6/597 nt)	fimE	regulator for fimA

GAAGTTATGTAATGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTT  >  NC_000913/3583178‑3583674                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       gaAGTTATGTAATGAAGTTAGTATTCTTTGTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAg                                                                                                                                                                                                                                                        <  2:327243/251‑1 (MQ=255)            aTGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGcc                                                                                                                                                                                                                                             <  2:71651/251‑1 (MQ=255)                       aTTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATAACATTAAAGGTGATGCTGCCAACTTACTGAt                                                                                                                                                                                                                                  >  1:498006/1‑251 (MQ=255)                            ttttAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGTGGATTATAGCATTTCAATATTTACTTATGCCATTAAATGTGATGCTGCCAACTTACTGATTTa                                                                                                                                                                                                                               <  2:30124/249‑1 (MQ=255)                                   aaTCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACACTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCAAACTTACTGATTTAGTGTatgat                                                                                                                                                                                                                      >  1:279189/1‑251 (MQ=255)                                      cAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATgg                                                                                                                                                                                                                    <  2:72037/250‑1 (MQ=255)                                       aaCCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGgtg                                                                                                                                                                                                                  >  1:372984/1‑251 (MQ=255)                                            gATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATGAAAGGTGAGGCTGCCAACTTACTGATTTAGTGTATGATGTTGttttt                                                                                                                                                                                                             <  2:54315/251‑1 (MQ=255)                                                    tACTTTTTTAAGAGCAATGGATGGATTCGACGTTAATGGATTACGATTATTTAGCCTCTCGATTCCTGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCt                                                                                                                                                                                                     >  2:270573/1‑251 (MQ=255)                                                                agCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCt                                                                                                                                                                                         >  1:270472/1‑251 (MQ=255)                                                                         tggatTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATc                                                                                                                                                                                >  2:485126/1‑251 (MQ=255)                                                                               cGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGt                                                                                                                                                                          <  2:209534/251‑1 (MQ=255)                                                                                             tACGATTATTTAGCCTCTCGGTTCCAGAACCTTCAGTTAAAAACCTTTTTTCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACACGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTTAGTTGCTCCAGTGCCTTCTGTTTCTATCATCTGTCCCTCCTGTTccgc                                                                                                                                                            >  2:212595/1‑248 (MQ=38)                                                                                               cGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCt                                                                                                                                                           >  1:514840/1‑250 (MQ=255)                                                                                                gATTATTTAGCCTCTCGATTCCAGAACCGTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGAAGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTGCTTATGACATTAAAGGTGATGCTGCCAACTTACTTATTTAGTGTATGATGTTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCTTTCTTGTTCAGCTAc                                                                                                                                                         <  2:365122/251‑1 (MQ=38)                                                                                                      ttAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACgg                                                                                                                                                   >  2:102467/1‑251 (MQ=255)                                                                                                                 aTTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGCTTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGAAATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGCGGTGCTCAAGTGGCTTCTGTTTATATCAGCTGTCCCTACTGTTAAGCTACTGACGGGGTGGTGCGt                                                                                                                                         >  2:466229/1‑250 (MQ=37)                                                                                                                        aaCCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATTACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCGGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaa                                                                                                                                 <  2:498006/251‑1 (MQ=255)                                                                                                                                   aaaaaCCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGGGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGTTGGTGCGTAACGGCAAAAGCAccgcc                                                                                                                       <  2:240080/250‑1 (MQ=255)                                                                                                                                          tttttGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATc                                                                                                                >  1:480272/1‑250 (MQ=255)                                                                                                                                             ttGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCTGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAgc                                                                                                             <  1:438122/250‑1 (MQ=255)                                                                                                                                              tGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAgcg                                                                                                            <  1:182830/250‑1 (MQ=255)                                                                                                                                                         aaTTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGctc                                                                                                >  1:462854/1‑251 (MQ=255)                                                                                                                                                              tataGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCAc                                                                                            <  1:153957/250‑1 (MQ=255)                                                                                                                                                               ataGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTg                                                                                          >  1:12229/1‑251 (MQ=255)                                                                                                                                                                  aaaaTAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTATTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCg                                                                                       <  1:171043/250‑1 (MQ=255)                                                                                                                                                                   aaaTAATGATGCTTTCATAAACCCTGATCTACAAGCACGGTTAGTGATCGGGGATTATAGCCTTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTCTGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCATCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACTGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGt                                                                                      >  2:354388/1‑251 (MQ=255)                                                                                                                                                                   aaaTAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGt                                                                                      <  1:454755/251‑1 (MQ=255)                                                                                                                                                                    aataatGATGATTTCATAAACCCTGATCTACAATAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTCTGCCATTCGAGGTGGTGCTGCAAACTTACTGATTTATAGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTATATCAGCTGCCACTCTTGTTCAGCTACTGACGGGGTGTTGCGTACCGTCAACAGCACCGCCGGATATCAGCGATCTCTCTGCTCTCACTGACGt                                                                                      >  2:4469/1‑250 (MQ=18)                                                                                                                                                                           atgatTTCATAAACCCTGATCTTCAAGAACTGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTTATTTAGTGTATGATGGTGTTTTTGAGGTGTGCCAGTGGCGTCTGTTTCTATCAGCTGTCCCTCCTGTTCATCTACGGACGGGGTGGTGCGTAACGGCAAAATCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATg                                                                              <  1:66920/251‑1 (MQ=18)                                                                                                                                                                                    aaaaCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATAGTTACTTATGACATTAAAGGTGATGCTGCCAACTGACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCa                                                                     <  1:190025/250‑1 (MQ=255)                                                                                                                                                                                           tGATCTACAAGCACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTCTCAGCTGTCCCTCCTGTTCCGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCCCCGCCGGACACCACCGCTACGTCTCCTCTCACTGCCGAAAAACCTGGCAACCGCAGTTCACt                                                              >  2:312836/1‑251 (MQ=18)                                                                                                                                                                                           tGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAACCATGGCAACTGCAGTTAACt                                                              >  2:47704/1‑251 (MQ=255)                                                                                                                                                                                            gATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACt                                                              <  1:215977/250‑1 (MQ=255)                                                                                                                                                                                             aTCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATCTTTACTTATGACATTAAAGGTGCTGCTGCCAACTTACTGATTTAGTGTATGATGGCGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCCGACCTCAGCGCCATCTCTCCTCTCCCTCCCGTAAAACCTGCCACCTGCaccccacttc                                                            >  1:504378/1‑242 (MQ=17)                                                                                                                                                                                                  cAAGACCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACCTTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACAccc                                                       >  2:37615/1‑250 (MQ=255)                                                                                                                                                                                                   aaGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGc                                                      >  2:62642/1‑251 (MQ=255)                                                                                                                                                                                                          gTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATctct                                                                                                    >  1:71147/1‑198 (MQ=255)                                                                                                                                                                                                          gTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATctct                                                                                                    <  2:71147/198‑1 (MQ=255)                                                                                                                                                                                                              gTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTGTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGTGGTGGGGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCGCTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCg                                           <  2:512578/251‑1 (MQ=255)                                                                                                                                                                                                                   cGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTAc                                       <  1:62642/250‑1 (MQ=255)                                                                                                                                                                                                                      ggATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATGTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCAc                                   >  1:212910/1‑251 (MQ=255)                                                                                                                                                                                                                       gATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTCCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGCTCAGCTACTGACGGGGTGGTGCGTAACGGCACCAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGcccc                                  >  1:352493/1‑248 (MQ=255)                                                                                                                                                                                                                               cATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATc                          <  2:124551/251‑1 (MQ=255)                                                                                                                                                                                                                                  ttCAATATTTCCTTATGACATTAAAGGTGCTGCTGCCCACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTCCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCAtt                       >  2:278234/1‑251 (MQ=255)                                                                                                                                                                                                                                    cAATATTTACTTATGACATTAAAGGTGATGCTGCCATCTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCGCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTg                      <  2:270472/250‑1 (MQ=255)                                                                                                                                                                                                                                    cAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTg                      >  1:446155/1‑250 (MQ=40)                                                                                                                                                                                                                                          ttACTTAAGACATTAAAGGAGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGAACGCACCAGAAAATCATTGATATGGc               >  2:349153/1‑251 (MQ=21)                                                                                                                                                                                                                                            aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCTCTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATg                 <  1:50807/247‑1 (MQ=25)                                                                                                                                                                                                                                            aCTTATGACATTAAAGGTGATGCTGACAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTCTCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGACAATCATTGATATgcgat             >  2:50807/1‑247 (MQ=255)                                                                                                                                                                                                                                                 ggACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGTCGGGTTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGGACGCTCCAGAAAATCATGGATAGGGCCATGAAt        <  1:458893/250‑1 (MQ=255)                                                                                                                                                                                                                                                     aTTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCg    >  1:87849/1‑251 (MQ=31)                                                                                                                                                                                                                                                      ttAAAGGTGCTGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAACGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCCACTGCAGTTCACTTCCACCGCTTCTCAACCCGGTACGCACCAGCAAATCATTGATATGGCCATGAATGGCGt   >  2:387325/1‑251 (MQ=9)                                                                                                                                                                                                                                                      ttAAAGGTGATGCGGCCAACTTCCTGATTTAGTGTATGATGGTGTTTTTGAGGAGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTCTCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAGTCATTGATATGGCCATGAATGTCGt   >  2:102754/1‑251 (MQ=255)                                                                                                                                                                                                                                                       tAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGtt  >  2:407588/1‑251 (MQ=31)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                       GAAGTTATGTAATGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGAATGGCGTT  >  NC_000913/3583178‑3583674

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.