New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 3000957 =51 (0.930)5 (0.110) 5/134 2.9 11.7% coding (89/1197 nt) dpaL diaminopropionate ammonia‑lyase
?NZ_CP009273 3000973 = 34 (0.760)coding (105/1197 nt) dpaL diaminopropionate ammonia‑lyase
Rejected: Frequency below/above cutoff threshold.

TTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGAAAATTCTC  >  NZ_CP009273/3000884‑3001057
                                                                                         |                                                                                    
ttGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTc                                                                                      >  1:895193/1‑90 (MQ=255)
    aGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTccacc                                                                                  >  1:14138/1‑90 (MQ=255)
      aTTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAg                                                                                <  1:1602903/90‑1 (MQ=255)
         gATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGaaa                                                                             >  1:1514858/1‑90 (MQ=255)
            atCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATa                                                                          >  1:232563/1‑90 (MQ=255)
            atCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGACAGTTCCACCAGAAAATa                                                                          >  1:329366/1‑90 (MQ=255)
              ttgcgagatcgtccagcgcacaaagcggtgttgggcgataaccagctattttctggtggAACTGGCGCGCCAGTTCCACCAGaaa                                                                             >  1:1169045‑M1/60‑85 (MQ=255)
              ttgcgagatcgtccagcgcacaaagcggtgttgggcgataaccagctattttctggtggAACTGGCGCGCCAGTTCCACCAGaaa                                                                             <  2:1169045‑M1/26‑1 (MQ=255)
               gCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCt                                                                       >  1:489787/1‑90 (MQ=255)
                ccGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTg                                                                      <  2:1173142/90‑1 (MQ=255)
                 cGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTgg                                                                     >  1:1333577/1‑90 (MQ=255)
                   aTAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGtt                                                                   <  1:769377/90‑1 (MQ=255)
                   aTAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGtt                                                                   <  1:1572393/90‑1 (MQ=255)
                     atcgtccagcacacaaagcggtgttgggcgataaccagctattttctggtggAACTGGCGCGCCCGTTACACCAGAAAATAGCTGGTTAt                                                                 >  2:62123‑M1/53‑90 (MQ=255)
                     aaCAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTAt                                                                 <  1:980036/90‑1 (MQ=255)
                     aaCAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTAt                                                                 <  2:1333577/90‑1 (MQ=255)
                     aaCAAATTTTTCAACGGCGAAACATCACCGCCCTTTTCGCAAAGCCAGTCCAAACTGGCGCGCCAGTACCACCAGAAAATAGCTGTTTAt                                                                 <  2:769119/90‑1 (MQ=255)
                         aaTTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGcc                                                             >  2:346756/1‑90 (MQ=255)
                          ccagcgcacaaagcggtgttgggcgataaccagctattttctggtggaactggcgcgccagttCCCCCAGAAAATAGCTGGTTATCGccc                                                            <  2:326495‑M2/27‑1 (MQ=255)
                          aTTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGccc                                                            <  1:1031456/90‑1 (MQ=255)
                          aTTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGccc                                                            <  1:1324535/90‑1 (MQ=255)
                          aTTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGccc                                                            <  2:744650/90‑1 (MQ=255)
                                    aagcggtgttgggcgataaccagctattttctggtgtGACGGGAGCGCCAGTTCCACCAGAAAAAAGATGGTTATCGCCCAACaccgttt                                                  <  2:45041‑M1/53‑4 (MQ=255)
                                    aaacggtgttgggcgataaccatcttttttctggtggaactggcgctcccgtcACACCAGAAAATAGCTGGTTATCGCCCAACACCGCtt                                                  <  2:45041‑M2/37‑1 (MQ=255)
                                     gCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCttt                                                 >  1:775397/1‑90 (MQ=255)
                                         aaCATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTgc                                             <  1:967734/90‑1 (MQ=255)
                                         aaCATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTgc                                             <  2:674225/90‑1 (MQ=255)
                                          aCATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTgcg                                            <  1:714436/90‑1 (MQ=255)
                                              gggcgataaccagctattttctgggggAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTgg                                        <  2:326495‑M1/63‑1 (MQ=255)
                                              cACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTgg                                        >  1:672757/1‑90 (MQ=255)
                                               aCCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGa                                       >  1:457384/1‑90 (MQ=255)
                                                  gCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGa                                    >  2:44483/1‑90 (MQ=255)
                                                  gCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGa                                    >  1:1045141/1‑90 (MQ=255)
                                                   ctctTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGAt                                   <  2:369839/90‑1 (MQ=255)
                                                   ataaccagctattttctggtgtaacgggcgcgccagttCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGTGCTGGACGAt                                   >  2:62123‑M2/39‑90 (MQ=255)
                                                       tttCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCg                               >  2:1564183/1‑90 (MQ=255)
                                                          cGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCaa                            <  2:996259/90‑1 (MQ=255)
                                                          cGCAAAGCCAGGCCAAACTGGCGCGACAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCaa                            <  2:329366/90‑1 (MQ=255)
                                                           gCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAAc                          <  2:868425/90‑1 (MQ=255)
                                                               tttctggtggaactggcgcgccagttCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCaa                            <  2:1169045‑M2/59‑1 (MQ=255)
                                                               tttctggtggaactggcgcgccagttCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCaa                            >  1:1169045‑M2/27‑85 (MQ=255)
                                                                  cAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCttttt                    <  1:135767/90‑1 (MQ=255)
                                                                   aGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTg                   >  2:1405223/1‑90 (MQ=255)
                                                                   aGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTg                   >  1:152645/1‑90 (MQ=255)
                                                                         AACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTga             <  1:869691‑M1/90‑1 (MQ=255)
                                                                             GGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTgaagaa         >  2:1426152‑M1/1‑90 (MQ=255)
                                                                             GGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTgaagaa         >  1:1617244‑M1/1‑90 (MQ=255)
                                                                              gcgcgccagttCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGaaa        >  1:19397‑M2/12‑90 (MQ=255)
                                                                                   ccagttCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGAAAATtct   <  1:1175550‑M2/84‑1 (MQ=255)
                                                                                    cagttCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGAAAATtctc  <  2:199427‑M2/85‑1 (MQ=255)
                                                                                         |                                                                                    
TTGAAGATTGATATCGCCGATAACAAATTTTTCAACGGCGAAACATCACCGCTCTTTTCGCAAAGCCAGGCCAAACTGGCGCGCCAGTTCCACCAGAAAATAGCTGGTTATCGCCCAACACCGCTTTGTGCGCTGGACGATCTCGCAAACCTTTTTGGTGTGAAGAAAATTCTC  >  NZ_CP009273/3000884‑3001057

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