New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 4318388 =20 (1.240)5 (0.320) 4/246 NT 23.7% coding (605/1221 nt) PP_3789 efflux transporter
?NC_002947 4318412 = 13 (0.840)coding (629/1221 nt) PP_3789 efflux transporter

GCTCGGGGTCCGGTTGGCGGAAGTCTCACTTACCCTGGCGCTTTTGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTGGGGTCATTCTTTAACACCGTCGTTCACGGTCGCTTCACGTTCTTCCTGTCCCTGTGGTTGCTGTATCAGTACCCTG  >  NC_002947/4318254‑4318531
                                                                                                                                      |                                                                                                                                                
gCTCGGGGTCCGGTTGGCGGAAGTCTCACTTACCCTGGCGCTTTTGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGc                                                                                                                                  >  2:79359/1‑150 (MQ=255)
    ggggGCCGGTTGGCGGAAGTCTCACTTACCCTGGCGCTTTTGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGAc                                                                                                                              >  2:422457/1‑150 (MQ=255)
            gTTGGCGGAAGTCTCACTTACCCTGGCGCTTTTGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGaa                                                                                                                      >  2:389510/1‑150 (MQ=255)
                  ggAAGTCTCACTTACCCTGGCGCTTTTGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCtt                                                                                                                >  1:266573/1‑150 (MQ=255)
                          acagcaaccacagggacaggaagaacgtgaagcgaccgtgaacgacggtgttaaagaatgaccccaaggtatagagaaccagatttcgatcgttgcccaagactttcaGGGCGCCCTCATTGGCGGACGcc                                                                                                                           <  1:183478‑M1/23‑1 (MQ=255)
                          acagcaaccacagggacaggaagaacgtgaagcgaccgtgaacgacggtgttaaagaatgaccccaaggtatagagaaccagatttcgatcgttgcccaagactttcaGGGCGCCCTCATTGGCGGACGcc                                                                                                                           >  2:183478‑M1/109‑131 (MQ=255)
                               aCCCTGGCGCTTTTGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGa                                                                                                   >  2:96661/1‑150 (MQ=255)
                                   tGGCGCTTTTGGTGGCAGCGGGTGCCTGCGGGGCAGCGATNGTGGTGNTGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATc                                                                                               >  1:352609/1‑150 (MQ=255)
                                         ttttGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGaaa                                                                                                                     <  1:137814/122‑1 (MQ=255)
                                         ttttGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGaaa                                                                                                                     >  2:138003/1‑122 (MQ=255)
                                         ttttGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGaaa                                                                                                                     >  2:137814/1‑122 (MQ=255)
                                         ttttGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGaaa                                                                                                                     <  1:138003/122‑1 (MQ=255)
                                         ttttGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTtc                                                                                         >  2:36454/1‑150 (MQ=255)
                                           ttGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTtctc                                                                                       <  2:266573/150‑1 (MQ=255)
                                                      gggTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTggg                                                                            <  1:36454/150‑1 (MQ=255)
                                                       ggTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTgggg                                                                           >  1:439566/1‑150 (MQ=255)
                                                         tGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTGGGGTc                                                                         <  1:96661/150‑1 (MQ=255)
                                                                   gCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAggg                                                                                                                                                 >  2:115272/1‑69 (MQ=255)
                                                                   gCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAggg                                                                                                                                                 <  1:115272/69‑1 (MQ=255)
                                                                                 gatgatGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTGGGGTCATTCTTTAGCACCGCCGTTCACgg                                                 >  2:309169/1‑150 (MQ=255)
                                                                                                cgcAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTggg                                                                            >  1:49739/1‑108 (MQ=255)
                                                                                                cgcAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTggg                                                                            <  2:49739/108‑1 (MQ=255)
                                                                                                      aaccagatttcgatcgttgcccaagactttcaGGGCGCCCTCATTGGCGGACGcc                                                                                                                           <  2:348422‑M1/23‑1 (MQ=255)
                                                                                                      aaccagatttcgatcgttgcccaagactttcaGGGCGCCCTCATTGGCGGACGcc                                                                                                                           >  1:348422‑M1/33‑55 (MQ=255)
                                                                                                       cAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTGGGGTCATTCTTTAACACCGTCGTTCACGGTCGCTTCACGTTCTTCCTGTcc                           <  1:389510/150‑1 (MQ=255)
                                                                                                                          ccaagactttcaGGGCGTCCGCCAATGAGGGCGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTGGGGTCATTCTTTAACACCGTCGTTCACGGTCGCTTCACGTTCTTCCTGTCCCTGTGGTTGCTGTATCAGt        <  2:243480‑M1/138‑1 (MQ=255)
                                                                                                                                aGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTGGGGTCATTCTTTAACACCGTCGTTCACGGTCGCTTCACGTTCTTCCTGTCCCTGTGGTTGCTGTATCAGTACCCTg  <  1:309169/150‑1 (MQ=255)
                                                                                                                                      |                                                                                                                                                
GCTCGGGGTCCGGTTGGCGGAAGTCTCACTTACCCTGGCGCTTTTGGTGGCAGCGGGTGCCTGCGGGGCAGCGATGGTGGTGATGATGTGCCTGGGCGCAAGTCAGGGCCAACCTCGCACCTCGCTGAAGACAGGGGCGCCCTCATTGGCGGACGCCCTGAAAGTCTTGGGCAACGATCGAAATCTGGTTCTCTATACCTTGGGGTCATTCTTTAACACCGTCGTTCACGGTCGCTTCACGTTCTTCCTGTCCCTGTGGTTGCTGTATCAGTACCCTG  >  NC_002947/4318254‑4318531

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