Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,755,755	IS5 (+) +4 bp	coding (58‑61/1413 nt)	pykF →	pyruvate kinase I

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 1397238	NA (NA)	24 (1.200)	20/284	0.1	95.9%	noncoding (1195/1195 nt)	IS5	repeat region
	?	NC_000913	1755755 =	1 (0.050)					coding (58/1413 nt)	pykF	pyruvate kinase I
*	?	NC_000913	= 1755758	1 (0.050)	19 (0.950)	18/284	0.2	94.9%	coding (61/1413 nt)	pykF	pyruvate kinase I
	?	NC_000913	= 2102943	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

CTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAA  >  NC_000913/2102812‑2103082                                                                                                                                    |                                                                                                                                            ctTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGCCCCCTGTTCTCTGGCTCCAGATGACAAACATGAACTCATATAAGGGACTAGTTCGCACCTTACTTAATCCCTaa                                                                                                                                  >  1:161573/1‑143 (MQ=11)       tgttttGACGGCGTTTACTGTTGAATTCGCTGTTGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTGTATGGCTCGAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATaaa                                                                                                                            <  1:260728/138‑1 (MQ=14)        ggTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAAt                                                                                                                           <  2:256143/143‑1 (MQ=40)                                       tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCCTATCAGGGACTTGTTCGCACCTTCCTTCATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAAGAAAtt                                                                                            >  1:144167/1‑143 (MQ=255)                                       tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATACATGACACAATAAAATCAGAAAGCACAGGAACtt                                                                                            >  2:181612/1‑143 (MQ=255)                                       tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAAtt                                                                                            >  1:100655/1‑143 (MQ=255)                                       tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAAtt                                                                                            >  2:68553/1‑143 (MQ=255)                                        cGGCGAAGGTAAGTTGATGACTCATGATTAATCCTGTTCTATTGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAAttt                                                                                           <  1:106079/143‑1 (MQ=255)                                          gCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTc                                                                                         <  1:251449/143‑1 (MQ=255)                                            gAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGAAACAAAAAAATCAGAAAGCACAGGAAATTTTCtt                                                                                       >  1:145826/1‑143 (MQ=255)                                               ggTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGACAGCACAGGAAATTTTCTTACg                                                                                    >  2:275297/1‑143 (MQ=255)                                                gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTTTTCGCACCTTCCTTAATCCCTCACATAAATGACACAATAAATTCAGAAAGCCGAGGAAATTTTCTTaaag                                                                                   >  1:126208/1‑140 (MQ=18)                                                   aGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATAAAACAATAAAAACAGAAAACAAAGGAAATTTTCTTACGGaat                                                                                >  1:282846/1‑143 (MQ=18)                                                             tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAACGCACAGGAAATTTTCTTATGGAATAATAATGTaa                                                                      >  1:365350/1‑143 (MQ=255)                                                             tCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAAAAAGCACAGGAAATATTCTTACGGAATAATCATGTaa                                                                      >  1:364927/1‑143 (MQ=255)                                                                   tgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGAttt                                                                >  1:153883/1‑143 (MQ=255)                                                                      aCCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGAAAAAATAATGTAAAGAttacac                                                             >  1:45339/1‑139 (MQ=255)                                                                                        aGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTTATTATTTTCTAACTTTGTTGTCAAAATAAAATCAGAAAGTACAGGAAATTTTCGTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCa                                           <  1:296675/143‑1 (MQ=25)                                                                                         gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAACTGACACTCTAACCTCAGTAAGCAGAAGAAATTTTCTTACGAAATAATCATGTAAAGATTTCGCTTCATATCGTAACACGCag                                          >  1:162623/1‑142 (MQ=25)                                                                                         gATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCACTAACATAAATAAAAAAAAAAAATAAAAAAGCACAGGAAATTTTATTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCAt                                          >  1:26548/1‑143 (MQ=37)                                                                                                  cATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTa                                 >  1:321211/1‑143 (MQ=255)                                                                                                       tctcATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTTTTACGGAATAATCATGAAAAGATTTCGCATCATCTCGAAACAAGCATAAAGTACTAAAAGc                            >  1:33577/1‑143 (MQ=39)                                                                                                       tctcATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAGTGTACTatgtgc                            >  1:370582/1‑138 (MQ=255)                                                                                                             atCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACg                      >  1:297772/1‑143 (MQ=255)                                                                                                                agagACTTGTTCGCACCTTCCTTAATCCCTAAATTTTCTGACACGATAAAATCAGAAAGCACATGAAATTTTCTTACAGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCtt                   <  2:266259/140‑1 (MQ=255)                                                                                                                aGGGACTTGTTCGCATCTTCCTTAATCTCTAACTTAAATGACACAATAAAATCAGAAAGCACAGGAAATGTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCtt                   <  1:271287/143‑1 (MQ=255)                                                                                                                  ggACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAACTCAGAAAGCACAGGAAATTTTCTTACGGAATTATTGTGGATAGATTTCGCTTCATCTAGGAACAGGCATAATGTACTAAAAGCATAACGCTTAc                 >  1:116147/1‑143 (MQ=255)                                                                                                                    aCTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAATTTTTCTTACGGAATAATCATGTAAAGATTTCGATTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACtt               >  1:202081/1‑143 (MQ=255)                                                                                                                            gCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTAAGGAATAAAAATGTAAAGATTTCGATTCATCTAGGAACAAGCATAATGTACTAAAAGCATAACGCTTACTTTCGACata       >  1:256072/1‑143 (MQ=255)                                                                                                                                 ttccttAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAAGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTaa  >  1:192616/1‑143 (MQ=255)                                                                                                                                 ttccttAAGTCCTATCTTAAATGACACAATAAAATCATAAAGTACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTaa  <  2:23809/143‑1 (MQ=255)                                                                                                                                    |                                                                                                                                            CTTTTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAATCCCTAACATAAATGACACAATAAAATCAGAAAGCACAGGAAATTTTCTTACGGAATAATCATGTAAAGATTTCGCTTCATCTCGGAACAGGCATAATGTACTAAAAGCATAACGCTTACTTTCGACATATTTAA  >  NC_000913/2102812‑2103082

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