Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	243,382	IS1 (+) +9 bp	intergenic (‑79/‑153)	fadE ← / → gmhA	acyl‑CoA dehydrogenase/D‑sedoheptulose 7‑phosphate isomerase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	19796 =	NA (NA)	64 (1.060)	57/482	0.1	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	243382 =	0 (0.000)					intergenic (‑79/‑161)	fadE/gmhA	acyl‑CoA dehydrogenase/D‑sedoheptulose 7‑phosphate isomerase
*	?	NC_000913	= 243390	0 (0.000)	65 (1.080)	59/482	0.1	100%	intergenic (‑87/‑153)	fadE/gmhA	acyl‑CoA dehydrogenase/D‑sedoheptulose 7‑phosphate isomerase
	?	NC_000913	= 258675	NA (NA)					noncoding (1/768 nt)	IS1	repeat region

GAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTG  >  NC_000913/19571‑20043                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         gAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAgt                                                                                                                                                                                                                                 <  1:124123/251‑1 (MQ=255)                                                                              gggAGATAGATTAGTTGTACATTACCACGATTTTGACTCAGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGCATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCTCCGATTTTGAGAACGACAGCGACTTCCGTccc                                                                                                                                                    >  2:347380/1‑251 (MQ=255)                                                                                                                                          agaCATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCTGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGc                                                                                        >  1:170236/1‑251 (MQ=255)                                                                                                                                                                                                        ggATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCCACTTATTGATAGTGTTTTCTGTTCAGATAATGCCCGATGAATTTGTCATGCAACTCCACCGATTTTGAGAACGACAGCGCCTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATCCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGctg                                                                       >  1:376245/1‑204 (MQ=14)                                                                                                                                                                                                        ggATGCTGAAAAGAGAGGTAATTGCTGGTAATGACTCCAACTTATTGAGAGTGTTTTATGTTCAGATAATGCCCGATGTCTTTGTCATGCAACCCCACCGATGTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCGGCCTCAGATCCAGGTCATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACTTGCat                                                                       <  2:376245/206‑2 (MQ=14)                                                                                                                                                                                                        acaagcgaatggcgttcggaagagatGGTAAGGACTCCAACTTAGTAATAGTGTTTTATGTTCAGATAAGGCCCGATTACTTTGTCATGCAGCGCCACCGATGTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATGCGCTGCGTATATCGCGTGCTGATTACGTGCAGCTTTCCCGTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCAt                         <  2:327428/226‑1 (MQ=11)                                                                                                                                                                                                         caagcgaatgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATc                        >  1:314102/25‑251 (MQ=2)                                                                                                                                                                                                            gcgaatgtctttctgaagagatGGTACTGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCat                     >  1:590497/22‑251 (MQ=255)                                                                                                                                                                                                              gaatgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGTTTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGAGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCTGGAGTCCT‑ACAGCGGCCAGCCATCCGTCATCCatat                   >  2:193308/20‑251 (MQ=11)                                                                                                                                                                                                              gaatgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTAATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCACTTCAGGCGGGATTCAT‑ACAGCGGCCAGACATCCGTCATCCatct                   >  1:255358/20‑249 (MQ=11)                                                                                                                                                                                                               atatatagtccctctcgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATc                  >  1:163712/19‑251 (MQ=2)                                                                                                                                                                                                               aatgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATc                  >  2:390538/19‑251 (MQ=2)                                                                                                                                                                                                                aaaaGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCTCCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATca                 <  1:347380/251‑1 (MQ=38)                                                                                                                                                                                                                 tgtctttctgaagagatGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATcac                >  1:126632/17‑251 (MQ=2)                                                                                                                                                                                                                     agtccctctcgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGTTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACg            <  1:27839/239‑1 (MQ=2)                                                                                                                                                                                                                         cctctcgttGGTAATGACTCCAACTTCTTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCCCAATTCGATGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTAAGGCGGGATTCATTAC‑GCGGCCAGCACTCCGTCCTCCATATCCCCCCGTCa         >  2:544541/9‑250 (MQ=11)                                                                                                                                                                                                                            ctcgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAgg     >  1:455342/6‑251 (MQ=11)                                                                                                                                                                                                                              cgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGcc                                                                                                                                                >  1:401685/4‑111 (MQ=11)                                                                                                                                                                                                                              cgttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCGCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGcc                                                                                                                                                <  2:401685/108‑1 (MQ=2)                                                                                                                                                                                                                               gttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGgcca                                     <  2:281558/214‑1 (MQ=11)                                                                                                                                                                                                                               gttGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGgcca                                     >  1:281558/3‑216 (MQ=11)                                                                                                                                                                                                                                atGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTg  >  1:135542/2‑250 (MQ=11)                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                         GAATCTGCAATACGGATTATTCTGCGCTTTTTAGTCCAGCGGTGCGTTAATCGGCAGCTCCCCCAAAGTTAAGGTGGGGGAGATAGATTAGTTGTACATTACCACGATTTTGACTCGGCTCATTATTTGCCCGCTTGAGACATTGTTTCCATATGTACGCGGGCGAATAAATAGAGGAATCTGATTACTTCCTTCATGGGGATGCTGAAAAGAGTAGTAATTGCTGGTAATGACTCCAACTTATTGATAGTGTTTTATGTTCAGATAATGCCCGATGACTTTGTCATGCAGCTCCACCGATTTTGAGAACGACAGCGACTTCCGTCCCAGCCGTGCCAGGTGCTGCCTCAGATTCAGGTTATGCCGCTCAATTCGCTGCGTATATCGCTTGCTGATTACGTGCAGCTTTCCCTTCAGGCGGGATTCAT‑ACAGCGGCCAGCCATCCGTCATCCATATCACCACGTCAAAGGGTG  >  NC_000913/19571‑20043

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