New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_pae_tpiA = 179295011 (0.880)3 (0.250) 3/246 1.3 20.4% coding (609/750 nt) btuD vitamin B12 ABC transporter ATPase
?NC_000913_3_pae_tpiA = 1792973 13 (1.100)coding (586/750 nt) btuD vitamin B12 ABC transporter ATPase

ATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_tpiA/1792797‑1792950
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caaccacACATTGCGTCATGCGCATCGGGCGTGGTTGCTAA  <  NC_000913_3_pae_tpiA/1792973‑1792940
                                                                                                                                                                                            
ATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAA                                        >  1:109799/1‑150
                                      GACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACACATTGCGTCATGCGCATCGGGCGTGGTTGCTAA  <  4:80639/150‑1
                                                                  CATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACACATTGCGTCATGCGCATCGGGCGTGGTTGCT    >  1:109084/1‑120
                                                                  CATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCACACATTGCGTCATGCGCATCGGGCGTGGTTGCT    <  2:109084/120‑1
                                                                                                                                                                                            
ATACTAATGAAATCAGATGGTCGAAATCAGCATTCTGTGACCTTCGATATCCAGACGGCGAAAATTCATCCCATAGGCCTGCGCCAGATTTGGCGGCGTGAGCACCTCTTCCCTGCGTCCACTGGCCAGCATTTTTCCACCTTTTAGCAACCAC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913_3_pae_tpiA/1792797‑1792950
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑caaccacACATTGCGTCATGCGCATCGGGCGTGGTTGCTAA  <  NC_000913_3_pae_tpiA/1792973‑1792940

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