New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4626303111 (1.030)4 (0.040) 3/168 NT 3.4% coding (392/1383 nt) radA DNA repair protein
?NC_000913 = 4626336 120 (1.180)coding (425/1383 nt) radA DNA repair protein
Rejected: Frequency below cutoff threshold.

GCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4626200‑4626303
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgaTGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGGTGACATACAGCGTTTTCATCTGCTGGGCCAGTTTGCACAGCGTTTGCAGTAGCAGCGT  <  NC_000913/4626336‑4626239
                                                                                                                                                                                                          
GCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATTCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCG                                                                                                        >  1:893231/1‑100
GCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCG                                                                                                        >  1:1060188/1‑100
GCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCG                                                                                                        >  2:850583/1‑100
 CCATTCTGATTGGCGGTACCCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGG                                                                                                       <  2:2271447/100‑1
 CCATTCTGATTGGCGGTAACCCAGGCGCGGGGAAATCCACGCTGCTGTTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGG                                                                                                       <  2:2769876/100‑1
  CATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGG                                                                                                       >  2:2660887/1‑99
   ATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCG                                                                                                     >  2:393059/1‑100
   ATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCG                                                                                                     >  1:64366/1‑100
    TTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGA                                                                                                    <  1:2355155/100‑1
                                            GCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGATGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGG                                                            <  1:1726605/100‑1
                                             CTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGATGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGGT                                                           <  1:434378/100‑1
                                               ACTGCAAACGCTTTGCAAAATGGCCCAGCAGATGAAAACGCTGTATGTCAACGGCGATGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGGTGA                                                         <  1:526851/100‑1
                                                CTGCAAACGCTGTGCAAACTGGCACAGCAGATGAAAACGCTGTATGTCACCGGCGATGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGGTGA                                                         <  2:1898179/99‑1
                                                                                                   GGCGATGAGCGCGCATTGCCACCTGGTGCAGCGACTCGTCGCCGGTGACATACAGCGTTTTCATCTGCTGGGCCAGTTTGCACAGCGTTTGCAGTAGCAG     >  2:2423848/1‑100
                                                                                                     CGATGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGGTGACATACAGCGTTTTCATCTGCTGGGCCAGTTTGCACAGCGTTTGCAGTAGCAGCG   >  2:1070777/1‑100
                                                                                                     CGATGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGGTGACATACAGCGTTTTCATCTGCTGGGCCAGTTTGCACAGCGTTTGCAGTAGCAGCG   >  2:3020001/1‑100
                                                                                                     CGATGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGGTGACATACAGCGTTTTCATCTGCTGGGCCAATTTGCACAGCGTTTGCAGTAGCAGCG   >  1:714375/1‑100
                                                                                                      GATGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGGTGACATACAGCGTTTTCATCTGCTGGGCCAGTTTGCACAGCGTTTGCAGTAGCAGCGT  >  1:232730/1‑100
                                                                                                                                                                                                          
GCCATTCTGATTGGCGGTAACCCTGGTGCGGGGAAATCCACGCTGCTACTGCAAACGCTGTGCAAACTGGCCCAGCAGATGAAAACGCTGTATGTCACCGGCGA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4626200‑4626303
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑ggcgaTGAGCGCGCATTGCCACCTGTTGCAGCGACTCTTCGCCGGTGACATACAGCGTTTTCATCTGCTGGGCCAGTTTGCACAGCGTTTGCAGTAGCAGCGT  <  NC_000913/4626336‑4626239

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