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)	22 (0.750)	22/488	0.3	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	4541425 =	0 (0.000)					coding (469/603 nt)	fimB	Type 1 fimbriae regulatory protein FimB
*	?	NC_000913	= 3584195	NA (NA)	36 (1.230)	35/488	0.1	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	= 4541432	0 (0.000)					coding (476/603 nt)	fimB	Type 1 fimbriae regulatory protein FimB

AGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGA  >  NC_000913/3583197‑3583666                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               aGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCAAACTTACt                                                                                                                                                                                                                             >  1:138809/1‑251 (MQ=255)              aGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTatg                                                                                                                                                                                                                >  1:140720/1‑251 (MQ=255)                      ccTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGtt                                                                                                                                                                                                        <  1:223000/251‑1 (MQ=255)                         gATTATCTTACTTTTTTAAGAGCAATGAATGGATTCGAAGTTCATGGATTACGATTATTTAGCCTCACGATTCCAGAACCTTCAGTTAACAACCTTTTTGCCGTAAATGAATTTTATAGAAATATTGATGATTTCTTAAACCCTGATATACAAATAAGGTTAGAAATCGGGGATTATAGAATTTCAATATTTCCTTATGACATTAACAGTTATGATGCCAACATACTGATTTatagtatgtgtgaggtttt                                                                                                                                                                                                     >  2:210866/1‑233 (MQ=255)                             atttttcatttttAAGAGCAAGGGATCGATTCGAATTTAATGGATTACGATTATTTAACCTCTCGTTTCCTAAACATTCAGTTAAAAACCTGTTTGCCGTAAATGAATTTGATAGAAATAATGATGATTTCATAATACCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTTAATATTTACTTATGACATTAAAGGTGATGCTGCGAACTTACTGCTTTAGTGTATGATGTTGTTTTTGAg                                                                                                                                                                                                  <  2:271767/242‑1 (MQ=255)                              tCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAACCCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTACAGGAGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGt                                                                                                                                                                                                >  1:239021/1‑251 (MQ=255)                                ttagttttttAAGAGCAATGGTAGGATTAGAAGGTCAGGAATTACAATTATTTAGCGTCTCGATTCAAGATCCTTCAGTTAAAAAACTTTTTGCCGTAAATGAATTTTATAGATATAAATAAGATTACATACACTATGATCTACAAGAACTGTTATTGATCGGGGCTTATATCATTTCAATATTTACTTATGACATTAAAGGTAATGCTGCAAACTTACTGATTTAGTGTAGTATGGTATTATTGAGGta                                                                                                                                                                                               >  2:3495/5‑249 (MQ=255)                                                         atTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAACCCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACACTGATCTACAAGAAAGGTTAGTGATCGGGGATTATATCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGACAAATTACTGATTTAGTGTATGATGGTGTTTTTGGGGTGATCCCGTGGCTTCTGTTTCTATCAGc                                                                                                                                                                     >  2:24384/1‑251 (MQ=255)                                                                 ttAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCcctc                                                                                                                                                             >  1:66368/1‑251 (MQ=255)                                                                                                 cacgaaactacAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCTGTGATTATAGCATTTCAATATTTTCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGATACTGACGGGGTGGTGCGTAAc                                                                                                                              <  2:198188/240‑1 (MQ=25)                                                                                                        cTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCaaa                                                                                                                        >  1:170541/1‑249 (MQ=255)                                                                                                        cTTCAGTTAAAAACCTTTATGCCGGAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAATAACGGTTAGTGCTCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGCGATGCTGCCAACTTACTGATTTAGTGTATGGTGGTGTTTTTGAGGTGCTCCAGTGACTTCTGTTGCTATGAGCTGTCCCTCCTGTTCAGCTCCTGCCGGGGTAGTGCGTTACGGccaa                                                                                                                        >  2:275276/1‑246 (MQ=25)                                                                                                         ttCAGTTAAAAACCTTATAGACGTAAATGACTTTTATAGAAATAATGATGATTTCATACACCCTGATCTACAAGAACGGTTAGTGATCGGGTATTAAAGCATTTCAATATTTACTTATGACATTAAAGGCGATGCTGCAAACTTAATGATTTAGAGTGTGATGGTGTTTTAGAGTTGCTCCAGTGGCTTATGTTTCAATCAGCAGTCCCTCATGATCAGCTACTGTCGGGGTGGTGCGTAActtcaaagt                                                                                                                      >  2:244776/1‑242 (MQ=255)                                                                                                             gTTAAAAACCTTTTTGCCGTCAATGACTTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCCATATTTACTTATGAAATTAAAGTAGATCCTGCCATCTTACTGATTTAGTGTATGATGGTGTTTTTGCGGTGCTCCAGTGGATTCTGTTTCTCTCAGCTGTACCTCCTGTTCCGCTACTGACGGGGCGGCGCGTCACGGCAAAAGccccc                                                                                                                 >  2:367029/1‑247 (MQ=255)                                                                                                                    aCCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACa                                                                                                          <  2:140720/251‑1 (MQ=255)                                                                                                                                tAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTAtct                                                                                              >  1:360008/1‑251 (MQ=255)                                                                                                                                                      tgatgaTTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAAc                                                                         >  1:368035/1‑250 (MQ=255)                                                                                                                                                       gatgatTTCATACACCCTGATCTACAAGAACGGTTAGTGATCGGGGCTTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                                                                                                                                                                         >  2:108924/1‑105 (MQ=255)                                                                                                                                                       gatgatTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGAtt                                                                                                                                                                                                                         <  1:108924/105‑1 (MQ=255)                                                                                                                                                                       ttGATATACAAGAACGGTTAGTGAGCGAGGATTATAGCATTTCACTATTTACTAATTAAATTAAAGGTGATGCTGCCAACTTACTGTTTTAGAGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTATGGTTATATCAGATGTCCCTCCTGTTGAGCTACTGACGGGGAGGTGAGTAACGGCAAAAGCACCCCCGGACATCAGCGCTATCGCTGCTCTCACTACCGAAAAAAATgggaaatgtagttaac                                                       >  2:259406/2‑237 (MQ=9)                                                                                                                                                                       cTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCAc                                                       >  1:91963/1‑251 (MQ=255)                                                                                                                                                                                   aaCGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTtc                                           <  1:244724/251‑1 (MQ=255)                                                                                                                                                                                            ggATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCgg                                  <  1:3366/250‑1 (MQ=255)                                                                                                                                                                                                cGGGGATTATAGCAGTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACGGACGGGGTGGTGCTTAACGGCAAAAGCACCGCCGGACATCAGCGCTTTCTCTGCTCTCACTGCCGTACAACATGGCAAATAAAGTTAACTTACACCGCTTCTCAACCCggaccg                              >  2:243856/1‑247 (MQ=14)                                                                                                                                                                                                                tCAATATTTACTTATGACATTAAAGGCGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCCCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCAtt               <  1:275276/250‑1 (MQ=255)                                                                                                                                                                                                                      tttACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATgg        >  1:175798/1‑251 (MQ=37)                                                                                                                                                                                                                         aCTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGcc      >  1:81178/1‑250 (MQ=35)                                                                                                                                                                                                                         aCTTATGACATTAAAGGTGATGCTGCAAACTTACTGATTTAGTGTATGATGGTGTTTTTGGGGTGCTCCAGTGGCGTGTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGGGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGcc      <  2:91963/250‑1 (MQ=255)                                                                                                                                                                                                                            tATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATg   <  1:24384/250‑1 (MQ=34)                                                                                                                                                                                                                             aTGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGa  >  1:145277/1‑250 (MQ=34)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               AGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAAACCCTGATCTACAAGAACGGTTAGTGATCGGGGATTATAGCATTTCAATATTTACTTATGACATTAAAGGTGATGCTGCCAACTTACTGATTTAGTGTATGATGGTGTTTTTGAGGTGCTCCAGTGGCTTCTGTTTCTATCAGCTGTCCCTCCTGTTCAGCTACTGACGGGGTGGTGCGTAACGGCAAAAGCACCGCCGGACATCAGCGCTATCTCTGCTCTCACTGCCGTAAAACATGGCAACTGCAGTTCACTTACACCGCTTCTCAACCCGGTACGCACCAGAAAATCATTGATATGGCCATGA  >  NC_000913/3583197‑3583666

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