New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913_3_vch_pgi 3645087 =17 (0.930)4 (0.230) 4/274 1.9 19.9% coding (96/2043 nt) prlC oligopeptidase A
?NC_000913_3_vch_pgi 3645120 = 16 (0.920)coding (63/2043 nt) prlC oligopeptidase A
Rejected: Frequency below/above cutoff threshold.

CGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAG  >  NC_000913_3_vch_pgi/3644976‑3645258
                                                                                                                                                |                                                                                                                                          
cGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATggt                                                                                                                                        >  2:188822/1‑147 (MQ=255)
       aagaTACACCCCAACACATCGTCCACTACCGCCAGCCGCCGGCACAGATAGTCCCAGGTGTACGGTGCCCCATACGCTACGACGCGCTCCACAATCTCCCTGCAGCCCGTCAACGCCTTAGCCCCGGCTGAACCGACAAGTCCCggaaa                                                                                                                                 >  2:101805/1‑146 (MQ=255)
                       caTCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGa                                                                                                                 <  2:55415/149‑1 (MQ=255)
                          ctctttaaccaggactgcgcgaatgacgaatccgttactgactccctttgaattgcctccgttttctaaaattctcccggaacatGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGga                                                                                                                                  >  1:21726‑M1/86‑129 (MQ=255)
                          ctctttaaccaggactgcgcgaatgacgaatccgttactgactccctttgaattgcctccgttttctaaaattctcccggaacatGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGga                                                                                                                                  <  2:21726‑M1/44‑1 (MQ=255)
                                                            ccAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTCTAAAAAACGGAGGCAAATCAAAGGGAGTCAGTAACGGCTTCGTCATTc                                                                            >  1:43444/1‑149 (MQ=255)
                                                                       gTGCCCCTTGCGCTAGTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACTGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTg                                                                 <  2:258154/149‑1 (MQ=255)
                                                                                 cgcTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAgg                                                       >  2:146155/1‑149 (MQ=255)
                                                                                  gcTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGt                                                      >  1:306730/1‑149 (MQ=255)
                                                                                           cgcTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAAGTTAAGGAAgc                                             <  1:15261/149‑1 (MQ=255)
                                                                                           cgcTCCACGTTTGCGCGGCAGTCGTTCAATCCCTTAGTCACGGCTGGAACGATATGTGCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACTGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAgc                                             <  1:14365/149‑1 (MQ=255)
                                                                                                    tcccggaacatgtcgttccagccgtgactaaggcattgaacgacATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAagag                                                        >  1:21726‑M2/45‑129 (MQ=255)
                                                                                                    tcccggaacatgtcgttccagccgtgactaaggcattgaacgacATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAagag                                                        <  2:21726‑M2/85‑1 (MQ=255)
                                                                                                              tgtcgttccagccgtgactaaggcattgaacgaca                                                                                                                                            >  1:48128‑M2/35‑35 (MQ=255)
                                                                                                              tgtcgttccagccgtgactaaggcattgaacgaca                                                                                                                                            <  2:48128‑M2/1‑1 (MQ=255)
                                                                                                              tGTCGTTCAATGCCTTAGTCACGGCTGGAACGACa                                                                                                                                            <  2:48128‑M1/34‑1 (MQ=255)
                                                                                                              tGTCGTTCAATGCCTTAGTCACGGCTGGAACGACa                                                                                                                                            >  1:48128‑M1/2‑35 (MQ=255)
                                                                                                                  GTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAt                      <  1:4493‑M1/149‑1 (MQ=255)
                                                                                                                   TTCAATTGCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTCACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAtg                     >  1:108518‑M1/1‑149 (MQ=255)
                                                                                                                      aaTGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAAtgtgt                  <  1:107412/149‑1 (MQ=255)
                                                                                                                           cTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTACTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATCATGTGTGTAGc             >  1:124841/1‑149 (MQ=255)
                                                                                                                                tCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATggg        >  1:195440/1‑149 (MQ=255)
                                                                                                                                 taCGGCTGGAACGATATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTGAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGGTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATgggg       <  2:112190/148‑1 (MQ=255)
                                                                                                                                      cTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAg  >  1:170372/1‑149 (MQ=255)
                                                                                                                                        ggAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTa                            <  1:7485/121‑1 (MQ=255)
                                                                                                                                        ggAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTa                            >  2:7485/1‑121 (MQ=255)
                                                                                                                                                |                                                                                                                                          
CGGGGAGAAGATACGCCCCAACACATCGTCCACTTCCGCCAGCGGCTGGCAGAGATTTTCCCAGGTGTACGGTGCCCCTTGCGCTACTACGCGCTCCACGTTTTCGCGGCAGTCGTTCAATGCCTTAGTCACGGCTGGAACGACATGTTCCGGGAGAATTTTAGAAAACGGAGGCAATTCAAAGGGAGTCAGTAACGGATTCGTCATTCGCGCAGTCCTGGTTAAAGAGGTTAAGGAAGCGCTCAACAGGCGCTTTACATAATGTGTGTAGCATGGGGTTAAG  >  NC_000913_3_vch_pgi/3644976‑3645258

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