New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 3483915 =11 (0.570)5 (0.270) 5/250 NT 37.5% coding (2052/3804 nt) PP_3091 membrane protein
?NC_002947 3483952 = 6 (0.320)coding (2015/3804 nt) PP_3091 membrane protein

CGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCAC  >  NC_002947/3483766‑3484055
                                                                                                                                                     |                                                                                                                                            
ctgtcactgtccgaagaaatgctcgggcgcctggacgagcagcaggtggcctcgatcatcgaggcctacaacaccctcaagctgtacctgttgctgaccgagccccaggcccatccagacccggccttcgtcgctgccagcctgccccag                                                                                                                                              >  2:115653‑M1/150‑150 (MQ=255)
      aCCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCggg                                                                                                                                        >  2:452343/1‑150 (MQ=255)
        cAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGt                                                                                                                                      >  2:186523/1‑150 (MQ=255)
                      tcgggcgcctggacgagcagcaggtggcctcgatcatcgaggcctacaacaccctcaagctgtacctgttgctgaccgagccccaggcccatccagacccggccttcgtcgctgccagcctgccccaGGCCGGGGTGCCTTCGCCTgcgg                                                                                                                        <  1:100861‑M1/23‑1 (MQ=255)
                          gAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTTCTCCAGCAACTGCACGTACAGCGTGGGATTTTTTAAGGCGTGGCCATCGTCGGCCGGGTTGCCTTCGCCTTCGGctct                                                                                                                    >  1:252132/1‑150 (MQ=255)
                                 ttGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAg                                                                                                             <  2:404518/150‑1 (MQ=255)
                                  tGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAgg                                                                                                            <  2:486933/150‑1 (MQ=255)
                                        agcaggtggcctcgatcatcgaggcctacaacaccctcaagctgtacctgttgctgaccgagccccaggcccatccagacccggccttcgtcgctgccagcctgccccaGGCCGTGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTggg                                                                                                      >  1:11855‑M1/110‑150 (MQ=255)
                                              gggTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGCCCCAGGCCTGGGCCAgcg                                                                                                <  1:171041/150‑3 (MQ=255)
                                                       tcatcgaggcctacaacaccctcaagctgtacctgttgctgaccgagccccaggcccatccagacccggccttcgtcgctgccagcctgccccaGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGcga                                                                                       >  2:100861‑M1/95‑150 (MQ=255)
                                                          tcgaggcctacaacaccctcaagctgtacctgttgctgaccgagccccaggcccatccagacccggccttcgtcgctgccagcctgccccaGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTgg                                                                                                       <  1:359918‑M1/40‑1 (MQ=255)
                                                          tcgaggcctacaacaccctcaagctgtacctgttgctgaccgagccccaggcccatccagacccggccttcgtcgctgccagcctgccccaGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTgg                                                                                                       >  2:359918‑M1/92‑131 (MQ=255)
                                                                                         cTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCATGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGt                                                     <  2:242696/150‑1 (MQ=255)
                                                                                         cTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCAGt                                                     <  1:452343/150‑1 (MQ=255)
                                                                                           ggCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCGTGGGGGGCTGGGGGAGGGGTGGGGGAGGCTGGGGGGGGGGGGGGGCGGGTCTGGATGGGCCTGGGGCTCGGTCa                                                   <  1:186523/150‑1 (MQ=255)
                                                                                                  cTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAg                                            >  1:286732/1‑150 (MQ=255)
                                                                                                                                            gCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCAc  >  1:462648/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                            
CGATCGACCAGCGAGCTGGAGATCATGAACGACTTGAGGTTCTGCCGGGTTTCGTTGATCAGCTGTTCGTTGCGCGGCAACGACGGCGCCTGGCCCTGCTCCAGCAACTGCACGTACAGCGGGGCATTTTCTGCGATCACGCCAGCGTCGGCCGGGGTGCCTTCGCCTGCGGCGCTGGCCCAGGCCTGGGGCAGGCTGGCAGCGACGAAGGCCGGGTCTGGATGGGCCTGGGGCTCGGTCAGCAACAGGTACAGCTTGAGGGTGTTGTAGGCCTCGATGATCGAGGCCAC  >  NC_002947/3483766‑3484055

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