New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 2353082 =19 (0.950)3 (0.170) 3/116 1.7 14.1% coding (436/1629 nt) glpA anaerobic sn‑glycerol‑3‑phosphate dehydrogenase, large FAD/NAD(P)‑binding subunit
?NC_000913 2353145 = 20 (1.140)coding (499/1629 nt) glpA anaerobic sn‑glycerol‑3‑phosphate dehydrogenase, large FAD/NAD(P)‑binding subunit
Rejected: Frequency below/above cutoff threshold.

GATGGCACCGTTGATCCATTTCGTCTGACCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTG  >  NC_000913/2353079‑2353211
                                                                  |                                                                  
gATGGCACCGTTGATCCATTTCGTCTGACCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGt                                                                   <  2:393944/68‑1 (MQ=255)
  tGGCACCGTTGATCCATTTCGTCTGACCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGTTa                                                                 <  2:410844/68‑1 (MQ=255)
                ttggcatccagcatgtttgctgcggtcagacgaaatggatcaacggtgccGTTATCCTTACCGCTCAt                                                   >  1:42935‑M2/51‑68 (MQ=255)
                cATTTCGTCTGACCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCAt                                                   >  1:169548/1‑68 (MQ=255)
                cATTTCGTCTGACCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCAt                                                   >  2:500643/1‑68 (MQ=255)
                   atCGTCTGACCGCAGCAAACATGCAGGATGCCAAAGAACCCGGTGCCGTTAACCTTACCGATCATGaa                                                >  1:653267/2‑68 (MQ=255)
                      gTCTGACCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTc                                             <  2:727580/68‑1 (MQ=255)
                      gTCTGACCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTc                                             >  2:411794/1‑68 (MQ=255)
                      gTCTGACCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTc                                             >  2:403676/1‑68 (MQ=255)
                           aCCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACggg                                        <  2:731198/68‑1 (MQ=255)
                                 gcaAACATGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGAt                                  <  2:300988/68‑1 (MQ=255)
                                   ctgcggtcagacgaaatggatcaacggtgccGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTc                                <  2:42935‑M2/37‑1 (MQ=255)
                                        tGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGaa                           >  1:603574/1‑68 (MQ=255)
                                         gCTGGATGCAAAAGAACACGGTGCCGTTATCCTTACAGCTCATGAAGTCACGGGGCTGATTCGTGAAg                          >  2:200607/1‑68 (MQ=255)
                                                   aaaGTACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACTTGGCTGATTCGTGTAGTCGCGACGgt                <  2:653268/68‑1 (MQ=255)
                                                   aaaGAACTCGGTGCCGTTATCCTTACCGCTCTTGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGgt                <  2:349657/68‑1 (MQ=255)
                                                   aaaGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGgt                >  2:33406/1‑68 (MQ=255)
                                                   aaaGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGgt                >  2:528713/1‑68 (MQ=255)
                                                   aaaGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGgt                >  1:630875/1‑68 (MQ=255)
                                                     aGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGgtgt              >  2:735320/1‑68 (MQ=255)
                                                        caacggtgccGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCg           >  1:161715‑M2/11‑68 (MQ=255)
                                                        acacGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCg           >  2:740492/1‑68 (MQ=255)
                                                        acacGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCg           >  1:49635/1‑68 (MQ=255)
                                                                 cGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCgtg  <  1:401811‑M2/67‑1 (MQ=255)
                                                                  |                                                                  
GATGGCACCGTTGATCCATTTCGTCTGACCGCAGCAAACATGCTGGATGCCAAAGAACACGGTGCCGTTATCCTTACCGCTCATGAAGTCACGGGGCTGATTCGTGAAGGCGCGACGGTGTGCGGTGTTCGTG  >  NC_000913/2353079‑2353211

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