Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,543,608	Δ1 bp :: IS186 (–) +6 bp :: Δ1 bp	coding (494‑499/549 nt)	fimA →	type 1 fimbriae major subunit

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2514274 =	NA (NA)	28 (0.660)	27/480	0.6	87.2%	noncoding (2/1345 nt)	IS186	repeat region
	?	NC_000913	4543608 =	4 (0.100)					coding (494/549 nt)	fimA	type 1 fimbriae major subunit
*	?	NC_000913	= 2515616	NA (NA)	41 (0.960)	37/480	0.2	90.9%	noncoding (1344/1345 nt)	IS186	repeat region
	?	NC_000913	= 4543613	4 (0.100)					coding (499/549 nt)	fimA	type 1 fimbriae major subunit

CAGAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGG  >  NC_000913/2514026‑2514520                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       cAGAAAATTGCCTCCACGCTCTCTCCGCTTGCTGAAGGCATCATTTCTGTGGTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAATCCGGGGATAATTccc                                                                                                                                                                                                                                                      <  2:367216/251‑1 (MQ=255)          gCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCt                                                                                                                                                                                                                                             >  1:167375/1‑251 (MQ=255)              ccACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACt                                                                                                                                                                                                                                         >  2:276783/1‑251 (MQ=255)                                     ggCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                                                                                                                                                                                                  <  1:122663/251‑1 (MQ=255)                                          catcTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                                                                                                                                                                                                             >  1:17211/1‑251 (MQ=255)                                               ctGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                                                                                         >  1:154754/1‑250 (MQ=255)                                               ctGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                                                                                                                                                                                                         >  1:452388/1‑250 (MQ=255)                                                tgtgTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGTGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                                                                                                                                                                                                        >  1:102155/1‑250 (MQ=255)                                                 gtgtTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCGTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGTGTTGTGGTATTACGCCTGATATGATTTa                                                                                                                                                                                                        <  2:461191/249‑1 (MQ=255)                                                     tCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGTGATTATCGCAGGTGCAGTTAAAGGCCTGAATTCAGAGCAAGGTAACGGGGTTTCTAGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCCGCACTGGTGCTGA‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGAATTACGCCTGATATGTTTTAGCGTGc                                                                                                                                                                                                  >  2:274760/1‑251 (MQ=255)                                                             tCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                                                                                                                                          >  2:359480/1‑251 (MQ=255)                                                                 tCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                                                                                                                                                      >  1:281395/1‑251 (MQ=255)                                                                         ttCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACgt                                                                                                                                                                               >  2:366930/1‑249 (MQ=255)                                                                           cgctTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGTGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATaa                                                                                                                                                                            <  2:386410/249‑1 (MQ=255)                                                                             ctgcaagtgggATTATCGCAGTGGCTGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAATCTTGTTTTCGGCTCTTCCCTGGTGAGTGTGCTTGCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCTGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATaa                                                                                                                                                                            <  1:49604/241‑1 (MQ=37)                                                                               tCAATCGGGATTATCGCAGGTGCAGGTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAATCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGTTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACt                                                                                                                                                                          <  2:138767/249‑1 (MQ=255)                                                                                 aaTCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAACAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTggac                                                                                                                                                                      >  2:231052/1‑249 (MQ=255)                                                                                   gcaggattatcgtagcgtccgTTAAAGGCCTGAATGAATAGTAAGGGAACGTCGTTTCTCTCTTCGGTCTTAAGCCTGTTTACGGCTCTACCCTGGTGATTGTGCTTTCTGCTTCAATCGCAGCCTTGGTTC‑GTGAAGACCATACATAATTAAGCCGGGGATAATTCCCATAAGCGCTAACGTAAGGTTTGTGGTATTACGCCTGATATGATTTAACGTGCCTATGAATTTCTCTCCCCATAACTGGTCa                                                                                                                                                                     <  2:42399/230‑1 (MQ=16)                                                                                    cGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAGTACTCTCACGATAACTGGTCAg                                                                                                                                                                    <  2:281395/250‑1 (MQ=255)                                                                                        aTTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAAAGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGACTCTACCCTTGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAACAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAtt                                                                                                                                                               >  1:271846/1‑251 (MQ=255)                                                                                              gCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATATCTGGTCAGCAATTCTGGcc                                                                                                                                                         <  1:117547/251‑1 (MQ=255)                                                                                                aGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGGTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGccc                                                                                                                                                        >  1:70299/1‑250 (MQ=255)                                                                                                 gatGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGAACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCat                                                                                                                                                      >  2:361745/3‑251 (MQ=255)                                                                                                       gTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt                                                                                                                                                <  1:359480/251‑1 (MQ=255)                                                                                                         tAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                                                                                                                              >  1:254289/1‑203 (MQ=255)                                                                                                         tAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTGTCTCGCTGCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                                                                                                                                              <  2:254289/203‑1 (MQ=255)                                                                                                            aGGCCTGAATGAAGAGCAAGGTAACGTAGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGcc                                                                                                                                           <  2:319708/251‑1 (MQ=255)                                                                                                               ccTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACACTGGTGAGTGTGCTGTCTGCGTAAATCGCAGCACTGGTGCTGT‑AAGACCAAACATAAAAAAGCCGGGGATAATTCACATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGACAGAAATACAGGCCCATATTGGTAAACCCg                                                                                                                                         >  2:192681/1‑250 (MQ=37)                                                                                                                      gaaGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTg                                                                                                                                 <  2:395246/251‑1 (MQ=255)                                                                                                                                   aCGTGGGTTCTCGCTTCGTTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATTAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCAc                                                                                                                    <  2:452388/251‑1 (MQ=255)                                                                                                                                      tGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCg                                                                                                                                         >  1:3814/1‑227 (MQ=255)                                                                                                                                      tGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCg                                                                                                                                         <  2:3814/227‑1 (MQ=255)                                                                                                                                            ctcGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCCTATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg                                                                                                           >  2:457647/1‑251 (MQ=255)                                                                                                                                              cGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCg                                                                                                           <  1:166474/249‑1 (MQ=255)                                                                                                                                                cTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGGGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGc                                                                                                       >  1:78475/1‑251 (MQ=255)                                                                                                                                                     gTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTa                                                                                                   >  1:83054/1‑250 (MQ=255)                                                                                                                                                         gAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACccg                                                                                              <  1:147687/251‑1 (MQ=255)                                                                                                                                                                         tctACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGt                                                                              >  1:316708/1‑251 (MQ=255)                                                                                                                                                                               cTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGAtgct                                                                        >  1:373747/1‑251 (MQ=255)                                                                                                                                                                                      gtttgctggcggCGTCAATCGCAGCACTGGTGCTGT‑AAGACCAGACATAAAAAAGGCGGGGATAATTCCCATAAGCGCTAACTTAAGGTTTGTGGTATTACTCCTGATATGATTTAACGTTCCGATGAATTACTCTCACGATATCGGGTCAGCAAGTCTGGCCCATATTGGTAAGCTCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGGTCTAACCCGCCGCCTCGAAATTCGTGAGGCTGCAActc                                                                 <  2:429290/240‑1 (MQ=17)                                                                                                                                                                                                    tCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAACAGCCGGGGATAATTCCCATAAGCGCTAACTTACGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGCTAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGGTGCTGCCACTCTGCTACGTCTgggg                                                   >  2:192680/1‑251 (MQ=255)                                                                                                                                                                                                    tCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGACTgggg                                                   >  2:81782/1‑251 (MQ=255)                                                                                                                                                                                                           cagcaCTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCtt                                            <  1:165690/251‑1 (MQ=255)                                                                                                                                                                                                                      cTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCggg                                 >  1:196283/1‑251 (MQ=255)                                                                                                                                                                                                                         t‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAATGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCgggggg                              >  1:5858/1‑251 (MQ=255)                                                                                                                                                                                                                                      taAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGt                  >  1:406419/1‑251 (MQ=37)                                                                                                                                                                                                                                       aaaataGCCGGGGATAATTCCCAAAAGCGCTAAATTAAGTGTTTTTGTATTTCGTCTGATACTATTTAACGCGCCGATGAATTACTCTCACTATAACTGGCTAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATTCTGCGTCACGTTATGTCGGGGTTCTACCCCGCTGCCGCTTAATTCGTGATGCTGCTACTCTGCTACGTCTGTGGCTGGATTTCGGCCCCGTGGGTACGCCATTACGt                  <  2:470346/250‑1 (MQ=2)                                                                                                                                                                                                                                                      aaTTCCCATAATCGCTAACTTCAGGGTTGTGGTATTACGCCTGTTATGATTTACCGTGCCGATGAATTACTCTCACGATCACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGACGAACTGGCTACTTCGGCACGTAATGCCGGGGCACTAACCCGCCGCCGCGAAATTCGTGGTGCTGCAACTCTGCTACGTCTGGGGCTGGATTACGGCTCCGGGGGGAGGTCATTACGTGAAGCCAcagccttgg  >  2:68118/1‑243 (MQ=11)                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                       CAGAAAATTGCTTCCACGCTCTCTCCGCGTGCTGAAGGCATCATCTCTGTGTTCCTGGTTTCCTTCGCTAACTTCTCTTCAATCGGGATTATCGCAGGTGCAGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGT‑AAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGCACGTAATGCCGGGGCTCTAACCCGCCGCCGCGAAATTCGTGATGCTGCAACTCTGCTACGTCTGGGGCTGGCTTACGGCCCCGGGGGGATGTCATTACGTGAAGTCACTGCATGGG  >  NC_000913/2514026‑2514520

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑
Reads not counted as support for junction
read_name Not counted due to insufficient overlap past the breakpoint.
read_name Not counted due to not crossing MOB target site duplication.