New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 185974227 (1.000)7 (0.280) 5/270 2.4 19.0% intergenic (+73/‑238) PP_1658/PP_1659 hypothetical protein/hypothetical protein
?NC_002947 = 1859773 35 (1.410)intergenic (+104/‑207) PP_1658/PP_1659 hypothetical protein/hypothetical protein
Rejected: Frequency below/above cutoff threshold.

AACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTCCGGCACAACGGCTAGCCTCCAG  >  NC_002947/1859595‑1859890
                                                                                                                                                   |                                                                                                                                                    
aaCAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACgc                                                                                                                                                    >  2:134742‑M1/1‑148 (MQ=255)
  cAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACccgc                                                                                                                                                  <  1:402134‑M1/150‑5 (MQ=255)
  cAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACAcgc                                                                                                                                                  >  1:329713/1‑150 (MQ=255)
  cAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACAcgc                                                                                                                                                  >  1:482919/1‑150 (MQ=255)
  cAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACAcgc                                                                                                                                                  >  1:482900/1‑150 (MQ=255)
       aTGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGtt                                                                                                                                             <  1:1812/150‑1 (MQ=255)
                     gACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCt                                                                                                                               >  2:455476/1‑150 (MQ=255)
                        gAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGc                                                                                                                            >  2:588605/1‑150 (MQ=255)
                                  cGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCgtcgc                                                                                                                  <  2:463243‑M2/150‑6 (MQ=255)
                                      cgcgGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACgcaaggccagcgtcttctagagaacgcgcgtgtcgcgaaa                                                                                                              <  1:576193‑M1/150‑41 (MQ=255)
                                      cgcgGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCgtcgcgaaa                                                                                                              >  1:406710‑M2/1‑141 (MQ=255)
                                          gCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGt                                                                                                          >  1:372591/1‑150 (MQ=255)
                                              aCAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACgcaaggccagcgtcttctagagaacgcgcgtgtcgcgaaagggccgca                                                                                                      <  1:134742‑M1/150‑49 (MQ=255)
                                                      ccTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACgcaaggccagcgtcttctagagaacgcgcgtgtcgc                                                                                                                  <  1:445988‑M1/130‑37 (MQ=255)
                                                      ccTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACgcaaggccagcgtcttctagagaacgcgcgtgtcgc                                                                                                                  >  2:445988‑M1/1‑94 (MQ=255)
                                                      ccTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACgcaaggccagcgtcttctagagaacgcgcgtgtcgc                                                                                                                  >  2:445980‑M1/1‑94 (MQ=255)
                                                      ccTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACgcaaggccagcgtcttctagagaacgcgcgtgtcgc                                                                                                                  <  1:445980‑M1/130‑37 (MQ=255)
                                                        tgatgaTCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGa                                                                                            >  1:546369/1‑150 (MQ=255)
                                                              tCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGgc                                                                                      <  2:69815/150‑1 (MQ=255)
                                                                    gtgGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGtt                                                                                >  1:687/1‑150 (MQ=255)
                                                                     tgGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGttt                                                                               <  2:101587/150‑1 (MQ=255)
                                                                             aCGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTg                                                                       <  1:427628/150‑1 (MQ=255)
                                                                             aCGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTg                                                                       <  1:595220/150‑1 (MQ=255)
                                                                             aCGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTg                                                                       <  1:427620/150‑1 (MQ=255)
                                                                               gAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCa                                                                     >  1:556258/1‑150 (MQ=255)
                                                                                         ccTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGACGCAGTTTGTCAGTTGCACTATCGTCtt                                                           <  1:603241/150‑1 (MQ=255)
                                                                                           ttGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTgg                                                         >  2:243249/1‑150 (MQ=255)
                                                                                                                      gaaaGGGCCGCAGAGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTa                              >  1:293810/4‑150 (MQ=255)
                                                                                                                      cGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTa                              <  1:588605/150‑1 (MQ=255)
                                                                                                                            gCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTc                        >  2:8270/1‑150 (MQ=255)
                                                                                                                               gCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCAGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCAc                                                                    <  2:46452/103‑1 (MQ=255)
                                                                                                                               gCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCAGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCAc                                                                    >  1:46452/1‑103 (MQ=255)
                                                                                                                                ctctGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTCCGGc                    <  2:556258/150‑1 (MQ=255)
                                                                                                                                   tGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCgtcgcgaaagggccgcagagcggccctgcgatcttgaatgaaccgcattaacgccaaggcgggattttcgtgcttaccacatcagatcatcagggatcttgt                 <  1:134742‑M2/150‑103 (MQ=255)
                                                                                                                                          cTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTCCGGCACAACGGCTa          >  1:381863/1‑150 (MQ=255)
                                                                                                                                           tttCGCGACgcaaggccagcgtcttctagagaacgcgcgtgtcgcgaaagggccgcagagcggccctgcgatcttgaatgaaccgcattaacgccaaggcgggattttcgtgcttaccacatcagatcatcagggatcttgtaagccgcg         >  1:406710‑M1/1‑9 (MQ=255)
                                                                                                                                           tttCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCgtcgcgaaagggccgcagagcggccctgcgatcttgaatgaaccgcattaacgccaaggcgggattttcgtgcttaccacatcagatcatcagggatcttgtaagccgcg         <  1:576193‑M2/150‑111 (MQ=255)
                                                                                                                                           tttCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTCCGGCACAACGGCTAg         >  2:601821/1‑150 (MQ=255)
                                                                                                                                               gcgACgcaaggccagcgtcttctagagaacgcgcgtgtcgcgaaagggccgcagagcggccctgcgatcttgaatgaaccgcattaacgccaaggcgggattttcgtgcttaccacatcagatcatcagggatcttgtaagccgcgtacg     <  2:463243‑M1/150‑146 (MQ=255)
                                                                                                                                               gcgACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCgtcgcgaaagggccgcagagcggccctgcgatcttgaatgaaccgcattaacgccaaggcgggattttcgtgcttaccacatcagatcatcagg                         >  2:445980‑M2/1‑36 (MQ=255)
                                                                                                                                               gcgACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCgtcgcgaaagggccgcagagcggccctgcgatcttgaatgaaccgcattaacgccaaggcgggattttcgtgcttaccacatcagatcatcagg                         >  2:445988‑M2/1‑36 (MQ=255)
                                                                                                                                               gcgACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCgtcgcgaaagggccgcagagcggccctgcgatcttgaatgaaccgcattaacgccaaggcgggattttcgtgcttaccacatcagatcatcagg                         <  1:445988‑M2/130‑95 (MQ=255)
                                                                                                                                               gcgACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCgtcgcgaaagggccgcagagcggccctgcgatcttgaatgaaccgcattaacgccaaggcgggattttcgtgcttaccacatcagatcatcagg                         <  1:445980‑M2/130‑95 (MQ=255)
                                                                                                                                                  acacGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTCCGGCACAACGGCTAGCCTCCAg  >  1:71066/1‑150 (MQ=255)
                                                                                                                                                  acacGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTCCGGCACAACGGCTAGCCTCCAg  >  1:71062/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
AACAGCCATGTCGAGGAAACGGACGAAGACGACCCGTACGCGGCTTACAAGATCCCTGATGATCTGATGTGGTAAGCACGAAAATCCCGCCTTGGCGTTAATGCGGTTCATTCAAGATCGCAGGGCCGCTCTGCGGCCCTTTCGCGACACGCGCGTTCTCTAGAAGACGCTGGCCTTGCGTCGCGAAAGGGTTCGTCGCTCACCGAACCGGCGCAGTTTGTCAGTTGCACTATCGTCTTGGTTCCCTCCGTTCGTCGGACATTTCTTAGGGTTCCGGCACAACGGCTAGCCTCCAG  >  NC_002947/1859595‑1859890

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