New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3783137 =24 (0.800)4 (0.140) 4/208 2.9 18.0% coding (525/1020 nt) gpsA glycerol‑3‑phosphate dehydrogenase (NAD+)
?NC_000913 3783176 = 14 (0.500)coding (486/1020 nt) gpsA glycerol‑3‑phosphate dehydrogenase (NAD+)
Rejected: Frequency below/above cutoff threshold.

TCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGCCTAAG  >  NC_000913/3783007‑3783281
                                                                                                                                  |                                                                                                                                                
tCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACtttt                                                                                                                                                <  3:59453/133‑1 (MQ=255)
tCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACtttt                                                                                                                                                >  4:59453/1‑133 (MQ=255)
   gcgcCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACtttt                                                                                                                                                <  1:81689/130‑1 (MQ=255)
   gcgcCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACtttt                                                                                                                                                >  2:81689/1‑130 (MQ=255)
     gcCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGctgctg                                                                                                                           <  3:43307/149‑1 (MQ=255)
                        cAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCg                                                                                                                 >  6:162818/1‑140 (MQ=255)
                        cAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCg                                                                                                                 <  5:162818/140‑1 (MQ=255)
                          aaaCCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCt                                                                                                       >  5:73943/1‑148 (MQ=255)
                                 tcaaattccgctggcggttatctctggcccaacgtttgcgaaagaactggcggcaggtttaccgacagctatttcgctggcctcgaccgatcagaccTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCaaa                                                                                                            <  3:81910‑M1/39‑1 (MQ=255)
                                 tcaaattccgctggcggttatctctggcccaacgtttgcgaaagaactggcggcaggtttaccgacagctatttcgctggcctcgaccgatcagaccTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCaaa                                                                                                            >  4:81910‑M1/98‑136 (MQ=255)
                                            aTCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTgg                                                                                                                          >  8:85701/1‑111 (MQ=255)
                                                     ccAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGatcat                                                                                                                   >  3:54147/1‑109 (MQ=255)
                                                     ccAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGatcat                                                                                                                   <  4:54147/109‑1 (MQ=255)
                                                      cAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTc                                                                          >  4:103677/1‑149 (MQ=255)
                                                            cccaacgtttgcgaaagaactggcggcaggtttaccgacagctatttcgctggcctcgaccgatcagaccTTTg                                                                                                                                               >  4:64056‑M1/71‑74 (MQ=255)
                                                                       ttCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCa                                                          <  3:103677/148‑1 (MQ=255)
                                                                                               cgacagctatttcgctggcctcgaccgatcagaccTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTc                                                                                                        >  5:4689‑M1/36‑78 (MQ=255)
                                                                                               cgacagctatttcgctggcctcgaccgatcagaccTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTc                                                                                                        <  6:4689‑M1/43‑1 (MQ=255)
                                                                                                                 cTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGGCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGaa                <  6:73943/148‑1 (MQ=255)
                                                                                                                  tGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAAtt              >  1:45562/1‑149 (MQ=255)
                                                                                                                          cgGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAg                     <  1:135496/134‑1 (MQ=255)
                                                                                                                          cgGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAg                     >  2:135496/1‑134 (MQ=255)
                                                                                                                           ggAAACTTATGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCATTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGcc      >  5:144030/1‑148 (MQ=255)
                                                                                                                             aaaCTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGCCTaa   <  2:45562/149‑1 (MQ=255)
                                                                                                                              aaCGTTTGCCGCAGTGCAGCAGCTGCGGGAGATCATCGGCAAAGGTCTGATCGGGCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGCCTAAg  <  6:144030/149‑1 (MQ=255)
                                                                                                                                 ttttGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAg                                                                                            <  7:39127/56‑1 (MQ=255)
                                                                                                                                 ttttGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAg                                                                                            >  8:39127/1‑56 (MQ=255)
                                                                                                                                 ttttGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCttt                                                   <  5:147714/97‑1 (MQ=255)
                                                                                                                                 ttttGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCttt                                                   >  6:147714/1‑97 (MQ=255)
                                                                                                                                  |                                                                                                                                                
TCAGCGCCGTACGCGCATTCGCACCAAAACCGATACCGTCGGACATCCCCGCACCAATGGCAATAACGTTTTTCACCGCGCCGCCAAGCTGCACGCCAATGAAATCCGGATTGCTGTAAACGCGGAAACTTTTGCCGCAGTGCAGCAGCTGCTGGAGATCATCGGCAAAGGTCTGATCGGTCGAGGCCAGCGAAATAGCTGTCGGTAAACCTGCCGCCAGTTCTTTCGCAAACGTTGGGCCAGAGATAACCGCCAGCGGAATTTGATCGCCTAAG  >  NC_000913/3783007‑3783281

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.