Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	3,720,136	IS3 (–) +3 bp	coding (88‑90/213 nt)	cspA →	cold shock protein CspA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2170171 =	NA (NA)	25 (0.740)	25/484	0.4	96.1%	noncoding (1/1258 nt)	IS3	repeat region
	?	NC_000913	3720136 =	1 (0.030)					coding (88/213 nt)	cspA	cold shock protein CspA
*	?	NC_000913	= 2171428	NA (NA)	35 (1.040)	35/484	0.1	97.2%	noncoding (1258/1258 nt)	IS3	repeat region
	?	NC_000913	= 3720138	1 (0.030)					coding (90/213 nt)	cspA	cold shock protein CspA

ATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAA  >  NC_000913/2169930‑2170390                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                            aTCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCCGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTAcc                                                                                                                                                                                                                    >  1:137663/1‑251 (MQ=255)       aCCGACAAACATTAACTCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCCAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGGCATTAATCCTTCGGCGATAAGCATTGTTGTACTTCCGTCGTCGAGAATAAACGTCTGATATTACCCCGCa                                                                                                                                                                                                               >  2:181932/1‑250 (MQ=255)         cGACAAACATTCAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAAt                                                                                                                                                                                                             >  1:297972/1‑250 (MQ=255)         cGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGAAGACCAGAGGAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTGGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGCCTGATCTTACCCAGAAAt                                                                                                                                                                                                             >  1:749/1‑250 (MQ=255)           aCAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAgg                                                                                                                                                                                                          >  1:80625/1‑250 (MQ=255)           aCAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGt                                                                                                                                                                                                          <  1:216738/251‑1 (MQ=255)                                         aaTGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGAGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTTTTTTAGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGGGATAAGCATTGTCGTACTTTCGCTGTCGGGAATAACGGGCTGATCTTACCCCGCAATAGGGGGCACGCGGCTAATGGAGTaaaacctccg                                                                                                                                                                            >  2:27130/1‑244 (MQ=255)                                                                                          ccGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAAcccc                                                                                                                           <  1:91044/251‑1 (MQ=255)                                                                                                  ccacAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGGAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAg                                                                                                                   <  1:181932/251‑1 (MQ=255)                                                                                                        acgacCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTGAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCg                                                                                                             >  2:107521/1‑251 (MQ=255)                                                                                                            ccAGATTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCCCCGTGATATTTTTGGCGGTCATTAATCCTTCTGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTTAACAGCATTCGCCTg                                                                                                         <  2:34579/251‑1 (MQ=255)                                                                                                                   aaTGTCTTTATTTTCGGAAAGGGAAAACGCTGCCGGGAGACTGTTAGTAATCACCGTGCTATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCAAGCAATAGTGGACCCGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACCAAAACAGTATCAACCAGTAAAAAACCACGTAAACAGCATTCGCCTGAATTTCg                                                                                                  >  2:313860/1‑251 (MQ=25)                                                                                                                          ttATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGa                                                                                            <  1:27130/250‑1 (MQ=255)                                                                                                                                           aaaCGCTGCCGGGAGACTGTTGGGAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGGGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTTACATGCCAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGACTTTCGCAGTGAAGCCCTGAAGCTTGCTg                                                                           >  2:263178/1‑250 (MQ=255)                                                                                                                                                    cGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTCCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCTCCTGAATTTCg                                                                                                  <  2:143340/218‑1 (MQ=255)                                                                                                                                                    cGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTCCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGAGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCg                                                                                                  >  1:143340/1‑218 (MQ=255)                                                                                                                                                          aCTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTAc                                                           <  1:344/251‑1 (MQ=255)                                                                                                                                                                   aaTCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCCTTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCCGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCCCATGACAAAAACAGTAACAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTAAAGCCCTGAAGCTTGCTGAAAGCATCGGGGTTACTGCCGCAGc                                                  >  2:253023/1‑251 (MQ=18)                                                                                                                                                                          gTGATCTTTTTGGCGGTCATTCATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGa                                            >  1:82616/1‑250 (MQ=255)                                                                                                                                                                                            aTTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATca                         <  2:297972/251‑1 (MQ=255)                                                                                                                                                                                            aTTAATCCTTCGGCGATAAGCATTGTAGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACCCAGCCTGTATGAATcc                         >  2:150759/1‑250 (MQ=255)                                                                                                                                                                                                                  ttGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGt    >  2:15776/1‑250 (MQ=255)                                                                                                                                                                                                                   tGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTGAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTa   <  1:107521/250‑1 (MQ=255)                                                                                                                                                                                                                    gTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTaa  >  1:288928/1‑250 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                            ATCAGCACCGACAAACATTAAATCGGCATTAATATCCTGTAATGAACGCTCGGCAATGGAACCATGCATCGAGCGCGTTTTATGGCGGACCGTGCCACCACAGACGACCAGAGTAATGTCTTTATTTTCGGAAAGGGCAAACGCTGCCGGGAGACTGTTGGTAATCACCGTGATATTTTTGGCGGTCATTAATCCTTCGGCGATAAGCATTGTCGTACTTCCGCTGTCGAGAATAACGGTCTGATCTTACCCAGCAATAGTGGACACGCGGCTAAGTGAGTAAACTCTCAGTCAGAGGTGACTCACATGACAAAAACAGTATCAACCAGTAAAAAACCCCGTAAACAGCATTCGCCTGAATTTCGCAGTGAAGCCCTGAAGCTTGCTGAACGCATCGGTGTTACTGCCGCAGCCCGTGAACTCAGCCTGTATGAATCACAACTCTACAACTGGCGCAGTAA  >  NC_000913/2169930‑2170390

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.