New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NZ_CP009273 = 349144039 (0.970)4 (0.110) 4/152 2.5 9.5% coding (254/1374 nt) cysG siroheme synthase CysG
?NZ_CP009273 = 3491445 41 (1.140)coding (259/1374 nt) cysG siroheme synthase CysG
Rejected: Frequency below/above cutoff threshold.

GTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTC  >  NZ_CP009273/3491357‑3491530
                                                                                        |                                                                                     
gTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCa                                                                                      <  1:630495‑M2/90‑2 (MQ=255)
      ccTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCagcgcgt                                                                                >  1:457490‑M2/1‑83 (MQ=255)
       cTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCagcgcgtc                                                                               <  2:207923‑M2/90‑9 (MQ=255)
        tcgtcgAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCagcgcgtca                                                                              <  1:954568‑M2/90‑10 (MQ=255)
         cgtcgAAGGGCCATTTGATGAAAGCCTTCTCGGCACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAg                                                                             >  1:294995/1‑90 (MQ=255)
         cgtcgAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAg                                                                             >  2:144073/1‑90 (MQ=255)
         cgtcgAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAg                                                                             >  2:1026487/1‑90 (MQ=255)
         cgtcgAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAg                                                                             >  2:418733/1‑90 (MQ=255)
         cgtcgAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAg                                                                             >  1:410833/1‑90 (MQ=255)
          gtcgAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGc                                                                            >  2:125675/1‑90 (MQ=255)
             gAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGaa                                                                         <  2:595047/90‑1 (MQ=255)
              aaGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAg                                                                        <  1:1033484/90‑1 (MQ=255)
                     aTTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGa                                                                 >  2:245883/1‑90 (MQ=255)
                        tgatgaAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGc                                                              <  2:706930/90‑1 (MQ=255)
                              aaGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCtcgtcg                                                        >  2:812926/1‑90 (MQ=255)
                              aaGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCtcgtcg                                                        <  1:1140330/90‑1 (MQ=255)
                                gCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCa                                                      <  2:267429/90‑1 (MQ=255)
                                  cTTCTCGACCCCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTACCCagcgcgtcatcatccgtcgctgcaatcgccagcca                                                    >  1:316920‑M2/1‑55 (MQ=255)
                                  cTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCAtc                                                    <  2:410833/90‑1 (MQ=255)
                                  cTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCAtc                                                    <  1:684906/90‑1 (MQ=255)
                                   ttCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCACTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCAtct                                                   >  2:116870/1‑90 (MQ=255)
                                    tctcGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCAtctt                                                  >  1:921134/1‑90 (MQ=255)
                                     ctcGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCagcgcgtcatcatccgtcgctgcaatcgccagccagca                                                 >  2:510604‑M2/1‑52 (MQ=255)
                                           acCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCa                                                      <  1:158215/79‑1 (MQ=255)
                                           acCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCa                                                      >  2:158215/1‑79 (MQ=255)
                                            ccTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCagcgcg                                                                                 >  2:325145‑M2/1‑45 (MQ=255)
                                            ccTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCagcgcg                                                                                 <  1:325145‑M2/51‑7 (MQ=255)
                                            ccTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACg                                          >  2:267138/1‑90 (MQ=255)
                                              tgctgGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTg                                        >  2:56015/1‑90 (MQ=255)
                                               gctggctgGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTgg                                       >  2:992451/1‑90 (MQ=255)
                                                   gctggcGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGa                                   >  2:760484/1‑90 (MQ=255)
                                                    ctggcGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGAt                                  >  1:659851/1‑90 (MQ=255)
                                                     tggcGATTGCAGCGACGGATGATGACGCGCTTACCCagcgcgtcatcacccgtcgctgcaatcgccagccagcaggtgtcgagaaggctt                                 >  2:239716‑M2/1‑36 (MQ=255)
                                                         gATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGcc                             <  1:1100331/90‑1 (MQ=255)
                                                              cAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCa                                                      >  2:48809/1‑60 (MQ=255)
                                                              cAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCa                                                      <  1:48809/60‑1 (MQ=255)
                                                                  gACGGATGATGCCGCCCTTACCCagcgcgtcatcatccgtggctgcaatcgccagccagcaggtgtcgagaaggctttcatcaaatggcc                    >  2:191460‑M2/1‑23 (MQ=255)
                                                                   aCGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAgccgcc                   <  1:826180/90‑1 (MQ=255)
                                                                   aCGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAgccgcc                   <  1:581535/90‑1 (MQ=255)
                                                                    cGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCa                  <  1:612544/90‑1 (MQ=255)
                                                                      gatgatGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGc                <  2:213254/90‑1 (MQ=255)
                                                                        tgatgaCGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCtt              >  1:181824/1‑90 (MQ=255)
                                                                             aCGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGc                                                              <  2:1096243/37‑1 (MQ=255)
                                                                             aCGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGc                                                              >  1:1096243/1‑37 (MQ=255)
                                                                               gcgcTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATg       >  1:241386/1‑90 (MQ=255)
                                                                                cgcTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGc      <  2:630495/90‑1 (MQ=255)
                                                                                cgcTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGc      <  2:78419/90‑1 (MQ=255)
                                                                                cgcTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGc      <  2:874563/90‑1 (MQ=255)
                                                                                cgcTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGc      <  2:954568/90‑1 (MQ=255)
                                                                                    tAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTc  <  2:1076682/90‑1 (MQ=255)
                                                                                    tAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCACATTTTCTGTAACGTGGTCGATGCGCAGAAAGCCGCCAGCTTTATTATGCCGTc  <  1:224657/90‑1 (MQ=255)
                                                                                        |                                                                                     
GTTAACCCTCGTCGAAGGGCCATTTGATGAAAGCCTTCTCGACACCTGCTGGCTGGCGATTGCAGCGACGGATGATGACGCGCTTAACCAGCGCGTCAGCGAAGCCGCTGAAGCTCGTCGCATCTTCTGTAACGTGGTCGATGCGCCGAAAGCCGCCAGCTTTATTATGCCGTC  >  NZ_CP009273/3491357‑3491530

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