Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,351,197	IS186 (–) +6 bp :: Δ1 bp	intergenic (‑158/+205)	fabI ← / ← ycjD	enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	608008 =	NA (NA)	36 (0.730)	31/250	0.8	100%	noncoding (2/1345 nt)	IS186	repeat region
	?	NC_000913	1351197 =	0 (0.000)					intergenic (‑158/+210)	fabI/ycjD	enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein
*	?	NC_000913	= 609351	NA (NA)	38 (0.780)	31/248	0.8	97.3%	noncoding (1345/1345 nt)	IS186	repeat region
	?	NC_000913	= 1351202	1 (0.020)					intergenic (‑163/+205)	fabI/ycjD	enoyl‑[acyl‑carrier‑protein] reductase, NADH‑dependent/DUF559 family endonuclease‑related protein

GCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCT  >  NC_000913/609213‑609495                                                                                                                                           |                                                                                                                                                 gCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATcc                                                                                                                                           <  1:213087/146‑1 (MQ=35) gCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATcc                                                                                                                                           >  2:213087/1‑146 (MQ=35)         gTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCg                                                                                                                                >  2:170297/1‑149 (MQ=255)                   gTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtcaggt                                                                                                                      >  1:196075/1‑149 (MQ=255)                        aCAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCg                                                                                                                                <  4:13109/134‑1 (MQ=255)                        aCAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCg                                                                                                                                >  3:13109/1‑134 (MQ=255)                        aCAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGa                                                                                                                 <  2:196075/149‑1 (MQ=255)                                ccTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtc                                                                                                                          <  2:52691/132‑1 (MQ=255)                                ccTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtc                                                                                                                          >  1:52691/1‑132 (MQ=255)                                  tGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtc                                                                                                                          <  8:13375/130‑1 (MQ=255)                                  tGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTgtc                                                                                                                          >  7:13375/1‑130 (MQ=255)                                     cAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCt                                                                                                    >  8:209174/1‑149 (MQ=255)                                         cTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCtt                                                                                                <  7:209174/149‑1 (MQ=255)                                          tACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCTCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCtt                                                                                                >  8:198712/1‑148 (MQ=255)                                               aaCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACt                                                                                          >  3:208474/1‑149 (MQ=255)                                                aCGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACtt                                                                                         >  4:119487/1‑149 (MQ=255)                                                         aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTg                                                                                <  4:156298/149‑1 (MQ=255)                                                         aagaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTg                                                                                <  1:170297/149‑1 (MQ=255)                                                           gaaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGt                                                                              >  4:33449/1‑149 (MQ=255)                                                            aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTg                                                                                <  2:60352/146‑1 (MQ=255)                                                            aaCGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTg                                                                                >  1:60352/1‑146 (MQ=255)                                                                      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCt                                                                   >  7:15539/1‑149 (MQ=255)                                                                         gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCg                                                                >  2:212603/1‑149 (MQ=255)                                                                              aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCg                                                            >  6:173172/1‑148 (MQ=255)                                                                              aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCg                                                            >  1:158477/1‑148 (MQ=255)                                                                               aCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGaa                                                          <  2:158477/149‑1 (MQ=255)                                                                                      aCCGCGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTAcc                                                   <  5:173172/149‑1 (MQ=255)                                                                                                 gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTc                                         >  8:209924/1‑148 (MQ=255)                                                                                                  cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCaa                                       <  3:33449/149‑1 (MQ=255)                                                                                                   aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTaa                                                                                                           >  8:135195/1‑80 (MQ=255)                                                                                                   aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTaa                                                                                                           <  7:135195/80‑1 (MQ=255)                                                                                                      cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTa                                    >  1:62416/1‑148 (MQ=255)                                                                                                       ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAAt                                  >  6:150900/1‑149 (MQ=255)                                                                                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGTAGTAATACCAGCATTTCCTGTTGTCAGGTCATGATGGTAAGTCACCTGCTTTACACTTTCGCCCGTTGTCAGTTATGGCTGCGGGCGAATAGTAACAGATGTTGTCAATTAATCgtgt                             <  4:208474/149‑1 (MQ=255)                                                                                                                    aCCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCa                                                  <  3:195982/120‑1 (MQ=255)                                                                                                                    aCCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCa                                                  >  4:195982/1‑120 (MQ=255)                                                                                                                     cccTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGc                    <  8:15539/149‑1 (MQ=255)                                                                                                                     cccTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGc                    <  1:212603/149‑1 (MQ=255)                                                                                                                      ccTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCg                   >  2:98508/1‑149 (MQ=255)                                                                                                                        ttAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCgt                               >  3:37268/1‑135 (MQ=255)                                                                                                                        ttAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCgt                               <  4:37268/135‑1 (MQ=255)                                                                                                                         tAAGTTAGCGCTTTTGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTa                >  8:204182/1‑149 (MQ=255)                                                                                                                          aaGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGtt                 <  5:150900/147‑1 (MQ=255)                                                                                                                             ttAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGAc            >  2:214627/1‑149 (MQ=255)                                                                                                                                gcgcTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTAt         >  5:10225/1‑149 (MQ=255)                                                                                                                                  gcTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTAtct       >  3:94824/1‑149 (MQ=255)                                                                                                                                    ttATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTaa                                   >  2:70050/1‑119 (MQ=255)                                                                                                                                    ttATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTaa                                   <  1:70050/119‑1 (MQ=255)                                                                                                                                       tGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTAtcttttct  <  7:209924/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                 GCTTACAGGTGGCCATTCGTGGGACAGTATCCCTGACAGCCTACAAAACGCAATTGAAGAACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGAGTAATCCCCGCATATCCGGTTGTCAGGTCAGGATGGTAAGGCACCTGCTTTACACTTTCGCCCGTGGTCAGTGATGGCTGCGGGCGAATCGTACCAGATGTTGTCAATTAATCGTGTTGGCACAGCGTTATGACTATCTTTTCT  >  NC_000913/609213‑609495

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