New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 349387123 (0.880)3 (0.130) 3/138 1.9 13.5% coding (603/1338 nt) frlA fructoselysine/psicoselysine transporter
?NZ_CP009273 = 3494083 18 (0.780)coding (815/1338 nt) frlA fructoselysine/psicoselysine transporter
Rejected: Frequency below/above cutoff threshold.

CTGCTGGCGCTGGTGATTTCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCGATGATCGTCATTCTTGGTTCGCTTTCCAGCTGCGTGATGT  >  NZ_CP009273/3494001‑3494172
                                                                                  |                                                                                         
ctgctgGCGCTGGTGATTTCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGgatcccc                                                                                    <  2:516208‑M2/90‑8 (MQ=255)
ctgctgGCGCTGGTGATTTCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGgatcccc                                                                                    <  1:282619‑M2/90‑8 (MQ=255)
ctgctgGCGCTGGTGATTTCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGgatcccc                                                                                    <  2:346415‑M2/90‑8 (MQ=255)
ctgctgGCGCTGGTGATTTCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGgatcccc                                                                                    <  2:57550‑M2/90‑8 (MQ=255)
ctgctgGCGCTGGTGATTTCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGgatcccc                                                                                    <  2:716206‑M2/90‑8 (MQ=255)
   ctgGCGCTGGTGATTTCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCa                                                                                 >  1:458763/1‑90 (MQ=255)
          tGGTGATTTCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGgca                                                                          >  2:114293/1‑90 (MQ=255)
                 ttCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGGCAGCACCGc                                                                   <  1:516208/90‑1 (MQ=255)
                 ttCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGGCAGCACCGc                                                                   <  1:716206/90‑1 (MQ=255)
                     ggCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGGCAGCACCGCTggg                                                               >  2:864232/1‑90 (MQ=255)
                     ggCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGGCAGCACCGCTggg                                                               >  2:699209/1‑90 (MQ=255)
                                          cTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCg                                          >  2:754371/1‑90 (MQ=255)
                                           tCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCga                                         <  2:458763/90‑1 (MQ=255)
                                                            ccGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGGCAGCAc                                                                      <  1:685878/44‑1 (MQ=255)
                                                            ccGATTTCCGACGCCCTGACCGGGATCCCCGCACTCGGCAGCAc                                                                      >  2:685878/1‑44 (MQ=255)
                                                                    cGACGCCCTGACCTGgatccccgccagcagcgccataaagctgcccgttgcgccaatcgcag                                            <  2:303136‑M2/62‑48 (MQ=255)
                                                                    cGACGCCCTGACCTGgatccccgccagcagcgccataaagctgcccgttgcgccaatcgcag                                            >  1:303136‑M2/1‑15 (MQ=255)
                                                                     gACGCCCTGACCTGGATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCGATGATCGTCATTCTTGGTTCGCTTTcc               <  2:282619/90‑1 (MQ=255)
                                                                        gCCCTGACCTGGATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCGATGATCGTCATTCTTGGTTCGCTTTCCAGc            >  2:516214/1‑90 (MQ=255)
                                                                          ccTGACCTGGATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCGATGATCGTCATTCTTGGTTCGCTTTCCAGCTg          <  1:21199/90‑1 (MQ=255)
                                                                                agggatccccgccagcagcgccataaagctgcccgttgcgccaatcgcagtggtggtaggggcggcaaaattctccgctttaaaccagaa    >  1:455484‑M2/3‑3 (MQ=255)
                                                                                  ggATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCgatg                                       >  1:117816/1‑53 (MQ=255)
                                                                                  ggATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCgatg                                       <  2:117816/53‑1 (MQ=255)
                                                                                  ggATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCGATGATCGTCATTCTTGGTTCGCTTTCCAGCTGCGTGATGt  >  1:68521/1‑90 (MQ=255)
                                                                                  ggATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCGATGATCGTCATTCTTGGTTCGCTTTCCAGCTGCGTGATGt  >  1:537778/1‑90 (MQ=255)
                                                                                  ggATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCGATGATCGTCATTCTTGGTTCGCTTTCCAGCTGCGTGATGt  >  2:764512/1‑90 (MQ=255)
                                                                                  |                                                                                         
CTGCTGGCGCTGGTGATTTCCGGCCTGATGCCCTTCGACAAACTCGCCAATTCTGAAACGCCGATTTCCGACGCCCTGACCTGGATCCCCGCACTCGGCAGCACCGCTGGGATCTTTGTTGCCATCACGGCGATGATCGTCATTCTTGGTTCGCTTTCCAGCTGCGTGATGT  >  NZ_CP009273/3494001‑3494172

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