Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 1,426,494 IS2 (+) +5 bp 15.3% pseudogene (41‑45/933 nt) tmpR → Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1426494 =105 (0.800)22 (0.160) 18/176 NT 16.5% pseudogene (41/933 nt) tmpR Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
?NC_000913 = 1651548 NA (NA)noncoding (706/706 nt) IS2 repeat region
* ? NC_000913 = 1426498105 (0.790)16 (0.110) 15/178 NT 12.6% pseudogene (45/933 nt) tmpR Rac prophage; pseudogene, tail protein family;Phage or Prophage Related; putative alpha helix protein
?NC_000913 = 2070271 NA (NA)noncoding (1/1331 nt) IS2 repeat region
Rejected: Frequency below cutoff threshold.

CCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1426405‑1426498
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGA  <  NC_000913/2070271‑2070173
                                                                                                                                                                                                 
CCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACT                                                                                               <  1:1508059/100‑1
 CTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACTG                                                                                              >  1:1186826/1‑100
                  GTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACTGGCCCCCTGAATCTCCAG                                                                             >  2:2333292/1‑100
                      GTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTAATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAA                                                                         >  1:972204/1‑100
                             GTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCAATTAGACTGGCCCCCTGAATCTCCAGACAACCAATAT                                                                  >  2:3693905/1‑100
                               TGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCA                                                                >  1:443745/1‑100
                                GTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCAC                                                               <  2:1916498/100‑1
                                   AGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTA                                                            <  2:1687450/100‑1
                                              TAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATA                                                 <  2:1551133/100‑1
                                                      AGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACTGGCCCGCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAAT                                         <  1:1834426/100‑1
                                                        AGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATAC                                       <  2:3276672/100‑1
                                                           AAACCTCGCGTCTCGCAATAATTATTGATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAG                                    <  1:3313856/100‑1
                                                                 CGCGTCTCGCAATAATTATTGATAGCACTTAGACTGGCACCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGT                                   <  1:3808035/95‑1
                                                                              AATTATTGATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGG                 >  2:598621/1‑100
                                                                                     GATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAG          >  2:1745355/1‑100
                                                                                     GATAGCACTTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAG          >  2:595775/1‑100
                                                                                             TTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGA  >  2:166711/1‑100
                                                                                             TTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGA  <  2:1191859/100‑1
                                                                                             TTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGA  >  2:1739917/1‑100
                                                                                             TTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGA  <  2:2788346/100‑1
                                                                                             TTAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGA  <  2:480651/100‑1
                                                                                                                                                                                                 
CCTGAAACCCACTTCGCGGTGGGTTTTTTGTTGTCAGGAGTTTTAATAAATGGCAGAGCAAACCTCGCGTCTCGCAATAATTATTGATAGCACT‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/1426405‑1426498
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑TAGACTGGCCCCCTGAATCTCCAGACAACCAATATCACTTAAATAAGTGATAGTCTTAATACTAGTTTTTAGACTAGTCATTGGAGAACAGATGATTGA  <  NC_000913/2070271‑2070173

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