Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	4,619,685	Δ1 bp :: IS186 (+) +6 bp :: Δ1 bp	coding (83‑88/1224 nt)	deoB →	phosphopentomutase

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	2514274 =	NA (NA)	34 (0.970)	30/242	0.3	100%	noncoding (2/1345 nt)	IS186	repeat region
	?	NC_000913	= 4619690	0 (0.000)					coding (88/1224 nt)	deoB	phosphopentomutase
*	?	NC_000913	= 2515616	NA (NA)	32 (0.920)	27/242	0.4	100%	noncoding (1344/1345 nt)	IS186	repeat region
	?	NC_000913	4619685 =	0 (0.000)					coding (83/1224 nt)	deoB	phosphopentomutase

GTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCC  >  NC_000913/2514128‑2514384                                                                                                                                                   |                                                                                                               gTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTcc                                                                                                               >  3:137176/1‑148 (MQ=255)                  agagCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTc                                                                                                                <  1:145243/130‑1 (MQ=255)                  agagCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTc                                                                                                                >  2:145243/1‑130 (MQ=255)                            aaCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCAt                                                                                                            >  1:91043/1‑124 (MQ=255)                            aaCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCAt                                                                                                            <  2:91043/124‑1 (MQ=255)                             aCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTaa                                                                                              >  7:27072/1‑137 (MQ=255)                             aCGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTaa                                                                                              <  8:27072/137‑1 (MQ=255)                                             cGGTCTGAAGCTGGTTTACGTCTTTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGTTGCTGTAAGACCATACATAAAAAAGCCGTGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                                                                    <  7:48921/147‑1 (MQ=255)                                                           tttACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                      <  4:157477/147‑1 (MQ=255)                                                           tttACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt                                                      >  3:157477/1‑147 (MQ=255)                                                             tACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACt                                                                                                 >  7:147240/1‑102 (MQ=255)                                                             tACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACt                                                                                                 <  8:147240/102‑1 (MQ=255)                                                                      aCCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                         >  5:105695/1‑149 (MQ=255)                                                                       cccTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                         >  2:150166/1‑148 (MQ=255)                                                                       cccTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCAc                                         <  4:137176/148‑1 (MQ=255)                                                                                        tCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTa                                                                                               >  3:90100/1‑77 (MQ=255)                                                                                        tCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTa                                                                                               <  4:90100/77‑1 (MQ=255)                                                                                             gTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                  <  6:105695/149‑1 (MQ=255)                                                                                                     cgtcgatGGTGCTGTAAGACCATACATAAAAAAGCCGGGGACAATACCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCGGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGGCAGCAATTCTGGCCCATATTg          <  8:85572/143‑1 (MQ=17)                                                                                                           tGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTaa      >  7:144067/1‑147 (MQ=255)                                                                                                           tGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTaa      <  8:144067/147‑1 (MQ=255)                                                                                                                            acataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACt                                  >  8:147927/1‑102 (MQ=255)                                                                                                                            acataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACt                                  <  7:147927/102‑1 (MQ=255)                                                                                                                              ataAAAAAGCCGTGGATAATTCCGATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACTTGCCGATGATTTACTCTGACGATAACTGGTCATCAATTCTGGc                  <  3:128806/116‑1 (MQ=14)                                                                                                                              ataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGt                               <  8:87376/103‑1 (MQ=255)                                                                                                                              ataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGt                               >  7:87376/1‑103 (MQ=255)                                                                                                                              ataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGc                  >  4:128806/1‑116 (MQ=255)                                                                                                                                aaaaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa                                                       <  7:74820/77‑1 (MQ=255)                                                                                                                                aaaaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa                                                       >  8:74820/1‑77 (MQ=255)                                                                                                                                  aaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCAATGAAt                                                  >  7:31348/1‑80 (MQ=255)                                                                                                                                  aaaaGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCAATGAAt                                                  <  8:31348/80‑1 (MQ=255)                                                                                                                                         ggggATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg          >  6:115600/1‑113 (MQ=35)                                                                                                                                         ggggATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTg          <  5:115600/113‑1 (MQ=35)                                                                                                                                               aaTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGccc  >  2:178693/1‑115 (MQ=32)                                                                                                                                               aaTTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGccc  <  1:178693/115‑1 (MQ=32)                                                                                                                                                   |                                                                                                               GTTAAAGGCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCC  >  NC_000913/2514128‑2514384

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