Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
JC	NC_000913	168,071	+AGAA	coding (588/2244 nt)	fhuA →	ferrichrome outer membrane transporter/phage receptor

New junction evidence
		seq id	position	reads (cov)	reads (cov)	score	skew	freq	annotation	gene	product
*	?	NC_000913	= 168064	0 (0.000)	30 (0.950) +GAAA	20/138	0.5	100%	coding (581/2244 nt)	fhuA	ferrichrome outer membrane transporter/phage receptor
	?	NC_000913	168065 =	0 (0.000)					coding (582/2244 nt)	fhuA	ferrichrome outer membrane transporter/phage receptor

GGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167987‑168064 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGA  >  NC_000913/168065‑168142                                                                               ||||                                                                               GGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGA                                                                                  >  1:29630/1‑80  GTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAA                                                                                <  1:29604/81‑1  GTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAA                                                                                <  1:29605/81‑1  GTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGA                                                                                  >  2:29629/1‑79    TTCCGTGCTTTACGGTAAAAGCAGTCGTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAG                                                                               >  1:29636/1‑80      CCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGA                                                                              >  2:29642/1‑79      CCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGA                                                                              >  2:29633/1‑79       CGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAG                                                                           <  1:29606/81‑1        GTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAG                                                                           >  1:29635/1‑80                    AAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTA                                                               >  1:29650/1‑80                        GCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCC                                                          >  2:29639/1‑81                        GCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCC                                                          <  2:29607/81‑1                                  CGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACA                                                <  1:29608/81‑1                                      CTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCT                                            <  1:29609/81‑1                                         TTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTG                                           <  1:29610/79‑1                                            AATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTC                                        >  2:29646/1‑79                                                  GTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGG                                <  2:29611/81‑1                                                    CAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTT                              >  2:29640/1‑81                                                     AGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTT                              <  2:29612/80‑1                                                     AGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTT                              <  2:29613/80‑1                                                        AAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGA                          >  1:29648/1‑81                                                        AAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGA                          >  1:29645/1‑81                                                        AAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGA                          >  2:29637/1‑81                                                        AAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGA                          >  1:29644/1‑81                                                          GCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACT                        <  2:29614/81‑1                                                             TCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTA                     <  2:29615/81‑1                                                              CCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAG                    >  2:29647/1‑81                                                                 ACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGA                 <  1:29616/81‑1                                                                    ACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTC              <  2:29617/81‑1                                                                     CCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCG             >  2:29643/1‑81                                                                           CGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGAT       <  1:29618/81‑1                                                                               GAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGAT    >  2:29675/1‑80                                                                                AAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGA  <  2:29619/81‑1                                                                               ||||                                                                               GGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167987‑168064 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGA  >  NC_000913/168065‑168142

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