Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,491,904	IS186 (+) +3 bp	intergenic (+54/+16)	ydcA → / ← hokB	putative periplasmic protein/toxic polypeptide, small

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 609351	NA (NA)	40 (0.860)	35/274	0.3	100%	noncoding (1345/1345 nt)	IS186	repeat region
	?	NC_000913	1491904 =	0 (0.000)					intergenic (+54/+18)	ydcA/hokB	putative periplasmic protein/toxic polypeptide, small
*	?	NC_000913	= 1491906	0 (0.000)	56 (1.200)	53/274	0.0	100%	intergenic (+56/+16)	ydcA/hokB	putative periplasmic protein/toxic polypeptide, small
	?	NC_000913	2514273 =	NA (NA)					noncoding (1/1345 nt)	IS186	repeat region

GCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCG  >  NC_000913/2514135‑2514402                                                                                                                                           |                                                                                                                                  gCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGCGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCa                                                                                                                               >  1:250975/1‑143 (MQ=255)   cTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTCAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATa                                                                                                                             >  2:769528/1‑143 (MQ=255)        tgaaGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCCATCGCAGCACTGGTGCTGTAAGACCATACATGAAAAAGCGGGGGATAATTCCCATGGgcgc                                                                                                                        >  1:323052/1‑143 (MQ=255)         gaaGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTTCCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCt                                                                                                                       <  1:27892/143‑1 (MQ=255)           cgagCAAGGTAACGTGGTTTCTCGCTTCGGTCGGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTaa                                                                                                                     <  2:455106/142‑1 (MQ=255)             agCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTCGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTACCATAAGCGCTAACt                                                                                                                   >  1:176566/1‑143 (MQ=255)                aaGTTAATGTGGTTTCTTGCTTCGGTATGATGCTGGTTTACGGCTTTAACCTTGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGTTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTaa                                                                                                                <  2:747561/143‑1 (MQ=255)                   gTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGCGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCAGACATAAAGAAGCCGGGGATAGTTCCCGTGAGCGCTGACTTAAgga                                                                                                             >  1:67684/1‑142 (MQ=255)                     aaCGTGGTTTCTCGTTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGtt                                                                                                           <  2:794909/143‑1 (MQ=255)                     aaCGTGGTTTCTCGCTTCGGTCGGACGCTGGTTTACGGCTCTGCCCCTGTGAGGGTGCTGTCTGCGTCAATCGCAGCACGGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGtt                                                                                                           <  1:426984/143‑1 (MQ=255)                     aaCGTGGTTTCTCGCTTCGGTATGAAGCTGGTTTCCGGCTCTACCCTGGCCAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGCAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGtt                                                                                                           <  2:337072/143‑1 (MQ=255)                       cggggTTTCTCGCTTCGGTCTGACGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTtgt                                                                                                         <  2:173484/140‑1 (MQ=255)                        gtggtgTCTCGCGTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTtgtg                                                                                                        <  1:166468/143‑1 (MQ=255)                            tttCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAACAAGCCGGGACTAATCCCCATAAGCGCTAACTTAAGGGTTGGGGTac                                                                                                    >  2:181818/1‑142 (MQ=255)                              tctcGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCCGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGTTAACTTAAGGGTTGTGGTATTa                                                                                                  <  1:64210/143‑1 (MQ=255)                              tctcGCTCCGTTCTGACGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                                                                                  <  2:703517/143‑1 (MQ=255)                                tcGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGGATTACg                                                                                                >  2:283910/1‑143 (MQ=255)                                tcGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCCGGTGAGTGTTCGCTCTGCGTCAATCGGAGCAATGGTGATGTAAGACCATACATTAAAATGCCGGGGATAATTCCCATAAGCGCTATCTTAAGGTGTGTgggttttct                                                                                                >  1:36228/1‑136 (MQ=255)                                  gCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCACTCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCTCATAAACGCTAACTTAAGGGGTGTAGTAttcacca                                                                                              >  2:265814/1‑138 (MQ=255)                                  gCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                                                                              >  1:433230/1‑143 (MQ=255)                                         tCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATACAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGGTGAGCGATTACGCCTGctagg                                                                                       >  2:723347/1‑141 (MQ=255)                                             aaGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATCAGCGCTAACTTAAGGGTTGTGGTACTACGCCTGATACGAttt                                                                                   >  2:671888/1‑143 (MQ=255)                                                  ggTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTCACGt                                                                              >  1:133390/1‑143 (MQ=255)                                                    tttGCGGCTCGAGCCTGGTTAGGGTGCTGTCTACGCCGATCGCATCACTGGTGCTGTAAGCCCATACATAAAAAAGCCTGGGATAATTCCCATAAACGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGc                                                                            <  2:790737/143‑1 (MQ=21)                                                           ctgtgcccGGGTGAGTGTGCCGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                     <  2:130797/138‑1 (MQ=255)                                                           ctctACCCGGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa                                                                     <  2:837465/143‑1 (MQ=255)                                                            tctACGCTGGTGGGTGTGCTGTCTGCGTGAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt                                                                    <  2:398001/143‑1 (MQ=255)                                                                         gtgtGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATa                                                       >  1:806767/1‑143 (MQ=255)                                                                            gcttGTCGGGGTCAATCGCAGCACTGGTGCCGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACt                                                    <  1:616995/140‑1 (MQ=255)                                                                                gTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTc                                                <  1:657472/143‑1 (MQ=255)                                                                                 tCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTTGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCa                                               >  2:874509/1‑143 (MQ=255)                                                                                    gCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTACCTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACGCTCACGATAACTGGTCAGCa                                            >  2:81104/1‑143 (MQ=255)                                                                                                 caCTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCata                               >  2:724041/1‑143 (MQ=255)                                                                                                           gTAAGACCATACATAAAGAAGCCGGGGATAATTCCTATAAGCGCTAACTTAAGGGTTTTGGTCTTCCGCCTGATATGGTTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTCGCTCGTATCGGTCAGcc                     >  1:106145/1‑143 (MQ=17)                                                                                                            tAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATACCTGGTCAGCAATTCTGGCCCATATTGGTAAGccc                    >  1:705467/1‑143 (MQ=255)                                                                                                            tAAGACCAGACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGACGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAAGTCTGGCCCAGATTCGTAAGgcc                    >  2:301802/1‑143 (MQ=17)                                                                                                                aCCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCgaag                <  2:133390/143‑1 (MQ=255)                                                                                                                   atacataTAGGAGCCGGGGATAATTCCCATAAGCGCTAGCTTAGGGGGTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAAc             <  2:576498/143‑1 (MQ=14)                                                                                                                    tacataAAAAAGCCGGGGAAAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGCTCCGCAATTCAGGCCCATATTGGTAAGCCCGAAGAACt            <  1:792750/143‑1 (MQ=255)                                                                                                                      catgAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTgg          <  2:572899/143‑1 (MQ=255)                                                                                                                       ataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCGCGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCACGACGAACTGGa         >  1:679/1‑143 (MQ=255)                                                                                                                       ataAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGa         >  2:152129/1‑143 (MQ=255)                                                                                                                        taAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCCTATTGGTAAGCCCGAAGAACTGGAt        >  2:23098/1‑143 (MQ=255)                                                                                                                              ggcgggggATAAGTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGACTTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCg  <  2:31795/139‑1 (MQ=12)                                                                                                                                           |                                                                                                                                  GCCTGAATGAAGAGCAAGGTAACGTGGTTTCTCGCTTCGGTCTGAAGCTGGTTTACGGCTCTACCCTGGTGAGTGTGCTGTCTGCGTCAATCGCAGCACTGGTGCTGTAAGACCATACATAAAAAAGCCGGGGATAATTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCACGATAACTGGTCAGCAATTCTGGCCCATATTGGTAAGCCCGAAGAACTGGATACTTCG  >  NC_000913/2514135‑2514402

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