New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 871513120 (1.770)4 (0.060) 4/250 NT 3.5% coding (547/921 nt) gsiC glutathione ABC transporter permease
?NC_000913 = 871541 107 (1.650)coding (575/921 nt) gsiC glutathione ABC transporter permease
Rejected: Frequency below/above cutoff threshold.

TGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/871364‑871513
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcggGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGAGCGTCAGGGAGGGTAAAATGTAGTGCTGCCAACTGTCTGCTCCCACGGTAGGCAGCCAGCCCAGTTCAACGGAGAATACCTGAATTAAAAGCATCCCCAGAGCAAATGCCG  <  NC_000913/871541‑871394
                                                                                                                                                                                                                                                                                                          
TGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCG                                                                                                                                                       >  8:74856/1‑149
   CCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGG                                                                                                                                                      >  5:237300/1‑147
                                                                                           TACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGA                                                                                                                 <  3:333360/96‑1
                                                                                           TACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGA                                                                                                                 >  4:333360/1‑96
                                                                                               GTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGAG                                                                                                                <  5:309022/93‑1
                                                                                               GTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGGGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGAG                                                                                                                >  6:309022/1‑93
                                                                                                                                                GCGCGGGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGAGCGTCAGGGAGGGTAAAATGTAGTGCTGCCAACTGTCTGCTCCCACGGTAGGCAGCCAGCCCAGTTCAACGGAGAATACCTGAATTAAAAGCATCCCCAGAGCAAA       <  7:11922/149‑1
                                                                                                                                                 CGCGGGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGAGCGTCAGGGAGGGTAAAATGTAGTGCTGCCAACTGTCTGCTCCCACGGTAGGCAGCCAGCCCAGTTCAACGGAGAATACCTGAATTAAAAGCATCCCCAGAGCAAAT      <  1:98441/149‑1
                                                                                                                                                    GGGTAAAGCGCGCCATCACTGCGGCCACTGCCGCGCCGAGCGTCAGGGAGGGTAAAATGTAGTGCTGCCAACTGTCTGCTCCCACGGTAGGCAGCCAGCCCAGTTCAACGGAGAATACCTGAATTAAAAGCATCCCCAGAGCAAATGCC   <  4:231474/149‑1
                                                                                                                                                     GGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGAGCGTCAGGGAGGGTAAAATGTAGTGCTGCCAACTGTCTGCTCCCACGGTAGGCAGCCAGCCCAGTTCAACGGAGAATACCTGAATTAAAAGCATCCCCAGAGCAAATGCCG  <  2:33663/149‑1
                                                                                                                                                     GGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGAGCGTCAGGGAGGGTAAAATGTAGTGCTGCCAACTGTCTGCTCCCACGGTAGGCAGCCAGCCCAGTTCAACGGAGAATACCTGAATTAAAAGCATCCCCAGAGCAAATGCCG  <  5:372772/149‑1
                                                                                                                                                                                                                                                                                                          
TGACCATTGCGGTGTCGGGGATCTCGTTTCCGGCATTTGCTCTGGGGATGCTTTTAATTCAGGTATTCTCCGTTGAACTGGGCTGGCTGCCTACCGTGGGAGCAGACAGTTGGCAGCACTACATTTTACCCTCCCTGACGCTCGGCGCGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/871364‑871513
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gcgcggGTAAAGCGCGCCATCACGGCGGCCACTGCCGCGCCGAGCGTCAGGGAGGGTAAAATGTAGTGCTGCCAACTGTCTGCTCCCACGGTAGGCAGCCAGCCCAGTTCAACGGAGAATACCTGAATTAAAAGCATCCCCAGAGCAAATGCCG  <  NC_000913/871541‑871394

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