New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 435495039 (0.860)4 (0.090) 4/242 NT 8.9% coding (106/681 nt) modB molybdate ABC transporter permease ModB
?NC_002947 = 4354978 44 (1.000)coding (134/681 nt) modB molybdate ABC transporter permease ModB

ATCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/4354824‑4354950
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcCCGCGCAGCCAAGAGCGCGTGCGCGCGAGCCACCAG  <  NC_002947/4354978‑4354943
                                                                                                                                                                   
ATCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCCCGCGCAGCCAAGAGCGCGTGCG               >  1:406226/1‑150
        GTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTG                                        <  1:907136/117‑1
        GTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTG                                        >  2:907136/1‑117
             GAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCCCGCGCAGCCAAGAGCGCGTGCGCGCGAGCCACCAG  <  2:406226/150‑1
                         CACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCCCGCGCAGCCAAGAGCGCGTGCGCGCGAGCCACCA   <  1:375275/137‑1
                         CACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGCCCGCGCAGCCAAGAGCGCGTGCGCGCGAGCCACCA   >  2:375275/1‑137
                                                                                                                                                                   
ATCCTATGGTTATGAACGCTGATGCCACTCGATGCCAGTGACCTGGGCGCCATCTGGCTGACCGTAAAACTGGCCAGCCTGACCACCCTGATCCTGCTGATCATCGGCACACCCGTCGCCTGGTGGC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_002947/4354824‑4354950
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcCCGCGCAGCCAAGAGCGCGTGCGCGCGAGCCACCAG  <  NC_002947/4354978‑4354943

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