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.

CTCTGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATG  >  NC_000913/629952‑630244
                                                                                                                                                  |                                                                                                                                                  
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)
      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)
                                                                 gTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGTCATACGCGTCGGGTCAACcgccagcacatttttggcgatatacagcccataaaaacggtaagcgatcagatagatacagaccgacg                                                                                <  2:284514‑M2/150‑69 (MQ=255)
                                                                 gTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCAc                                                                                >  2:222407/1‑150 (MQ=255)
                                                                 gTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCAc                                                                                >  1:55174/1‑150 (MQ=255)
                                                                 gTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCAc                                                                                >  1:83362/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)
                                                                                          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)
                                                                                                                                gTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGc                 <  2:163628/150‑1 (MQ=255)
                                                                                                                                              tATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatg  >  2:116445/1‑151 (MQ=255)
                                                                                                                                              tATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGgatg  >  2:41537/1‑151 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
CTCTGGGATACATTGCTTTAAATCGTGGGGAACAGATCAACGCGCTGTGGATTGTGGTGGCGTCGGTCTGTATCTATCTGATCGCTTACCGTTTTTATGGGCTGTATATCGCCAAAAATGTGCTGGCGGTTGACCCGACGCGTATGACGCCAGCGGTGCGCCATAACGACGGGCTGGACTATGTGCCGACGGACAAGAAAGTGCTGTTCGGTCACCATTTTGCGGCCATTGCCGGAGCAGGTCCGCTGGTGGGGCCGGTACTGGCGGCGCAAATGGGCTACCTGCCGGGGATG  >  NC_000913/629952‑630244

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