New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 3749840 =23 (0.820)3 (0.110) 3/232 NT 13.1% coding (60/372 nt) PP_3313 heat shock protein
?NC_002947 3749877 = 18 (0.680)coding (23/372 nt) PP_3313 heat shock protein

GCCTCCAGCTGTTCGGACAACTGCGGCATGATCTCCCCTTCCAGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTGGATGACTATACGCTTGGGCTGCAGCGCTTCCCGCCGCGGCAGGGAAATGCGCA  >  NC_002947/3749694‑3749985
                                                                                                                                                  |                                                                                                                                                 
accgacaaggtcgaggcccactaccgcgacggcgtcctgcgcatttccctgccgcggcgggaagcgctgcagcccaagcgtatagtcatccagtaacagcaaaggtttgaaggaggaatcaccatgaacgacaagcaagcagtcgtgcgc                                                                                                                                                >  1:238559‑M1/147‑150 (MQ=255)
     cAGCTGTTCGGACAACTGCGGCATGATCTCCCCTTCCAGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGgcag                                                                                                                                           <  1:243363/150‑1 (MQ=255)
       gCTGTTCGGACAACTGCGGCATGATCTCCCCTTCCAGCACCAGCGTGTCACCATCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGgcagca                                                                                                                                         <  2:250785/150‑1 (MQ=255)
               gACAACTGCGGCATGATCTCCCCTTCCAGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTgg                                                                                                                                 <  1:557073/150‑1 (MQ=255)
                    ccgCGGCATGATCTCCCCGTCCAGCACCAGCGTGTCACCCTCCACAGTGACCACGAGATGTTCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTc                                                                                                                            <  1:312711/148‑1 (MQ=255)
                         gCATGATCTCCCCTTCCAGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTg                                                                                                                                  >  1:525401/1‑139 (MQ=255)
                         gCATGATCTCCCCTTCCAGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTg                                                                                                                                  <  2:525401/139‑1 (MQ=255)
                             gATCTCCCCTTCCAGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTcgcc                                                                                                                   >  1:216341/1‑150 (MQ=255)
                                 tcCCCTTCCAGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCCAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCACCAGTGTTTGGGTTTCAGTTTCGCCgcgc                                                                                                               >  2:251560/1‑150 (MQ=255)
                                   cccTTCCAGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATGCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCAc                                                                                                             >  1:548528/1‑150 (MQ=255)
                                    ctgcgcatttccctgccgcggcgggaagcgctgcagcccaagcgtatagtcatccagtaacagcaaaggtttgaaggaggaatcaccatgaacgacaagcaagcagtcgtGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACg                                                                                                            >  2:337817‑M1/111‑150 (MQ=255)
                                          aGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACtgct                                                                                                      >  2:37029/1‑150 (MQ=255)
                                                ctgccgcggcgggaagcgctgcagcccaagcgtatagtcatccagtaacagcaaaggtttgaaggaggaatcaccatgaacgacaagcaagcagtcgtGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACt                                                                                                         <  2:70253‑M1/43‑1 (MQ=255)
                                                ctgccgcggcgggaagcgctgcagcccaagcgtatagtcatccagtaacagcaaaggtttgaaggaggaatcaccatgaacgacaagcaagcagtcgtGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACt                                                                                                         >  1:70253‑M1/99‑141 (MQ=255)
                                                 gCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGt                                                                                               >  1:17984/1‑150 (MQ=255)
                                                            tCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTgg                                                                                                                                 >  2:660544/1‑105 (MQ=255)
                                                            tCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTgg                                                                                                                                 <  1:660544/105‑1 (MQ=255)
                                                                    gACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATtcctcc                                                                            <  1:37029/150‑1 (MQ=255)
                                                                        aCGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCa                                                                        >  2:745132/1‑150 (MQ=255)
                                                                               gatCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCACACCttt                                                                 >  2:207198/1‑150 (MQ=255)
                                                                                         ccccccGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTg                                                       <  2:216341/150‑1 (MQ=255)
                                                                                                         gCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTGGAt                                                    >  1:79476/1‑137 (MQ=255)
                                                                                                         gCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTGGAt                                                    <  2:79476/137‑1 (MQ=255)
                                                                                                                  tGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTGGATGACTATACGCTTGGGCTGCAgc                              <  2:548528/150‑1 (MQ=255)
                                                                                                                      ttCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTGGATGACTATACGCTTGGGCTGCAGCGCtt                          >  2:452201/1‑150 (MQ=255)
                                                                                                                              tCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTGGATGACTATACGCTTGGGCTGCAGCGCTTCCcgccgc                  <  1:745132/150‑1 (MQ=255)
                                                                                                                                            aTCGACGCGCAGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTGGATGACTATACGCTTGGGCTGCAGCGCTTCCCGCCGCGGCAGGGAAATgcg    >  2:627458/1‑150 (MQ=255)
                                                                                                                                              cGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTGGATGACTATACGCTTGGGCTGCAGCGCTTCCCGCCGCGGCAGGGAAATGCGCa  <  2:17984/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCCTCCAGCTGTTCGGACAACTGCGGCATGATCTCCCCTTCCAGCACCAGCGTGTCACCCTCCACATTGACCACGAGATGATCCTTGGGCCCCCCGGGAAAGTCGGCAAACAGCTGGATTCCGTTTTCGTCTTCGAAGACATCGACGCGCGGCAGCAGTGTTTGGGTTTCAGTTTCGCCGCGCACGACTGCTTGCTTGTCGTTCATGGTGATTCCTCCTTCAAACCTTTGCTGTTACTGGATGACTATACGCTTGGGCTGCAGCGCTTCCCGCCGCGGCAGGGAAATGCGCA  >  NC_002947/3749694‑3749985

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