Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,105,376	IS2 (–) +5 bp	intergenic (+83/‑34)	csgC → / → ymdA	curlin chaperone/uncharacterized protein YmdA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	381260 =	NA (NA)	24 (0.590)	24/484	0.5	100%	noncoding (1/1331 nt)	IS2	repeat region
	?	NC_000913	1105376 =	0 (0.000)					intergenic (+83/‑38)	csgC/ymdA	curlin chaperone/uncharacterized protein YmdA
*	?	NC_000913	= 1105380	0 (0.000)	9 (0.220)	9/486	1.7	100%	intergenic (+87/‑34)	csgC/ymdA	curlin chaperone/uncharacterized protein YmdA
	?	NC_000913	1467910 =	NA (NA)					noncoding (1331/1331 nt)	IS2	repeat region

CGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCC  >  NC_000913/1467670‑1468136                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                   cgcTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGccc                                                                                                                                                                                                                          >  1:405851/1‑251 (MQ=255)                                  gcACAGGATATCGTTATACCCGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCCGGGGCGCTACGCAAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAAccc                                                                                                                                                                                         >  2:311339/1‑251 (MQ=255)                                   cacaGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAAccc                                                                                                                                                                                         >  1:15328/1‑250 (MQ=255)                                               ttATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGc                                                                                                                                                                            >  1:128751/1‑251 (MQ=255)                                                       gCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGc                                                                                                                                                                      >  1:264234/1‑249 (MQ=255)                                                             ttGTCGTTAGTGGTGAGTCGTCTGTATATCTGAATGGAGATACGACCATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGATATTCGGGTACAGGATAAAGCTTTGCTGGCAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCCTGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTc                                                                                                                                                              >  1:281326/1‑251 (MQ=255)                                                                                  cgGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCa                                                                                                                                          <  2:128751/250‑1 (MQ=255)                                                                                   ggTAGATCTGAATGGAGATACGACAATCAGCGGAGAATTACTTCTGGGGTTTGCCGGTGTTATTCGGGTACAGGATAAATCGTTGCTTGTAATTGGCAGTGGCGCTACGCTAACAATGCAGGATAGTGACAGTTTCGAACATCATGGGACAAGAACCCAGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAgc                                                                                                                                        <  1:415313/251‑1 (MQ=25)                                                                                         tCTGACTGGAGATACGACAATCAGCGGAGGATTCCCTGTGGGGTTTGACGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCCGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACt                                                                                                                                   >  2:254477/1‑250 (MQ=255)                                                                                         tCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAAACCAGCGCACt                                                                                                                                   >  1:247485/1‑250 (MQ=255)                                                                                              aTGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATGCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGGGGCTCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATGACAAGCCCGc                                                                                                                                                                            <  1:384106/204‑1 (MQ=255)                                                                                              aTGGAGATAAGACACTCCGCGGCAAATGCCCTCTGGGGTTTGCCGGGGTTATTAGGGTCCAGGATAACGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACCATGCAGGATATTGACAGTTTTGAACATCATGGGACCAGACCCCTGGATTTGCCCATATATTTCCAGACATCTGTTCTCACTTACACCATTACAAGCCCGc                                                                                                                                                                            >  2:384106/1‑204 (MQ=255)                                                                                                      aCGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTTGCGAGCGATACCCCAGCGCACTATGCGGATGCCAtt                                                                                                                     <  2:73283/251‑1 (MQ=255)                                                                                                      aCGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCAtt                                                                                                                     >  2:15278/1‑251 (MQ=255)                                                                                                        gACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTCCAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCg                                                                                                                   >  2:5604/1‑251 (MQ=255)                                                                                                                              gtGGGGTTTGGCTGGGTTATTCGGGTACAGGATAAAGCTTGGCTGGAAATTGGCAGTGGCGCTGCGCTAACAATGCAGGATAGTGACAGTTTTGAACATCATGGTACAAGATCCCTGGATTTGCGCCTGTAGTTGCAGGCATCTGTTATCACTTAGCCGATTGCGAGCCCGCTGCGGCTGATATTCCCGTGGCGAGCGTTAACCTAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCctct                                                                                              <  1:180342/249‑1 (MQ=21)                                                                                                                                   gTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTGGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCTCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAAGGCTCGAACGCCTCTGCAAgg                                                                                        <  2:207132/251‑1 (MQ=255)                                                                                                                                                           ggATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCt                                                                <  2:405851/251‑1 (MQ=255)                                                                                                                                                                aaGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGttt                                                           >  2:157010/1‑251 (MQ=255)                                                                                                                                                                              aTTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGCATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCCTTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGAACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTATGGGCATGATACTg                                              >  1:381631/1‑250 (MQ=255)                                                                                                                                                                                  gcggTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGCACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCTCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGACGGCCTCTGCAAGGTTCTTTGCTGCCGTTACCCCGTCTGGTTTTGGCATGATCCTGCTGt                                          >  2:340141/4‑250 (MQ=17)                                                                                                                                                                                          gcgacgctaACAATGAAGGATATTGACAGGTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCCTTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCCAGGATCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCt                                 >  2:17225/4‑251 (MQ=255)                                                                                                                                                                                                  aaCAATGCAGGCTATTGACAGCTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGtt                         >  2:247950/1‑251 (MQ=255)                                                                                                                                                                                                           ggATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGa                  >  2:137505/1‑249 (MQ=255)                                                                                                                                                                                                                          tGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTcc  <  1:409604/250‑1 (MQ=40)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                   CGCTTGCAGATGGGTCGCAGGTTACTTCCTTGCGCACAGGATATAGTTATACCAGCGTTATTGTCGTTAGTGGTGAGTCGTCGGTATATCTGAATGGAGATACGACAATCAGCGGAGAATTCCCTCTGGGGTTTGCCGGGGTTATTCGGGTACAGGATAAAGCTTTGCTGGAAATTGGCAGTGGCGCTACGCTAACAATGCAGGATATTGACAGTTTTGAACATCATGGGACAAGAACCCTGGATTTGCCCCTATATTTCCAGACATCTGTTATCACTTAACCCATTACAAGCCCGCTGCCGCAGATATTCCCGTGGCGAGCGATAACCCAGCGCACTATGCGGATGCCATTCGTTATAATGCTCGAACGCCTCTGCAAGGTTCTTTGCTGCCGTTAACCCGTCTGGTTTGGGCATGATACTGATGTAGTCACGCTTTATCGTTTTCACGAAGCTCTCTGCTATTCC  >  NC_000913/1467670‑1468136

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.