New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 1040657 =20 (0.570)4 (0.130) 4/150 2.4 15.1% coding (603/1140 nt) gfcE polysaccharide export protein
?NZ_CP009273 1040766 = 27 (0.850)coding (494/1140 nt) gfcE polysaccharide export protein
Rejected: Frequency below/above cutoff threshold.

TGATGGTCAGAGGAATATTGGTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTGGGTAAGTGTTTTCCCGGC  >  NZ_CP009273/1040677‑1040851
                                                                                         |                                                                                     
tgatgGTCAGAGGAATATTGGTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGa                                                                                       >  1:1076262/1‑90 (MQ=255)
   tgGTCAGAGGAATATTGGTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGc                                                                                    <  1:435767/90‑1 (MQ=255)
    gtacgtcaccggcgaagtcagcaaatcgggtcagcagccgattaccaatattcctctgaccatcatggatgcgattaacgccgcgACGCt                                                                                   <  2:226490‑M2/5‑1 (MQ=255)
    ggTCAGAGGAATATTGGTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCt                                                                                   >  1:243873/1‑90 (MQ=255)
    ggTCAGAGGAATATTGGTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCt                                                                                   >  1:89735/1‑90 (MQ=255)
         gagGAATATTGGTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACAt                                                                              >  2:625868/1‑90 (MQ=255)
              atatTGGTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACt                                                                         >  2:313982/1‑90 (MQ=255)
                 ttGGTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGc                                                                      <  2:133007/90‑1 (MQ=255)
                  aagtcagcaaatcgggtcagcagccgattaccaatattcctctgaccatcatggatgcgattaacgccgcgACGCTGACATCGACTTGCg                                                                     <  1:372536‑M2/19‑1 (MQ=255)
                   ggTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCgg                                                                    >  2:87399/1‑90 (MQ=255)
                   ggTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCgg                                                                    >  1:849314/1‑90 (MQ=255)
                   ggTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCgg                                                                    >  1:636365/1‑90 (MQ=255)
                   agtcagcaaatcgggtcagcagccgattaccaatattcctctgaccatcatggatgcgattaacgccgcgACGCTGACATCGACTTGCgg                                                                    <  1:441930‑M2/20‑1 (MQ=255)
                     tAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGc                                                                  <  1:402821/90‑1 (MQ=255)
                     tAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGc                                                                  <  2:243873/90‑1 (MQ=255)
                        tCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCttt                                                               >  2:1086944/1‑90 (MQ=255)
                            ctgctgACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAAt                                                           <  2:1128979/90‑1 (MQ=255)
                               gggtcagcagccgattaccaatattcctctgaccatcatggatgcgattaacgccgcgACGCTGACATCGACTTGCGGGCTTTCAATGAc                                                        <  1:757988‑M2/32‑1 (MQ=255)
                                            cTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGAcgc                                           <  2:791404/90‑1 (MQ=255)
                                            cTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGAcgc                                           <  2:1076262/90‑1 (MQ=255)
                                            cTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGAcgc                                           <  1:313982/90‑1 (MQ=255)
                                            cTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGAcgc                                           <  2:210983/90‑1 (MQ=255)
                                               aCTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCt                                        <  1:862690/90‑1 (MQ=255)
                                                           aCGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTg                            <  1:625868/90‑1 (MQ=255)
                                                            cGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGc                           <  2:482369/90‑1 (MQ=255)
                                                                  cTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCt                     <  1:732407/90‑1 (MQ=255)
                                                                    ttttGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTaa                                    <  2:369676/73‑1 (MQ=255)
                                                                    ttttGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTaa                                    >  1:369676/1‑73 (MQ=255)
                                                                    ttttGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTgg                   >  2:71237/1‑90 (MQ=255)
                                                                       tGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTGGGTa                <  2:347818/90‑1 (MQ=255)
                                                                            tCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTGGGTAAgtgt           >  1:396235/1‑90 (MQ=255)
                                                                                    ccgcgACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAAtct                                 <  1:1108953‑M2/55‑1 (MQ=255)
                                                                                    ccgcgACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAAtct                                 >  2:1108953‑M2/6‑60 (MQ=255)
                                                                                     cgcgACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTGGGTAAGTGTTTTCCCGGc  <  1:71237‑M2/86‑1 (MQ=255)
                                                                                         aCGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTGGGTAAGTGTTTTcc      >  2:1106898/1‑82 (MQ=255)
                                                                                         aCGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTGGGTAAGTGTTTTcc      <  1:795658/82‑1 (MQ=255)
                                                                                         aCGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTGGGTAAGTGTTTTcc      >  2:795658/1‑82 (MQ=255)
                                                                                         aCGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTGGGTAAGTGTTTTcc      <  1:1106898/82‑1 (MQ=255)
                                                                                         |                                                                                     
TGATGGTCAGAGGAATATTGGTAATCGGCTGCTGACCCGATTTGCTGACTTCGCCGGTGACGTACGCTTTTTGCGATCGGAACGCCGCGACGCTGACATCGACTTGCGGGCTTTCAATGACTGAATCGAGACGCGCTGTAATCTCGTTGCGTACCTGGGTAAGTGTTTTCCCGGC  >  NZ_CP009273/1040677‑1040851

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.