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 →	major type 1 subunit fimbrin (pilin)

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2514274 =	NA (NA)	27 (0.800)	27/264	0.4	100%	noncoding (2/1345 nt)	IS186	repeat region
	?	NC_000913	4543608 =	0 (0.000)					coding (494/549 nt)	fimA	major type 1 subunit fimbrin (pilin)
*	?	NC_000913	= 2515616	NA (NA)	25 (0.740)	20/264	0.8	100%	noncoding (1344/1345 nt)	IS186	repeat region
	?	NC_000913	= 4543613	0 (0.000)					coding (499/549 nt)	fimA	major type 1 subunit fimbrin (pilin)

AGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGC  >  NC_000913/2514127‑2514404                                                                                                                                                    |                                                                                                                                   aGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTcc                                                                                                                                   <  3:108581/149‑1 (MQ=255)     aaaGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCAt                                                                                                                                >  2:93578/1‑148 (MQ=255)     aaaGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCAt                                                                                                                                >  8:18949/1‑148 (MQ=255)        ggCCTGAATGAAGAGCAAGGTAACGTGGTTGCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAgc                                                                                                                            <  4:95843/149‑1 (MQ=255)          ccTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAgcgc                                                                                                                          <  4:129081/149‑1 (MQ=255)              aatgaaGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAAc                                                                                                                      >  4:40559/1‑149 (MQ=255)                              aCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                                                                                      <  8:71476/149‑1 (MQ=255)                                 tGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                                                                                   <  7:18949/149‑1 (MQ=255)                                          gCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatat                                                                                          >  4:62004/1‑149 (MQ=255)                                          gCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                                                                            <  2:156100/147‑1 (MQ=255)                                                     aaGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAAc                                                                                  <  3:40559/146‑1 (MQ=255)                                                       gCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATCAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                                                                             <  3:12371/149‑1 (MQ=255)                                                       gCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                                                                             >  7:8704/1‑149 (MQ=255)                                                           gTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                                          <  3:62004/148‑1 (MQ=255)                                                                  gCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct                                                                  <  1:93578/149‑1 (MQ=255)                                                                       aCCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                                             >  2:35291/1‑149 (MQ=255)                                                                                      cTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGTCGATGAATTACTCTCACGATAACTGGTCAGc                                               <  2:36663/148‑1 (MQ=255)                                                                                      cTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGTCGATGAATTACTCTCACGATAACTGGTCAGc                                               >  1:36663/1‑148 (MQ=255)                                                                                           tGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCt                                         >  1:162378/1‑149 (MQ=255)                                                                                              gTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                      <  8:8704/149‑1 (MQ=255)                                                                                              gTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGACGGGGATAATTCCCATAAGCGCTAACTTAACGGTAGTGGTAATACGCCTGATATGATTTAACGTGCCGATGAATTACACTCACGATAACTGGTCAGCAATTCTgga                                      >  3:125088/1‑148 (MQ=21)                                                                                                 aaTCGCAGCACTGGTCCTGTAAGACCATACATAAAAAAGCCGGGGATAAATCCCATAAGCGCTAACTTAAGAGTTGTGGTATTACGCCAGATATGAATAAACGAGCCGATGAATAACTCTAACGATAACTGGTCAGCAATTCTGGCCCa                                   >  6:125795/1‑149 (MQ=18)                                                                                                  aTCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCa                                   >  8:58755/1‑148 (MQ=255)                                                                                                  aTCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCa                                   >  7:53474/1‑148 (MQ=255)                                                                                                     gcagcaCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                       <  4:12824/141‑1 (MQ=255)                                                                                                     gcagcaCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTgg                                       >  3:12824/1‑141 (MQ=255)                                                                                                      cagcaCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                      <  5:17551/141‑1 (MQ=255)                                                                                                      cagcaCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                                      >  6:17551/1‑141 (MQ=255)                                                                                                      cagcaCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                              <  8:53474/149‑1 (MQ=255)                                                                                                          aCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTaa                                                                                                                  >  4:95427/1‑61 (MQ=255)                                                                                                            tGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGc                        >  3:121264/1‑149 (MQ=255)                                                                                                              gTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGccc                      <  7:58755/149‑1 (MQ=255)                                                                                                                  tGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt                            <  1:116471/139‑1 (MQ=255)                                                                                                                  tGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGt                            >  2:116471/1‑139 (MQ=255)                                                                                                                   gTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaag                  >  3:13031/1‑148 (MQ=255)                                                                                                                       gACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTg             <  4:13031/149‑1 (MQ=255)                                                                                                                         ccATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGa           <  1:35291/149‑1 (MQ=255)                                                                                                                             acataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg                              >  8:31723/1‑126 (MQ=255)                                                                                                                                  aaaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGc  <  2:162378/149‑1 (MQ=40)                                                                                                                                                    |                                                                                                                                   AGTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCGGC  >  NC_000913/2514127‑2514404

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