Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC JC	NC_000913	1,416,518	IS5 (–) +4 bp	coding (866‑869/2601 nt)	recE ←	Rac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	1416518 =	0 (0.000)	43 (0.680)	41/298	1.1	100%	coding (869/2601 nt)	recE	Rac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
	?	NC_000913	= 2290113	NA (NA)					noncoding (1/1195 nt)	IS5	repeat region
*	?	NC_000913	= 1416521	0 (0.000)	86 (1.360)	55/298	0.3	100%	coding (866/2601 nt)	recE	Rac prophage; exonuclease VIII, 5' to 3' specific dsDNA exonuclease
	?	NC_000913	2101749 =	NA (NA)					noncoding (1195/1195 nt)	IS5	repeat region

TGCTAATGCTGGCGGCGGAAACCTCACCGATCGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1416661‑1416518 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAA  <  NC_000913/2290113‑2289984                                                                                                                                                                                                                                                                                    TGCTAATGCTGGCGGCGGAAACCTCACCGATCGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGG                                                                                                                              >  2:75634/1‑150  GCTAATGCTGGCGGCGGAAACCTCACCGATCGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTG                                                                                                                            <  2:978060/151‑1           GGCGGCGGAAACCTCACCGATCGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAG                                                                                                                    <  2:649774/150‑1           GGCGGCGGAAACCTCACCGATCGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAG                                                                                                                    <  1:682076/150‑1               GCGGAAACCTCACCGATCGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCT                                                                                                               <  1:895774/151‑1                CGGAAACCTCACCGATCGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCTACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTG                                                                                                              <  1:701887/151‑1                  GAAACCTCACCGATCGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGAT                                                                                                            <  1:982637/151‑1                                CGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTT                                                                                              >  2:273149/1‑151                                CGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTT                                                                                               <  1:857659/150‑1                                   GGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTC                                                                                           >  1:408933/1‑151                                    GCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCA                                                                                          <  1:1007726/151‑1                                     CGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCAT                                                                                         <  1:333463/151‑1                                       AAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCT                                                                                       <  1:87681/151‑1                                          GTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGG                                                                                     >  1:670618/1‑150                                            TTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAG                                                                                   <  1:966850/150‑1                                                   ACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAA                                                                           <  1:20018/151‑1                                                    CGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAA                                                                           <  1:718700/150‑1                                                      ATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAAC                                                                          >  1:329205/1‑149                                                         TGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGT                                                                     >  2:749035/1‑151                                                          GACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGG                                                                      <  1:765177/149‑1                                                          GACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAG                                                                        <  1:126540/147‑1                                                                 CTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATG                                                             >  1:324084/1‑151                                                                    GCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGT                                                          >  1:740353/1‑151                                                                        GCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATC                                                      >  2:309188/1‑151                                                                            CGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACT                                                  >  1:952454/1‑151                                                                               AGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTAC                                               <  1:700151/151‑1                                                                                 CCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCT                                             <  2:779739/151‑1                                                                                      GGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCC                                        <  1:942994/151‑1                                                                                         GTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGAC                                     >  1:156346/1‑151                                                                                             TGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTG                                 >  2:894933/1‑151                                                                                              GGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGA                                <  1:1005938/151‑1                                                                                              GGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTG                                 <  2:12307/150‑1                                                                                                CCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAAT                              <  2:789313/151‑1                                                                                                  CGTTCAATGGATCTGGACATCTAT‑ACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCCTGTTTGTCATCTGGAGCCATAGAACCGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTC                            >  2:780938/1‑150                                                                                                   GTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCA                           <  2:795878/151‑1                                                                                                     TCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGC                         >  2:55468/1‑151                                                                                                      CAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCA                        <  2:809968/151‑1                                                                                                        ATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGT                      >  2:1027921/1‑151                                                                                                                GGACATCTATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGC                <  2:584356/149‑1                                                                                                                        ATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCA      <  2:840624/151‑1                                                                                                                        ATAACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGA         <  2:408826/148‑1                                                                                                                          AACCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGA    <  2:230100/151‑1                                                                                                                            CCTTCATCCGGCACACGCTAAGGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAA  <  1:168238/151‑1                                                                                                                                                                                                                                                                                    TGCTAATGCTGGCGGCGGAAACCTCACCGATCGCGGCGAAGGTTTCGTACACGATCTGACGTCACTGGCGCGCGACGTAGCCACTGGCGTACTGGCCCGTTCAATGGATCTGGACATCTATAACCTTCATCCGGCACACGCTAA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1416661‑1416518 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GGAAGGTGCGAACAAGTCCCTGATATGAGATCATGTTTGTCATCTGGAGCCATAGAACAGGGTTCATCATGAGTCATCAACTTACCTTCGCCGACAGTGAATTCAGCAGTAAGCGCCGTCAGACCAGAAA  <  NC_000913/2290113‑2289984

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