New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 480308245 (1.510)6 (0.210) 6/266 2.6 12.8% coding (277/972 nt) PP_4231 xanthine dehydrogenase accessory factor
?NC_002947 = 4803104 39 (1.370)coding (255/972 nt) PP_4231 xanthine dehydrogenase accessory factor
Rejected: Frequency below/above cutoff threshold.

TCACCTCACGCAGTACCCGGCGCTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAA  >  NC_002947/4802959‑4803242
                                                                                                                                                 |                                                                                                                                           
tcacctcacGCAGTACCCGGCGCTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGctgc                                                                                                                                         <  2:102350‑M2/150‑5 (MQ=255)
 cacctcacGCAGTACCCGGCGCTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACCCCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGctgcc                                                                                                                                        >  2:249472‑M2/1‑145 (MQ=255)
  acctcacGCAGTACCCGGCGCTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACAGCAGGCGAACGCGCGAACGCGGGTCGctgccg                                                                                                                                       <  1:426160‑M2/150‑7 (MQ=255)
   cctcacGCAGTACCCGGCGCTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGctgccgt                                                                                                                                      <  1:300310‑M2/150‑8 (MQ=255)
     tcacGCAGTACCCGGCGCTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCg                                                                                                                                    >  1:670017/1‑150 (MQ=255)
           aGTACCCGGCGCTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATaa                                                                                                                              <  1:313195/150‑1 (MQ=255)
             tACCCGGCGCTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACg                                                                                                                            >  1:261811/1‑150 (MQ=255)
                    cgcTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATg                                                                                                                     <  2:480871/150‑1 (MQ=255)
                           ccAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCg                                                                                                             >  1:201429/1‑150 (MQ=255)
                                  cgcACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCgg                                                                                                      <  1:353891/150‑1 (MQ=255)
                                   gcACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCggg                                                                                                     >  1:421182/1‑150 (MQ=255)
                                     aCGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGaa                                                                                                   <  1:102350/150‑1 (MQ=255)
                                          ccAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTccc                                                                                              >  2:218835/1‑150 (MQ=255)
                                            aaCTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCg                                                                                            >  2:536769/1‑150 (MQ=255)
                                            aaCTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCg                                                                                            >  2:266949/1‑150 (MQ=255)
                                             aCTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGa                                                                                           <  2:300656/150‑1 (MQ=255)
                                              cTCCCGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAg                                                                                          <  2:613265/150‑1 (MQ=255)
                                                         gcgcCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGtt                                                                               >  2:639931/1‑150 (MQ=255)
                                                           gcCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCACCCGGCTCCGGGCATTCCCCGAGGGCCACGCGTTCg                                                                             >  2:290163/1‑150 (MQ=255)
                                                            ccTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACAGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCga                                                                            <  2:426160/150‑1 (MQ=255)
                                                                cACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGaaa                                                                        <  1:53647/150‑1 (MQ=255)
                                                                   cTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAgtc                                                                     >  1:641034/1‑150 (MQ=255)
                                                                         ttCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTc                                                               <  2:670017/150‑1 (MQ=255)
                                                                            aaacGGCAGGTTATCCACCAGCACTAACAGCACACAGCGATACGGCAGGCGAGGGCGCTAACGCGGGTCGCTGACGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGGGTTCGAGAAAGTCGTCTTCGac                                                            >  2:118199/4‑150 (MQ=255)
                                                                            aaCCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGac                                                            >  2:439239/1‑150 (MQ=255)
                                                                               cGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGAcaca                                                         <  2:261811/150‑1 (MQ=255)
                                                                               cGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGAcaca                                                         <  2:300310/150‑1 (MQ=255)
                                                                                      ttCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCg                                                  >  2:470091/1‑150 (MQ=255)
                                                                                       tctcCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCgg                                                 >  1:468840/1‑150 (MQ=255)
                                                                                            accaccACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGc                                            <  1:10776/150‑1 (MQ=255)
                                                                                            accaGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGc                                            >  1:684387/1‑150 (MQ=255)
                                                                                            accaGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGc                                            <  1:470091/150‑1 (MQ=255)
                                                                                                    ttCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAg                                    >  1:343910/1‑150 (MQ=255)
                                                                                                           acacCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTgg                             <  1:639931/150‑1 (MQ=255)
                                                                                                               cGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGc                         <  1:218835/150‑1 (MQ=255)
                                                                                                                  cacaCGGCAGGCGAACGCGCCAACGCGGGGCGCTGCCGTAGCCATAACGCACGATGGCAACCGGCCCGGGGAATTCCCCGAGCGCCACGCGTTCGCGAAAGTCGCCTTCGACACAGCCACCCGACAGCGAGCCCCGCCACTGGCCGCtcg                      >  2:108604/1‑148 (MQ=255)
                                                                                                                   acacGGCAGGCGAACGCGCGAACGCGGGTCGctgccgtgtggcggtgtgctggaagtgctggtggagaacctgccggttgaatgtgaggtgcaggcgcatctgcgtgagttgtaacgtgc                                                    <  2:406170‑M2/120‑90 (MQ=255)
                                                                                                                   acacGGCAGGCGAACGCGCGAACGCGGGTCGctgccgtgtggcggtgtgctggaagtgctggtggagaacctgccggttgaatgtgaggtgcaggcgcatctgcgtgagttgtaacgtgc                                                    >  1:406170‑M2/1‑31 (MQ=255)
                                                                                                                    cacGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCg                    <  2:468840/150‑1 (MQ=255)
                                                                                                                     acGGCAGGCGAACGCGCGAACGCGGGTCGctgccgtgtggcggtgtgctggaagtgctggtggagaacctgccggttgaatgtgag                                                                                    >  2:322247‑M2/1‑29 (MQ=255)
                                                                                                                     acGGCAGGCGAACGCGCGAACGCGGGTCGctgccgtgtggcggtgtgctggaagtgctggtggagaacctgccggttgaatgtgag                                                                                    <  1:322247‑M2/86‑58 (MQ=255)
                                                                                                                     acGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGACACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGt                   >  2:681910/1‑150 (MQ=255)
                                                                                                                         cAGGCGAACGCGCGAACGCGGGTCGctgccgtgtggcggtgtgctggaagtgctggtggagaacctgccggttgaatgtgaggtgcaggcgcatctgcgtgagttggaacgtgcgctgctgggccagcgccgggtactgcgtgaggtgag                <  1:630289‑M2/150‑126 (MQ=255)
                                                                                                                          aGGCGAACGCGCGAACGCGGGTCGctgccgtgtggcggtgtgctggaagtgctggtggagaacctgccggttgaatgtgaggtgcaggcgcatctgcgtgagttggaacgtgcgctgctgggccagcgccgggtactgcgtgaggtgagc               <  1:213557‑M2/150‑127 (MQ=255)
                                                                                                                              gaacgcgcgaacgcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCca          <  2:343910/150‑1 (MQ=255)
                                                                                                                               ctcgcgcgaacgcgGGTCGCTGCAGTCGTAATCACGCTTGAGGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGGTGAACGCcag         <  2:695014/148‑1 (MQ=255)
                                                                                                                                  gcgcgaacgcgGGTCGCTGCCGTCGCCATACAGAGGAAGAACAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGGTGAACGCcagcag      >  1:695014/1‑150 (MQ=255)
                                                                                                                                     cgaacgcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCTGCTGGCGGTGAACGCcagcggcca   <  2:421182/150‑6 (MQ=255)
                                                                                                                                      gaacgcgGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGaa  <  1:439239/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                           
TCACCTCACGCAGTACCCGGCGCTGGCCCAGCAGCGCACGTTCCAACTCACGCAGATGCGCCTGCACCTCACATTCAACCGGCAGGTTCTCCACCAGCACTTCCAGCACACCGCCACACGGCAGGCGAACGCGCGAACGCGGGTCGCTGCCGTCGCCATAACGCACGATGGCAACCGGCTCCGGGAATTCCCCGAGGGCCACGCGTTCGAGAAAGTCGTCTTCGACACAGCCACCGGACAGCGAGCCCAGCCACTGGCCGCTGGCGTTCACCGCCAGCAGCGAA  >  NC_002947/4802959‑4803242

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.