Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,631,612	IS5 (+) +4 bp	coding (976‑979/1623 nt)	yhiJ ←	DUF4049 domain‑containing protein YhiJ

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	27 (0.970)	26/488	0.2	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	3631612 =	0 (0.000)					coding (979/1623 nt)	yhiJ	DUF4049 domain‑containing protein YhiJ
*	?	NC_000913	= 2290113	NA (NA)	33 (1.190)	33/488	0.1	100%	noncoding (1/1195 nt)	IS5	repeat region
	?	NC_000913	= 3631615	0 (0.000)					coding (976/1623 nt)	yhiJ	DUF4049 domain‑containing protein YhiJ

ACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTT  >  NC_000913/2288674‑2289153                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            aCAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATCCCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCGTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCGTTTTCCATGACAACACCAGCAACTGTCTCTCCATACTCAACAGTTGATGCCAATGCTATGGGTTAATATaag‑gaagg                                                                                                                                                                                                                                       <  2:44003/251‑1 (MQ=255)                                                cattTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGAGAAATTCTTCTCGGCTGACTCAGtcatttgat                                                                                                                                                                                         <  2:20449/250‑4 (MQ=255)                                                         gTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTTCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGt                                                                                                                                                                               >  1:144164/1‑251 (MQ=255)                                                              gATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCACATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCAACATACTCAACAGTTGCTGCCAATGCTATGGGGGAATATAAG‑GAAGGGGCGAATAAGCGGGGAAATTCTTCTCGGCTTGCTCAGACATTTCCTTTCTTTCTGTTTg                                                                                                                                                                            >  1:143091/1‑249 (MQ=255)                                                                tACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATCTTTTTCCAGGACAACAACATCAACTGTCTCACCATACTCAAAAGTTGATGCCAAAGCTATGGGGGGATATAAG‑GAAGGTGCGAATAAGCGGGGGAATTCTTCTCGGCTGGATCAGTCATTTCCTTTTTTCATGTTTGTGcc                                                                                                                                                                        >  1:32391/1‑251 (MQ=255)                                                                 aCCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGcc                                                                                                                                                                        <  2:32391/250‑1 (MQ=255)                                                                                      tattaCTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGc                                                                                                                                                   <  1:8683/250‑1 (MQ=255)                                                                                                                                        ttCGTATTATTGGCGGTCCCCTAAACCAACTGGATATCTTTTTCCAGGGCAACACCATCAACTGGCTCAACATACTCAAAAGGTGATGCCAATGCTATTGGTGAATATAAG‑GAAGGTGCGAATAAGCCGGGAAATTCTTCTCGGCTGACTCCGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc                                                                                                >  1:73153/1‑251 (MQ=18)                                                                                                                                                  tGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttt                                                                                                                                                                  <  1:29298/175‑1 (MQ=255)                                                                                                                                                  tGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGAtttt                                                                                                                                                                  >  2:29298/1‑175 (MQ=255)                                                                                                                                                                       ggATATCTTTTTCCAGGACAACACCATCAACTGGCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTGGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTCTTTAGTGAGAGCTCTCCCACTGACGTa                                                                 <  2:2706/251‑1 (MQ=255)                                                                                                                                                                            tCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCAtt                                                            >  1:239369/1‑251 (MQ=255)                                                                                                                                                                               ttttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCAttt                                                           <  2:100605/249‑1 (MQ=255)                                                                                                                                                                                 ttCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGACATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTc                                                       >  2:163362/1‑251 (MQ=255)                                                                                                                                                                                  tCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGCGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTcc                                                      >  2:82136/1‑251 (MQ=255)                                                                                                                                                                                                                    cAACAGTTGATGTCAATGCTGTGGGTGCATATAAG‑GAAGGTGCGAATAAGCGCGGAAATTCTTCTCGGCTGACTCCCTCATTTCATTTCTTCATGTTTGAGCCGCTTTTTTCTCCCGTAGATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTCTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCg                    >  2:32452/1‑251 (MQ=11)                                                                                                                                                                                                                                       aatgttcatgtataAGAGAAGGTGAGAATAACCGGGGAACTTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATATTTGAGCCGATTTTTTTTCCCGTTAATCCCTTGAATCAACCTATTTAGACCGTTTCCTACCCCTTTAAGGCGTTATCCCCAGTTTTTACTCAGATCTCTCACACTGACGTATCATTTTGTCCGCACCAAACAGGTTGGCCCGCGTGAATAACATCCCCCGTTGGTTATCGgttt  >  2:141650/11‑251 (MQ=2)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                            ACAATTTGATAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCATCAACTGTCTCACCATACTCAACAGTTGATGCCAATGCTATGGGTGAATATAAG‑GAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTT  >  NC_000913/2288674‑2289153

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