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 = 1680640 (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

CGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGA  >  NC_000913/167986‑168142
                                                                               |                                                                             
ttccgtgctttacggtaaaagcagtcgtggcggcctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaag                                                                               >  1:29636‑M2/80‑80 (MQ=255)
  ccgtgctttacggtaaaagcagtcctggcggcctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaaga                                                                              >  2:29642‑M2/78‑79 (MQ=255)
  ccgtgctttacggtaaaagcagtcctggcggcctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaaga                                                                              >  2:29633‑M2/78‑79 (MQ=255)
   cgtgctttacggtaaaagcagtcctggcggcctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaagaaag                                                                           <  1:29606‑M2/5‑1 (MQ=255)
    gtgctttacggtaaaagcagtcctggcggcctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaagaaag                                                                           >  1:29635‑M2/76‑80 (MQ=255)
                aaaagcagtcctggcggcctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTa                                                               >  1:29650‑M2/64‑80 (MQ=255)
                    gcagtcctggcggcctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGcc                                                          >  2:29639‑M2/60‑81 (MQ=255)
                    gcagtcctggcggcctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGcc                                                          <  2:29607‑M2/22‑1 (MQ=255)
                              cggcctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACa                                                <  1:29608‑M2/32‑1 (MQ=255)
                                  ctgttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCt                                            <  1:29609‑M2/36‑1 (MQ=255)
                                     ttgaatatggtcagcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTg                                           <  1:29610‑M2/37‑1 (MQ=255)
                                        aatatggtcagcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTc                                        >  2:29646‑M2/40‑79 (MQ=255)
                                              gtcagcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTgg                                <  2:29611‑M2/48‑1 (MQ=255)
                                                cagcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGtt                              >  2:29640‑M2/32‑81 (MQ=255)
                                                 agcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGtt                              <  2:29612‑M2/50‑1 (MQ=255)
                                                 agcaagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGtt                              <  2:29613‑M2/50‑1 (MQ=255)
                                                    aagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGa                          >  1:29648‑M2/28‑81 (MQ=255)
                                                    aagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGa                          >  2:29637‑M2/28‑81 (MQ=255)
                                                    aagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGa                          >  1:29645‑M2/28‑81 (MQ=255)
                                                    aagcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGa                          >  1:29644‑M2/28‑81 (MQ=255)
                                                      gcgtccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACt                        <  2:29614‑M2/56‑1 (MQ=255)
                                                         tccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTa                     <  2:29615‑M2/59‑1 (MQ=255)
                                                          ccgaccaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAg                    >  2:29647‑M2/22‑81 (MQ=255)
                                                             accaccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGa                 <  1:29616‑M2/63‑1 (MQ=255)
                                                                accgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTc              <  2:29617‑M2/66‑1 (MQ=255)
                                                                 ccgaaccgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCg             >  2:29643‑M2/15‑81 (MQ=255)
                                                                       cgctgaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGgat       <  1:29618‑M2/73‑1 (MQ=255)
                                                                           gaaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGgatgat    >  2:29675‑M2/5‑80 (MQ=255)
                                                                            aaaGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGAtgatga  <  2:29619‑M2/78‑1 (MQ=255)
                                                                               |                                                                             
CGGTTTCCGTGCTTTACGGTAAAAGCAGTCCTGGCGGCCTGTTGAATATGGTCAGCAAGCGTCCGACCACCGAACCGCTGAAAGAAGTTCAGTTTAAAGCCGGTACTGACAGCCTGTTCCAGACTGGTTTTGACTTTAGCGATTCGTTGGATGATGA  >  NC_000913/167986‑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.