Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,630,598	IS5 (+) +4 bp	coding (313‑316/1461 nt)	ydfI ←	putative NAD‑dependent D‑mannonate oxidoreductase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1630598 =	0 (0.000)	30 (0.680)	27/304	0.6	100%	coding (316/1461 nt)	ydfI	putative NAD‑dependent D‑mannonate oxidoreductase
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region
*	?	NC_000913	= 1630601	0 (0.000)	30 (0.680)	28/302	0.5	100%	coding (313/1461 nt)	ydfI	putative NAD‑dependent D‑mannonate oxidoreductase
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region

TATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGC  >  NC_000913/2101537‑2101966                                                                                                                                                                                                                     |                                                                                                                                                                                                                          tATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAg                                                                                                                                                                                                            <  1:588671/228‑1 (MQ=255)    cGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCCGAATAGATACTTTCTTTAGGAAGTTTTGTTTTATTAAGGAAGGTGCGACTaaac                                                                                                                                                                                                           >  2:651615/1‑224 (MQ=255)    cGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGc                                                                                                                                                                                                           >  2:242862/1‑226 (MQ=255)                           gAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACt                                                                                                                                                                                    <  1:243236/226‑1 (MQ=255)                                         gccagGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCatttcatttc                                                                                                                                                                     >  2:397736/1‑227 (MQ=255)                                                                                            ccTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTaagg                                                                                                                                                                                                                          >  2:161928/1‑123 (MQ=255)                                                                                             cTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTaagga                                                                                                                                                                                                                         <  1:162163/123‑1 (MQ=255)                                                                                                                                                         cgcACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAttctt                                                                                                                                                                                               >  1:169179/1‑89 (MQ=255)                                                                                                                                                          gcACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAAttcttc                                                                                                                                                                                              <  2:168925/89‑1 (MQ=255)                                                                                                                                                              cAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTa                                                >  1:114945/1‑227 (MQ=255)                                                                                                                                                                           tATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAATGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGcc                                  >  2:411864/1‑228 (MQ=255)                                                                                                                                                                               aGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGaa                              <  1:398391/228‑1 (MQ=255)                                                                                                                                                                                   tAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAg                                                                      >  2:225509/1‑184 (MQ=255)                                                                                                                                                                                    aGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGt                                                                     <  1:225854/184‑1 (MQ=255)                                                                                                                                                                                           ttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCg                                                                                                                                                                                           >  2:390176/1‑59 (MQ=255)                                                                                                                                                                                           ttctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTgg                                        >  2:423030/1‑206 (MQ=255)                                                                                                                                                                                            tctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCgg                                                                                                                                                                                          <  1:390820/59‑1 (MQ=255)                                                                                                                                                                                            tctttAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGt                                       <  1:423721/206‑1 (MQ=255)                                                                                                                                                                                                           gTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGc  <  2:114773/228‑1 (MQ=33)                                                                                                                                                                                                                     |                                                                                                                                                                                                                          TATCGCATGAAAAGCGGGAACATAATGAAGTCTGACACCAGCCAGGCTAAGCCTCAAGCACAGGTCAATATCTTCACAGTACATAAAGTAACCTTGATCGAAGCCATTTACACGCACAAAATCTGAAAAACGTACCAGCATAAATGATCCTGCGCACCAATCAACAACCGTATCAGAATAGATACTTTCTTTAGGAATTTTTGTTTTATTAAGGAAGGTGCGAATAAGCGGGGAAATTCTTCTCGGCTGACTCAGTCATTTCATTTCTTCATGTTTGAGCCGATTTTTTCTCCCGTAAATGCCTTGAATCAGCCTATTTAGACCGTTTCTTCGCCATTTAAGGCGTTATCCCCAGTTTTTAGTGAGATCTCTCCCACTGACGTATCATTTGGTCCGCCCGAAACAGGTTGGCCAGCGTGAATAACATCGC  >  NC_000913/2101537‑2101966

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