New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 63007920 (0.980)5 (0.260) 5/284 1.5 19.6% coding (186/2106 nt) cstA carbon starvation protein involved in peptide utilization; APC peptide transporter family protein
?NC_000913 = 630098 22 (1.140)coding (205/2106 nt) cstA carbon starvation protein involved in peptide utilization; APC peptide transporter family protein
Rejected: Frequency below cutoff threshold.

CTGTAATGGGAGCATTTGCTCTGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATG  >  NC_000913/629934‑630226
                                                                                                                                                 |                                                                                                                                                   
ctGTAATGGGAGCATTTGCTCTGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTgg                                                                                                                                                >  2:230440/1‑150 (MQ=255)
                  ctctGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGAcgc                                                                                                                               >  1:211039‑M2/1‑147 (MQ=255)
                  ctctGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGAcgc                                                                                                                               <  2:11782‑M2/150‑4 (MQ=255)
                  ctctGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGAcgc                                                                                                                               >  1:82555‑M2/1‑147 (MQ=255)
                      gggATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGAcgccagca                                                                                                                          <  2:265363‑M2/151‑9 (MQ=255)
                      gggATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCg                                                                                                                          >  1:187860/1‑151 (MQ=255)
                      gggATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCg                                                                                                                          >  2:88671/1‑151 (MQ=255)
                        gATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGttctacgcgtcgggtcaaccgccagccca                                                                                                                        >  2:87798‑M1/1‑122 (MQ=255)
                        gATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCgg                                                                                                                         <  1:107221/150‑1 (MQ=255)
                                    tAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAacga                                                                                                            <  1:8210/151‑1 (MQ=255)
                                          gTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACggg                                                                                                        >  1:17131/1‑149 (MQ=255)
                                               gAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGAc                                                                                                  <  1:26751/150‑1 (MQ=255)
                                                aaCAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACt                                                                                                 <  1:36899/150‑1 (MQ=255)
                                                             gcTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGa                                                                                    >  1:115284/1‑150 (MQ=255)
                                                                gtgGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAg                                                                                >  2:96936/1‑151 (MQ=255)
                                                                     ttGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGTCATACGCGTCGGGTCAACcgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatag                                                                           <  1:97454‑M2/151‑56 (MQ=255)
                                                                       gtggtgGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGTCATACGCGTCGGGTCAACcgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagat                                                                         <  1:265365‑M2/151‑58 (MQ=255)
                                                                       gtggtgGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTg                                                                          >  2:92236/1‑150 (MQ=255)
                                                                              cGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGgttgacccgacgcgtatgacgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagac                                                                   <  2:284514‑M1/150‑83 (MQ=255)
                                                                                   gTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGTCATACGCGTCGGGTCAACcgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagaccgacg                                                              <  2:284514‑M2/150‑69 (MQ=255)
                                                                                   gTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCAc                                                              >  1:83362/1‑150 (MQ=255)
                                                                                   gTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCAc                                                              >  1:55174/1‑150 (MQ=255)
                                                                                   gTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCAc                                                              >  2:222407/1‑150 (MQ=255)
                                                                                     cTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCAt                                                           >  1:215909/1‑151 (MQ=255)
                                                                                       gTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCAttt                                                         >  2:305300/1‑151 (MQ=255)
                                                                                        tatctatcTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAATAAAGTGCTGTTCGGTCACCAttt                                                         <  2:159548/150‑1 (MQ=255)
                                                                                         atctatctGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGgttgacccgacgcgtatgacgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagaccgacgccaccac                                                       <  1:265365‑M1/151‑95 (MQ=255)
                                                                                           ctatctGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGgttgacccgacgcgtatgacgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagaccgacgccaccacaa                                                     <  1:97454‑M1/151‑97 (MQ=255)
                                                                                                            gTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAgg                                    >  2:183968/1‑151 (MQ=255)
                                                                                                               tttATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTc                                  >  1:233121/1‑150 (MQ=255)
                                                                                                                         gTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTgggg                       >  2:268908/1‑151 (MQ=255)
                                                                                                                                        tgggCTGGCGGTTGACCCGACGCGTAGAAcgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagaccgacgccaccacaatccacagcgcgttgatctgttccccacgatttaaagcaatgtatc        >  2:87798‑M2/1‑29 (MQ=255)
                                                                                                                                          tgCTGGCGtcatacgcgtcgggtcaaccgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagaccgacgccaccacaatccacagcgcgttgatctgttccccacgatttaaagcaatgtatccc      <  2:265363‑M1/151‑144 (MQ=255)
                                                                                                                                               gCGtcatacgcgtcgggtcaaccgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagaccgacgccaccacaatccacagcgcgttgatctgttccccacgatttaaagcaatgtatcccagag  >  1:211039‑M1/1‑3 (MQ=255)
                                                                                                                                               gCGtcatacgcgtcgggtcaaccgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagaccgacgccaccacaatccacagcgcgttgatctgttccccacgatttaaagcaatgtatcccagag  <  2:11782‑M1/150‑148 (MQ=255)
                                                                                                                                               gCGtcatacgcgtcgggtcaaccgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagaccgacgccaccacaatccacagcgcgttgatctgttccccacgatttaaagcaatgtatcccagag  >  1:82555‑M1/1‑3 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CTGTAATGGGAGCATTTGCTCTGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATG  >  NC_000913/629934‑630226

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