New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_hsa_tpiA = 2267299 (0.230)5 (0.130) 5/272 4.6 36.4% noncoding (971/2904 nt) rrlH 23S ribosomal RNA of rrnH operon
?NC_000913_3_hsa_tpiA 2728240 = NA (NA)noncoding (945/2904 nt) rrlG 23S ribosomal RNA of rrnG operon
Rejected: Coverage evenness skew score above cutoff.

TTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAAGAGGGAAACAACCCAGACCGCC  >  NC_000913_3_hsa_tpiA/226584‑226753
                                                                                                                                                 |                        
tttAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGa                         >  1:179044/1‑147 (MQ=32)
tttAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGa                         <  2:179044/147‑1 (MQ=32)
    ggTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGACAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAAGAGGGa                  <  1:463069/150‑1 (MQ=17)
           ccTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAAGAGGGAAACAAcc           <  1:582363/150‑1 (MQ=32)
             tCGTGAACTCATCGCCGAGGGTAGAGTACTGTTTCGGCAAGGGGGTCATCGCGACTTACTAACCCGATGGAAACAGCGACTACTGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGACGTGAAGAGGCAAACAACGCa         >  1:572966/1‑150 (MQ=11)
               gTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGCCACACGGCGGGTCATAACCTCCGTCGTGAACAGCGGAACAACCCaac       >  2:193492/1‑148 (MQ=11)
               gTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCACCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAAGAGGGAAACAACCCAGa       >  1:471510/1‑150 (MQ=17)
                   aCTCATTTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAAGAGGGAAACAACCCAGACCGc   <  1:282947/150‑1 (MQ=14)
                    cTCATCTCCGGGGGGAGAGCACTGTGTCGGCAAGGGGGTCATCACGACTTACCAACCCGATGCAAACTGAGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGCGCTAACGTCCGTCGTGAAGGCGGAAACAACCCAGACCGcc  >  1:533628/1‑150 (MQ=11)
                                                                                                                                                 |                        
TTTAGGTAGCGCCTCGTGAACTCATCTCCGGGGGTAGAGCACTGTTTCGGCAAGGGGGTCATCCCGACTTACCAACCCGATGCAAACTGCGAATACCGGAGAATGTTATCACGGGAGACACACGGCGGGTGCTAACGTCCGTCGTGAAGAGGGAAACAACCCAGACCGCC  >  NC_000913_3_hsa_tpiA/226584‑226753

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.