New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 460985734 (1.020)4 (0.130) 4/268 3.3 11.5% coding (344/876 nt) PP_4081 type IV/VI secretion system protein
?NC_002947 = 4609877 30 (0.970)coding (364/876 nt) PP_4081 type IV/VI secretion system protein
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

CTTGTTCGGGCTGGTAATGCGCCTGCGCACGCTCGATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCAGCATCTTCCACGACGAGACCTGGGGTGGGGAA  >  NC_002947/4609744‑4610008
                                                                                                                                     |                                                                                                                                   
cTTGTTCGGGCTGGTAATGCGCCTGCGCACGCTCGATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTc                                                                                                                        >  6:52329/1‑147 (MQ=255)
 ttGTTCGGGCTGGTAATGCGCCTGCGCACGCTCGATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCt                                                                                                                      >  2:98200/1‑148 (MQ=255)
           tGGTAATGCGCCTGCGCACGCTCGATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGgcctgtg                                                                                                            >  6:87389/1‑148 (MQ=255)
               aaTGCGCCTGCGCACGCTCGATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGgcctgtgcct                                                                                                         >  7:88594/1‑147 (MQ=255)
                          gcACGCTCGATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTTCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGa                                                                                             <  5:129370/148‑1 (MQ=255)
                          gcACGCTCGATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGAt                                                                                               <  1:98200/146‑1 (MQ=255)
                                tCGATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCgg                                                                                        >  6:6203/1‑147 (MQ=255)
                                  gATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGa                                                                                     >  2:110819/1‑148 (MQ=255)
                                                cATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCTTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTgg                                                                       <  5:52329/148‑1 (MQ=255)
                                                cATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTgg                                                                       <  6:94338/148‑1 (MQ=255)
                                                  tCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGAATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGATGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTggg                                                                      >  2:174058/1‑147 (MQ=255)
                                                   cAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGc                                                                    <  1:31394/148‑1 (MQ=255)
                                                             cATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGtt                                                          <  5:87389/148‑1 (MQ=255)
                                                              aTCAGCAAGTGTGTAGTCAGCTCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGGGATGGAGCGGCCCGGGGGCAAAAGCAGTTg                                                         >  4:7388/1‑148 (MQ=21)
                                                              aTCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTg                                                         >  5:68133/1‑148 (MQ=255)
                                                                   cAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGa                                                    <  7:73765/148‑1 (MQ=255)
                                                                           tATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAg                                            >  7:145643/1‑148 (MQ=255)
                                                                             ttCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGcc                                          >  8:94715/1‑148 (MQ=255)
                                                                                 gATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTGGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCgct                                       <  8:83606/147‑1 (MQ=255)
                                                                                      aCAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCa                                  <  4:26127/147‑1 (MQ=255)
                                                                                       cAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCAGc                                <  5:169486/148‑1 (MQ=255)
                                                                                          cATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCAGCATc                             >  1:18558/1‑148 (MQ=255)
                                                                                           aTCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCAGCATCt                            <  8:98794/148‑1 (MQ=255)
                                                                                                        tACGCCAGCACGGCTATGAGCCCGCGCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACCTCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGGTGCTGGAGCCACGAGCCGCTGCTCAGCATCTTCCACGACGAGAc               <  5:20596/148‑1 (MQ=255)
                                                                                                           gCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTCCCTCGATGCAGCGGTGATGGAGCGGCGCTGGGGCAACATCAGTTGCTGGAGCCACGCGCCGCTGCTCGGCCTCTTCCACGGCGAGACCtg            >  6:165402/1‑148 (MQ=17)
                                                                                                           gCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGGGCCACGAGCCGCTGCTCAGCATCTTCCACGACGAGACCtg            >  3:52648/1‑148 (MQ=255)
                                                                                                           gCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCAGCATCTTCCACGACGAGACCtg            <  6:68133/148‑1 (MQ=255)
                                                                                                            ccAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCAGCATCTTCCACGACGAGACCtgg           <  3:138141/148‑1 (MQ=255)
                                                                                                                 aCGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCAGCATCTTCCACGACGAGACCtggggtgg      <  8:88594/148‑1 (MQ=255)
                                                                                                                      tATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCAGCATCTTCCACGACGAGACCTGGGGTGGGGaa  >  4:44249/1‑147 (MQ=255)
                                                                                                                                     |                                                                                                                                   
CTTGTTCGGGCTGGTAATGCGCCTGCGCACGCTCGATGAACTGCCCAACATCAAAGATGTGCATCAGCAAGTGCGTATTCAGATCGACAGCATCCGGGAAGAAATACGCCAGCACGGCTATGAGCCCGCTCAGTTGCTGGCGTACTCCTACGGCCTGTGCCTGTACATCGATGAAGCGGTGATGGAGCGGCCCTGGGGCAAAAGCAGTTGCTGGAGCCACGAGCCGCTGCTCAGCATCTTCCACGACGAGACCTGGGGTGGGGAA  >  NC_002947/4609744‑4610008

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