Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	2,906,143	IS5 (+) +4 bp	coding (433‑436/873 nt)	ygcG →	TPM domain‑containing protein YgcG

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2288919 =	NA (NA)	11 (0.410)	11/490	0.9	100%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	2906143 =	0 (0.000)					coding (433/873 nt)	ygcG	TPM domain‑containing protein YgcG
*	?	NC_000913	= 2906146	0 (0.000)	16 (0.600)	16/490	0.6	100%	coding (436/873 nt)	ygcG	TPM domain‑containing protein YgcG
	?	NC_000913	3130146 =	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

TAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCA‑TCAACTGTCTCACCATA‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCC  >  NC_000913/2288683‑2289164                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      tAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCATCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATTAAACAGCTGGCTTTCTTTTTCCAAGAAAACAAAC‑TCAACTGACTGACCATA‑CTCACAAGATGATGCCAATAATTAGGGTGTATATAAGGCAGGTGCGAATa                                                                                                                                                                                                                                          >  1:122544/1‑250 (MQ=255)    ccGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTCTTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGGTATCTTTTTTCCGGACAAC‑CCACTCCACCGTCTTACCCTAACCAAACA‑TTGTTTCCAATTCTATGGGTGGATATTAGGAAGGGGCGAATAAGCg                                                                                                                                                                                                                                      >  2:100098/1‑251 (MQ=255)         aTTTAATACTGAGTATTCTGCGGCTCAATTAATTTCCATGTACTGATACCCTCCGTTAATTTCCGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGCTTTTCGTATTATTGGCGGGCTCATAACCCAGCTGGGTATTTTTTTCCAGGACAAAACCC‑TCAACTGACTCACCATA‑ATCAACAGATATTGCCAACGCTCTGGGTGAACATAACGAAGGTGCGCATCAGAGGGGa                                                                                                                                                                                                                                  >  2:163292/1‑250 (MQ=255)                 cTGAGTATTCTGCGCCGCCATTCATTTCAATGTACTGATAACATCCGTTAATTTCAGTATTATTACTTACCCCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGCCCTCTACACAAGCTGGATAACTTTTTTCAGGGCAAAACCA‑TCAACCGTCTCACCATA‑CCCAATAGTTGAAGACAATTCTATGGGGGAATATAAGGCAAGTGCCGATAAGCGGGGGAAttctact                                                                                                                                                                                                                         >  1:224845/1‑248 (MQ=255)                                                          ccTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCTGGATATTTTTTTCCCGGACAACACCA‑TCAACTGTCTCCCCCTA‑CTCCACAGTTGATGCCAATGCTATGGGGGAATATAAGGGAGGGGCGGATAAGCGGGGAAATTCTTCTCGGGTGACTCAGTCATTTACTTTCTTCCTGTTTGatgcgtt                                                                                                                                                                                >  2:122718/1‑245 (MQ=255)                                                                                                   tataTGCTGTTCTCCGCCAGGATTTATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCA‑TCAACTGTCTCACCATT‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGTGGAAATTCGTCTCGGCTGTCTCAGTCATTTCATTTCTTCATGTGTTAGCCGAGTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGAcc                                                                                                                                        <  1:155409/250‑1 (MQ=37)                                                                                                        gCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCCTAAACCAGCCGGGTATTTTTTTTCAGGACAAAACCA‑TCAACCGTCTCACCATA‑CTCAACAGTTGATGCCAATGCTATGGGGGAATATAAGGAAGGTGCGAATAAGCGGGGGAATTCTTCTCGGCTGACTCCGTCATTTTATTTCTTTATGTTTGAGCCGGTTTTTTTTCCCGTAAATGCCTTTAATCAGCCTATTTAGACCGTttct                                                                                                                                  >  2:189740/1‑251 (MQ=255)                                                                                                                                tCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCA‑TCAACTGTCTCACCATA‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATcccc                                                                                                           <  1:122718/250‑1 (MQ=255)                                                                                                                                              gTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCA‑TCAACTGTCTCACCATA‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTg                                                                                                >  2:87760/1‑247 (MQ=255)                                                                                                                                                  cATAAACCAGCTGGATATCTTTTTCCAGGACAACACCA‑TCAACTGTCTCACCATA‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGAtctct                                                                                        <  1:189740/251‑1 (MQ=255)                                                                                                                                                           gCTGGATATCTTTTTCCAGGACAACACCA‑TCAACTGTCTCACCATA‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGAc                                                                               <  1:214356/251‑1 (MQ=255)                                                                                                                                                                                    acCA‑TCAACTGTCTCACCATA‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAg                                                      <  1:187046/251‑1 (MQ=255)                                                                                                                                                                                       a‑TCAACTGTCTCACCATA‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGtt                                                   >  2:169828/1‑251 (MQ=255)                                                                                                                                                                                                      cATA‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATa                                     <  1:260431/251‑1 (MQ=255)                                                                                                                                                                                                                  gTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGt                          >  2:7722/1‑251 (MQ=255)                                                                                                                                                                                                                            aTGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGTCTCAGTCATTTCATTTCTTCATGTTTGTGCCGATTGTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTGCTTCGCCATTTAAGGCGTTATGCCCAGTTTTTAGTGAGATCTCTCCCTCTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGt                <  1:163292/251‑1 (MQ=255)                                                                                                                                                                                                                                           ataAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAAcccc  <  1:169828/250‑1 (MQ=11)                                                                                                                                                                                                                                           ataAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAAcccc  <  1:100098/250‑1 (MQ=11)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                      TAACCGTCATTTAATACTGAGTATTCTGCGGCGCCATTCATTTCAATGTACTGATACCCTCCGTTAATTTCAGTATTATTACTTACACCAAATTCTTTTATATGCTGTTCTCCGCCAGGATTGATTTTCGTATTATTGGCGGTCCCATAAACCAGCTGGATATCTTTTTCCAGGACAACACCA‑TCAACTGTCTCACCATA‑CTCAACAGTTGATGCCAATGCTATGGGTGAATATAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGCCAGTTGGTTATCGTTTTTCAGCAACCCC  >  NC_000913/2288683‑2289164

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