New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 3947913 =27 (0.510)7 (0.160) 6/228 5.2 25.6% coding (136/1797 nt) PP_3483 type II secretion system protein
?NC_002947 3947998 = 19 (0.450)coding (51/1797 nt) PP_3483 type II secretion system protein
Rejected: Coverage evenness skew score above cutoff.

TGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCTGTGAGCAGTTGATCGTCGT  >  NC_002947/3947858‑3948137
                                                                                                                                             |                                                                                                                                           
tGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACAACACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACCTCGGTg                                                                                                                                      >  8:278132/1‑148 (MQ=255)
   aGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATACTGTGGGAGCGGGCATGCCCGCGAAGACGTCCGTGCAgg                                                                                                                                  >  2:88638/1‑149 (MQ=255)
                    gaggTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTGCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGGGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCggg                                                                                                                  <  3:30169/148‑1 (MQ=255)
                               ggTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCGCCCACAGGTTCATCACAGGTATGTGGGGTGCATACATGCTGTGGGAGCGGGCCTGGCCGCGAAGACCTCGGTGCAGGCGTTGCAGATTGCGGGGTTGGATCGAAt                                                                                                      >  4:163422/1‑149 (MQ=255)
                                      gtaaaTTTTTGGTGGGTTCTTCGGGGGCATGCCCGCTCCCACAGGTACATCACCGGTCTGTGGGCTGAATACATCCTGTGGGCGGGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTg                                                                                               <  2:27579/147‑1 (MQ=39)
                                                               ggCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATTCATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAggg                                                                      >  1:222997/1‑149 (MQ=255)
                                                                    gCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATCGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGa                                                                                                                         <  1:198120/94‑1 (MQ=255)
                                                                    gCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGa                                                                                                                         >  2:198120/1‑93 (MQ=255)
                                                                    gCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGcgcg                                                                 >  7:247466/1‑149 (MQ=255)
                                                                            ccACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAg                                                           <  5:57038/147‑1 (MQ=255)
                                                                                  gTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTc                                                   >  5:118614/1‑149 (MQ=255)
                                                                                       tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCg                                                                                                                    >  3:73794/1‑79 (MQ=255)
                                                                                       tCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCg                                                                                                                    <  4:73794/79‑1 (MQ=255)
                                                                                          caGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAggg                                           >  4:140070/1‑149 (MQ=255)
                                                                                             gTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGAAGGCTTTGCAGATTACGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTg                                        >  1:243039/1‑149 (MQ=255)
                                                                                                tgtgGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGc                                       <  2:128851/147‑1 (MQ=255)
                                                                                                tgtgGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGc                                       >  1:128851/1‑147 (MQ=255)
                                                                                                    ggCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGgcgc                                 <  3:42310/149‑1 (MQ=255)
                                                                                                              taccTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTccc                       <  1:43141/147‑1 (MQ=255)
                                                                                                              aTCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTccc                       >  2:37950/1‑149 (MQ=255)
                                                                                                                    tgGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCtgt                   <  1:139461/147‑1 (MQ=255)
                                                                                                                       gAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCTGTGAg                >  2:297593/1‑147 (MQ=255)
                                                                                                                        aGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCTGTGAGCAg             <  1:297593/149‑1 (MQ=255)
                                                                                                                          cGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCTGTGAGCAGtt           >  3:40472/1‑149 (MQ=255)
                                                                                                                              cATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCTGTGAGCAGTTGAtc       <  4:88295/149‑1 (MQ=255)
                                                                                                                                    ccGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCTGTGAGCAGTTGAtcgtcgt  >  6:223754/1‑149 (MQ=255)
                                                                                                                                     cgcgAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAgg                                                                       >  7:194839/1‑79 (MQ=255)
                                                                                                                                     cgcgAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAgg                                                                       <  8:194839/79‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                           
TGCAGACGGGCCTGGGCCAGGAGGTCGCGGGGGTAGGCGGAAATTTTTGCTGGGTTCTTCGCGGGCATGCCCGCTCCCACAGGTACATCACAGGTCTGTGGGCTGCATACATCCTGTGGGAGCGGGCATGCCCGCGAAGACGTCGGTGCAGGCTTTGCAGATTCCGGGGTTGCATCGAATGCTTCTGGCATGGCCTGCGACTCCGCTCAGGGGCGCGTCGCAGGCTGGCTCGGCCAGGGGTGCGGGCGCCGCTATTCCCCTGTGAGCAGTTGATCGTCGT  >  NC_002947/3947858‑3948137

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