Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,366,110	IS1 (+) +8 bp	coding (657‑664/1302 nt)	dcuA ←	C4‑dicarboxylate transporter DcuA

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	3583428 =	NA (NA)	11 (0.870)	9/310	0.2	100%	noncoding (1/768 nt)	IS1	repeat region
	?	NC_000913	= 4366117	0 (0.000)					coding (657/1302 nt)	dcuA	C4‑dicarboxylate transporter DcuA
*	?	NC_000913	= 3584195	NA (NA)	7 (0.560)	6/308	0.5	100%	noncoding (768/768 nt)	IS1	repeat region
	?	NC_000913	4366110 =	0 (0.000)					coding (664/1302 nt)	dcuA	C4‑dicarboxylate transporter DcuA

CAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCGACTTTACAGATAACGGTGGTCAGAGTTGCAACGCTGAGCATGATGATCAGGATTGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGAAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATG  >  NC_000913/4365884‑4366346                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      gatgacggatggctggcctctgtatgaagcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagctaattgagcggcataacctgaatctgaggcagcacctggcacggctgtgacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaatcactatcaatgaggtggagtcattaccCGGATTTGATTTCAATCTGCTTTTc                                                                                                                                                                                                                      <  1:146026‑M2/25‑1 (MQ=255)        gatggctggctgctgtatgaatcccgcctgaagggaaagctgcacgtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgatgcagcacctggcacggctgggacggaaggcgctgtcgttcgcaaaatcggtggatctgcatgacaaagtcatcgggcattatctgaacatataacactagcaagaagttggagtcattcccCGGATTTGATTTCAATCTGCTTTTCACCgcgc                                                                                                                                                                                                               <  1:263879‑M2/32‑1 (MQ=255)                                 gTTGAGGATGTTGTCGGTATCGACTTTACAGATAACGGTGGTCAGAGGTGCAACGCTGAGCATGATGATCAGGATTGCGTTGGGGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGgg                                                                                                                                                                                                                                     >  1:270619‑M2/1‑202 (MQ=255)                                           ttgtCGGTATCGACTTTACAGATAACGGTGGTCAGAGTTGCAACGCTGAGCATGATGATCAGGATTGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtg                                                                                                                                                                             >  1:101111‑M2/1‑192 (MQ=255)                                                     gtaatcagcaagcgatatacgcagcgaattgagcggcataacctgaatctgaggcagcacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGGATTTGATTTCAATCGGCTTGTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCg                                                                                                                                                                  <  1:204368‑M2/77‑1 (MQ=255)                                                                        ngtcAGAGTTGCAACGCTGAGCATGATGATCAGGATTGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttct                                                                                                                                                                 <  1:56603‑M2/232‑71 (MQ=255)                                                                                                 atgatCAGGATTGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACATCCAGACGTACGTTTTTGCACCGGATTTGggtgatgctgccaacttactgatttagtggatgatggtgttgttgaggtgctccagtggcttctgtttctatcagctttccctcctgttcagctactgacgggttggtgcgta                                                                                                                      >  1:239834‑M2/1‑138 (MQ=255)                                                                                                   ntctgaggcagcacctggcacggctgggacggaagtcgcagtcgttctcaaaatcgaagaagctgcaagacaaagtcatagggccttatctgaacataaaacactatcaataagttggagacattacaCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAAc                                                                                                                                                                                                       >  1:32767‑M2/129‑168 (MQ=255)                                                                                                              cacctggcacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGAAGGACATCACCAGAACCGcc                                                                                                         >  1:43833‑M2/118‑251 (MQ=255)                                                                                                              cGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtgtctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcac                                                                                                         >  1:212623‑M2/1‑125 (MQ=255)                                                                                                                     cacggctgggacggaagtcgctgtcgttctcaaaatcggtggagctgcatgacaaagtcatcgggcattatctgaacataaaacactatcaataagttggagtcattaccCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCa                                                                                                                                >  1:242157‑M2/111‑221 (MQ=255)                                                                                                                                     ttttCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGggtgatgctgc                                                                                                                                                                                                                            >  1:128887‑M2/1‑102 (MQ=255)                                                                                                                                                      gCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaa                                                                  <  1:100464‑M2/250‑166 (MQ=255)                                                                                                                                                                    tgatTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGgg                                                                                                                                                                                                                                     >  1:43825‑M2/1‑71 (MQ=255)                                                                                                                                                                     gatTGCATAGATAACCAAGCCAACTACGACCAGAAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggc                                                                                                                 <  1:91845‑M2/188‑119 (MQ=255)                                                                                                                                                                              gATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaact                                                      <  1:18237‑M2/238‑178 (MQ=255)                                                                                                                                                                                                                ataagtgggagtcattaccCGGATTTGATTTCAATCTGCTTTTCCCCGCGCAGTACAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTAAACAGCATAGTGACCAGAAAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTgg       <  1:235744‑M2/232‑1 (MQ=255)                                                                                                                                                                                                                      gacgTTTTTGCACCGGATTTGggtgatgctgccaacttactgatttagtgtatgatggtgtttttgaggtgctccagtggcttctgtttctatcagctgtccctcctgttcagctactgacggggtggtgcgtaacggcaaaagcaccgccggacatcagcgctatctctgctctcactgccgtaaaacatggcaactgcagttcacttacaccgcttctcaacccggtacgcaccagaaaatcattgat  <  1:125028‑M2/250‑230 (MQ=255)                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                      CAGGCGCTCATACCTGCTTTGAAGGTGCTGGAGTTGAGGATGTTGTCGGTATCGACTTTACAGATAACGGTGGTCAGAGTTGCAACGCTGAGCATGATGATCAGGATTGCGTTGGTGGTGTTCATCAGCGGTTTTTCAACCAGACCCATGCTTGGGCTGTTGATGATTGCATAGATAACCACGCCAACTACGCCCAGCAGGAACAGCCAGACGGACGTTTTTGCACCGGATTTGATTTCAATCTGCTTTTCACCGCGCAGTTCAACCAGGCCCTCTTCCAGACGCTTGCGATAAATCGGATCGTCAGAGAGTTTGGAGTTGAACAGCATAGTGACCAGGAAGGACATCACCAGAACCGCCAGCAGGGTGGACGGGATGACCACGGAGAGCAGATGGAGGTAGCTGATGCCATGACCTTCCATCACGGAAGACATGTAAACCACTGCCGCTGAGATTGGCGATG  >  NC_000913/4365884‑4366346

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