New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 1119638 =44 (1.200)9 (0.250) 9/264 2.6 17.4% coding (626/1494 nt) pepA cytosol aminopeptidase
?NC_002947 1119667 = 43 (1.210)coding (597/1494 nt) pepA cytosol aminopeptidase
Rejected: Frequency below/above cutoff threshold.

GGCCTTTTTGCCGCCCTGGTAATTGAGCACGATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCGGCCTGGCCGGCTTTGTCGGC  >  NC_002947/1119490‑1119783
                                                                                                                                                    |                                                                                                                                                  
gccgcgcgcgctgaaaaaggtcaccctgctggccgacaaagccggccaggccgaagtcgagcgtgcggtcaagcatgccagcgccatcgccagcggcatggccttcacccgcgacctcggcaacctgccgacctacctttgccacccaAg                                                                                                                                                   >  1:891149‑M1/149‑150 (MQ=255)
           cgccCTGGTAATTGAGCACGATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCt                                                                                                                                        >  2:592972/1‑150 (MQ=255)
           cgccCTGGTAATTGAGCACGATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCt                                                                                                                                        >  1:267657/1‑150 (MQ=255)
            gccCTGGTAATTGAGCACGATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGGTCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTg                                                                                                                                       >  2:105851/1‑150 (MQ=255)
              ccTGGTAATTGAGCACGATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGtt                                                                                                                                     >  1:899426/1‑150 (MQ=255)
                 ggTAATTGAGCACGATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGGCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCgg                                                                                                                                 >  1:277640/1‑150 (MQ=255)
                           cACGATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCt                                                                                                                       >  2:887994/1‑150 (MQ=255)
                           cACGATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCt                                                                                                                       >  2:124076/1‑150 (MQ=255)
                              gATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTgg                                                                                                                    <  1:242552/150‑1 (MQ=255)
                               gccgacaaagccggccaggccgaagtcgagcgtgcggtcaagcatgccagcgccatcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccacccaAGCTACCTGGCCTGTTCGGCCAGGt                                                                                                                           >  1:663150‑M1/118‑142 (MQ=255)
                               gccgacaaagccggccaggccgaagtcgagcgtgcggtcaagcatgccagcgccatcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccacccaAGCTACCTGGCCTGTTCGGCCAGGt                                                                                                                           <  2:663150‑M1/25‑1 (MQ=255)
                                tCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGt                                                                                                                  >  1:794796/1‑150 (MQ=255)
                                    gCGTGGCGGCTGGTCGCTTCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCa                                                                                                              >  1:421030/1‑150 (MQ=255)
                                     aaagccggccaggccgaagtcgagcgtgcggtcaagcatgccagcgccatcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccacccaAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTg                                                                                                                 <  1:414164‑M1/35‑1 (MQ=255)
                                     aaagccggccaggccgaagtcgagcgtgcggtcaagcatgccagcgccatcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccacccaAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTg                                                                                                                 >  2:414164‑M1/112‑146 (MQ=255)
                                         ccggccaggccgaagtcgagcgtgcggtcaagcatgccagcgccatcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccaccccAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTggg                                                                                                                   <  2:155516‑M1/33‑1 (MQ=255)
                                         ccggccaggccgaagtcgagcgtgcggtcaagcatgccagcgccatcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccaccccAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTggg                                                                                                                   >  1:155516‑M1/108‑140 (MQ=255)
                                         ccggccaggccgaagtcgagcgtgcggtcaagcatgccagcgccatcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccacccaAGCt                                                                                                                                                 <  2:567197‑M1/4‑1 (MQ=255)
                                         ccggccaggccgaagtcgagcgtgcggtcaagcatgccagcgccatcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccacccaAGCt                                                                                                                                                 >  1:567197‑M1/108‑111 (MQ=255)
                                             cTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTggg                                                                                                     <  1:124076/150‑1 (MQ=255)
                                                gTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGggcgg                                                                                                  >  1:95117/1‑150 (MQ=255)
                                                 tCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGggcggc                                                                                                 >  2:812664/1‑150 (MQ=255)
                                                        ccTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGc                                                                                          >  1:380438/1‑150 (MQ=255)
                                                         cTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCtt                                                                                                                      <  2:367136/121‑1 (MQ=255)
                                                         cTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCtt                                                                                                                      >  1:367136/1‑121 (MQ=255)
                                                            gCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAg                                                                                      <  2:614533/150‑1 (MQ=255)
                                                                aCGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTcg                                                                                  >  1:335438/1‑150 (MQ=255)
                                                                   gCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCgg                                                                               <  1:105851/150‑1 (MQ=255)
                                                                     gTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGt                                                                             >  1:852916/1‑150 (MQ=255)
                                                                         aaCGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAg                                                                         >  1:82673/1‑150 (MQ=255)
                                                                                     cccAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTg                                                             >  2:802123/1‑150 (MQ=255)
                                                                                      tcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccacccaAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGt                                                                                                                  <  2:595886‑M1/34‑1 (MQ=255)
                                                                                      tcgccagcggcatggccttcacccgcgacctcggcaacctgccgcccaacctttgccacccaAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGt                                                                                                                  >  1:595886‑M1/63‑96 (MQ=255)
                                                                                      ccAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTgg                                                            <  1:802123/150‑1 (MQ=255)
                                                                                      ccAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTgg                                                            >  1:266023/1‑150 (MQ=255)
                                                                                      ccAGGTCCTTAATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTgg                                                            >  1:446144/1‑150 (MQ=255)
                                                                                                 aTCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTg                                                 <  1:121383/150‑1 (MQ=255)
                                                                                                  tCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCACGGTGGGCAAAGGTTGGGCGGCAGGTTGCCGACGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTgg                                                >  1:384884/1‑150 (MQ=255)
                                                                                                    tttttCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCa                                              >  1:71387/1‑150 (MQ=255)
                                                                                                         ctcGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCtt                                         <  2:794796/150‑1 (MQ=255)
                                                                                                          tcGTCCAGCACTTCAACCTTGAGGGCATTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGggcggc                                                                                                 >  2:520977/1‑93 (MQ=255)
                                                                                                          tcGTCCAGCACTTCAACCTTGAGGGCATTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGggcggc                                                                                                 <  1:520977/93‑1 (MQ=255)
                                                                                                             tCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGAcc                                     >  2:737484/1‑150 (MQ=255)
                                                                                                              ccAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCg                                    >  2:349544/1‑150 (MQ=255)
                                                                                                                    cTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCAt                                                                    >  1:623589/1‑112 (MQ=255)
                                                                                                                    cTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCAt                                                                    <  2:623589/112‑1 (MQ=255)
                                                                                                                     ttCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTc                             >  2:660876/1‑150 (MQ=255)
                                                                                                                       cAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGGTGGCATGCTTGACCGCACGCTcta                           <  2:335438/150‑3 (MQ=255)
                                                                                                                            ttGAGGGCCTTCGGCGCCTTGCCCAGCTCCTGGGCCGGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCg                      <  1:592972/150‑1 (MQ=255)
                                                                                                                              gAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCGGc                    <  2:277640/150‑1 (MQ=255)
                                                                                                                                   ccTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCGGCCTGGc               >  1:548905/1‑150 (MQ=255)
                                                                                                                                        tgCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCGGCCTGGCCGGCt          <  2:899426/150‑1 (MQ=255)
                                                                                                                                         tttgccacccaAGCTACCTGGCCGAACAGGCCAAGGAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCGGCCTGGCCGGCtt         <  2:507952‑M1/139‑1 (MQ=255)
                                                                                                                                              ttGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCGGCCTGGCCGGCTTTGTCg    <  2:421030/150‑1 (MQ=255)
                                                                                                                                              ttGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCGGCCTGGCCGGCTTTGTCg    <  1:887994/150‑1 (MQ=255)
                                                                                                                                                gCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCGGCCTGGCCGGCTTTGTCGGc  <  1:737484/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GGCCTTTTTGCCGCCCTGGTAATTGAGCACGATCAGGCGTGGCGGCTGGTCGCTGCCCTGGCCCACGGCGTAGAACGCGCCCATGCCCAGGTCCTTGATCTTTTTCTCGTCCAGCACTTCAACCTTGAGGGCCTTGTGCGCCTTGCCCAGCTCCTTGGCCTGTTCGGCCAGGTAGCTTGGGTGGCAAAGGTTGGGCGGCAGGTTGCCGAGGTCGCGGGTGAAGGCCATGCCGCTGGCGATGGCGCTGGCATGCTTGACCGCACGCTCGACTTCGGCCTGGCCGGCTTTGTCGGC  >  NC_002947/1119490‑1119783

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.