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)	28 (0.820) +GAAA	22/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

CCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167985‑168064 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACG  >  NC_000913/168065‑168144                                                                                 ||||                                                                                 CCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTG                                                                                     <  2:31811/81‑1  CGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGA                                                                                    <  2:31812/81‑1        CCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAA                                                                               <  1:31813/80‑1         CGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGA                                                                                >  2:31856/1‑78                   GGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTT                                                                   <  2:31814/81‑1                    GTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTT                                                                  <  2:31815/81‑1                       AAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAA                                                               <  2:31816/81‑1                          GCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCC                                                            >  1:31852/1‑81                           CAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCG                                                           >  1:31843/1‑81                               CCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTAC                                                       >  1:31845/1‑81                                   GCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGAC                                                   <  1:31818/81‑1                                   GCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTG                                                     <  1:31817/79‑1                                     GGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAG                                                 >  2:31865/1‑81                                           TTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTT                                           >  2:31846/1‑81                                             GAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTC                                          <  1:31820/80‑1                                             GAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTC                                          <  2:31819/80‑1                                               ATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAG                                       <  2:31821/81‑1                                                   GGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTG                                   >  2:31868/1‑81                                                   GGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTG                                   >  2:31860/1‑81                                                          AAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGA                            >  1:31869/1‑81                                                          AAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGA                            >  2:31854/1‑81                                                          AAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGA                            >  2:31848/1‑81                                                          AAGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGA                            <  2:31822/81‑1                                                           AGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGAC                           <  2:31824/81‑1                                                           AGCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGAC                           <  1:31823/81‑1                                                            GCGTCCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACT                          <  1:31825/81‑1                                                                CCGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAG                      >  1:31855/1‑81                                                                 CGACCACCGAACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGC                     <  2:31826/81‑1                                                                          AACCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTG            <  2:31827/81‑1                                                                            CCGCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGA          >  2:31871/1‑81                                                                              GCTGAAAGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGAT         >  1:31859/1‑80                                                                                    AGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACG  >  1:31861/1‑81                                                                                 ||||                                                                                 CCGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/167985‑168064 ‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGACG  >  NC_000913/168065‑168144

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