New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 9037414 (1.070)3 (0.260) 3/242 1.2 18.0% intergenic (+49/+73) PP_0079/PP_0080 transcriptional regulator/hypothetical protein
?NC_002947 = 90397 15 (1.300)intergenic (+72/+50) PP_0079/PP_0080 transcriptional regulator/hypothetical protein
Rejected: Frequency below/above cutoff threshold.

GCTCGGAACGCGGCTATTACGTAGTCTTACCGCCGCGCAAGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCGCAGCAGCCCCGGCAATCTAGACCACGAAGTTCAGGTGCTCGCTCCGGTGCAAGTCCAGCTCCAGCATGTCCACCATCCCCGCTGCGATCCTGG  >  NC_002947/90225‑90523
                                                                                                                                                     |                                                                                                                                                     
gcTCGGAACGCGGCTATTACGTAGTCTTACCGCCGCGCAAGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACA                                                                                                                                                       >  2:126916‑M1/1‑150 (MQ=255)
gcTCGGAACGCGGCTATTACGTAGTCTTACCGCCGCGCAAGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACA                                                                                                                                                       >  1:189785‑M1/1‑150 (MQ=255)
         gcgGCTATTACGTAGTCTTACCGCCGCGCAAGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTg                                                                                                                                              >  1:269774/1‑150 (MQ=255)
          cgGCTATTACGTAGTCTTACCGCCGCGCAAGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGc                                                                                                                                             <  2:198911/150‑1 (MQ=255)
                        tCTTACCGCCGCGCAAGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCAt                                                                                                                               <  1:306037‑M2/150‑2 (MQ=255)
                          ttACCGCCGCGCAAGCGGCGTGGGGCGTTGATCGAACGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCAtgt                                                                                                                             <  1:231558‑M2/150‑4 (MQ=255)
                               gccgcGCAAGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAtgggccggcgtgtgcagtgccactgtgggag                                                                                                                        <  1:296277‑M1/150‑32 (MQ=255)
                                   cgcAAGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCAtgtgggagcggc                                                                                                                    <  2:206687‑M2/150‑13 (MQ=255)
                                      aaGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCAtgtgggagcggcctt                                                                                                                 <  2:134631‑M2/150‑16 (MQ=255)
                                                gggCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGAt                                                                                                       <  1:126916/150‑1 (MQ=255)
                                                     ttGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCt                                                                                                  <  2:269774/150‑1 (MQ=255)
                                                        aTCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAtgggccggcgtgtgcagtgccactgtggga                                                                                                                         >  2:25258‑M1/1‑94 (MQ=255)
                                                        aTCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAtgggccggcgtgtgcagtgccactgtggga                                                                                                                         <  1:25258‑M1/124‑31 (MQ=255)
                                                        aTCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCAtgtgggagcggcct                                                                                                                  >  1:52314‑M2/1‑117 (MQ=255)
                                                        aTCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCAcgtgggagcggcct                                                                                                                  <  2:52314‑M2/131‑15 (MQ=255)
                                                           gAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCgcag                                                                                            >  2:131636/1‑150 (MQ=255)
                                                            aGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCgcagc                                                                                           <  2:189785/150‑1 (MQ=255)
                                                             gCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCgcagca                                                                                          <  1:131636/150‑1 (MQ=255)
                                                                                gAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCGCAGCAGCCCCGGCAATCTAGACCa                                                                       >  1:84085/1‑150 (MQ=255)
                                                                                          gcAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCGCAGCAGCACCGGCAATCTAGACCACGAAGTTCAg                                                             <  1:134631/150‑1 (MQ=255)
                                                                                                      aTTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCGCAGCAGCCCCGGCAATCTAGACCACGAAGTTCAGGTGCTCGCTCCg                                                 <  2:306037/150‑1 (MQ=255)
                                                                                                                      gCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCGCAGCAGCCCCGGCAATCTAGACCACGAAGTTCAGGTGCTCGCTCCGGTGCAAGtccagctcc                                 <  2:8621/150‑1 (MQ=255)
                                                                                                                              tCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCGCATTTGGGAGCGGCCTTGCGTCGCGATGGGCTGCGCAGCAGCCCCGGCAATCTAGACCACGAAGTTCAGGTGCTCGCTCCGGTGCAAGTCCAGCTCCAGCATGTc                         >  1:151141/1‑150 (MQ=255)
                                                                                                                                       aaGGCCGCTCCCACAtgggccggcgtgtgcagtgccactgtgggagcggccttgtgtcgcgatgggccgcaaagcggccccaaaatctcaaaggctgcgctcctgctccaaccaatccacaaaccgctcgatcaacgccccacgccgctt                <  2:134631‑M1/150‑136 (MQ=255)
                                                                                                                                        aGGCCGCTCCCACGtgggccggcgtgtgcagtgccactgtgggagcggccttgtgtcgcgatgggccgcaaagcggccccaaaatctcaaaggctgcgctcctgctccaaccaatccacaaaccgctcgat                                  <  2:52314‑M1/131‑118 (MQ=255)
                                                                                                                                        aGGCCGCTCCCACAtgggccggcgtgtgcagtgccactgtgggagcggccttgtgtcgcgatgggccgcaaagcggccccaaaatctcaaaggctgcgctcctgctccaaccaatccacaaaccgctcgat                                  >  1:52314‑M1/1‑14 (MQ=255)
                                                                                                                                          gCCGCTCCCACAtgggccggcgtgtgcagtgccactgtgggagcggccttgtgtcgcgatgggccgcaaagcggccccaaaatctcaaaggctgcgctcctgctccaaccaatccacaaaccgctcgatcaacgccccacgccgcttgcg             <  2:206687‑M1/150‑139 (MQ=255)
                                                                                                                                              cTCCCACAGTGGCACTGCACACGCCGGCCCAtgtgggagcggccttgtgtcgcgatgggccgcaaagcggccccaaaatctcaaaggctgcgctcctgctccaaccaatccacaaaccgctcgatcaacgccccacgccgcttgcgcggc         <  1:296277‑M2/150‑120 (MQ=255)
                                                                                                                                               tCCCACAGTGGCACTGCACACGCCGGCCCAtgtgggagcggccttgtgtcgcgatgggccgcaaagcggccccaaaatctcaaaggctgcgctcctgctccaaccaatccacaaaccgctcgat                                  <  1:25258‑M2/124‑95 (MQ=255)
                                                                                                                                               tCCCACAGTGGCACTGCACACGCCGGCCCAtgtgggagcggccttgtgtcgcgatgggccgcaaagcggccccaaaatctcaaaggctgcgctcctgctccaaccaatccacaaaccgctcgat                                  >  2:25258‑M2/1‑30 (MQ=255)
                                                                                                                                                 ccACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCGCAGCAGCCCCGGCAATCTAGACCACGAAGTTCAGGTGCTCGCTCCGGTGCAAGtc                                        >  1:232711/1‑116 (MQ=255)
                                                                                                                                                 ccACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCGCAGCAGCCCCGGCAATCTAGACCACGAAGTTCAGGTGCTCGCTCCGGTGCAAGtc                                        <  2:232711/116‑1 (MQ=255)
                                                                                                                                                   acatgggccggcgtgtgcagtgccactgtgggagcggccttgtgtcgcgatgggccgcaaagcggccccaaaatctcaaaggctgcgctcctgctccaaccaatccacaaaccgttcgatcaacgccccacgccgcttgcgcggcggtaa    <  1:231558‑M1/150‑148 (MQ=255)
                                                                                                                                                     atgggccggcgtgtgcagtgccactgtgggagcggccttgtgtcgcgatgggccgcaaagcggccccaaaatctcaaaggctgcgctcctgctccaaccaatccacaaaccgctcgatcaacgccccacgccgcttgcgcggcggtaaga  <  1:306037‑M1/150‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCGGAACGCGGCTATTACGTAGTCTTACCGCCGCGCAAGCGGCGTGGGGCGTTGATCGAGCGGTTTGTGGATTGGTTGGAGCAGGAGCGCAGCCTTTGAGATTTTGGGGCCGCTTTGCGGCCCATCGCGACACAAGGCCGCTCCCACAGTGGCACTGCACACGCCGGCCCATGTGGGAGCGGCCTTGCGTCGCGATGGGCTGCGCAGCAGCCCCGGCAATCTAGACCACGAAGTTCAGGTGCTCGCTCCGGTGCAAGTCCAGCTCCAGCATGTCCACCATCCCCGCTGCGATCCTGG  >  NC_002947/90225‑90523

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