New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 606302833 (1.040)3 (0.100) 3/238 NT 8.7% intergenic (‑85/‑140) rubA/ubiC rubredoxin‑1/chorismate pyruvate‑lyase
?NC_002947 = 6063060 31 (1.020)intergenic (‑117/‑108) rubA/ubiC rubredoxin‑1/chorismate pyruvate‑lyase

CCTTCGGCTTCGTCGTAGATCAGGCCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACG  >  NC_002947/6062891‑6063177
                                                                                                                                         |                                                                                                                                                     
ccTTCGGCTTCGTCGTAGATCAGGCCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCggg                                                                                                                                           <  1:522319/150‑1 (MQ=255)
  ttCGGCTTCGTCGTAGATCAGGCCACAGACAATACATTGCCCCTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCt                                                                                                                                         >  1:96729/1‑150 (MQ=255)
  ttCGGCTTCGTCGTAGATCAGGCCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCt                                                                                                                                         >  1:38944/1‑150 (MQ=255)
   tCGGCTTCGTCGTAGATCAGGCCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCtt                                                                                                                                        <  1:736746/150‑1 (MQ=255)
                    cAGGCCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCa                                                                                                                       <  2:376919/150‑1 (MQ=255)
                     aGGCCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAg                                                                                                                      <  2:592688‑M2/150‑2 (MQ=255)
                       gCCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcg                                                                                                                    <  2:596817‑M2/150‑4 (MQ=255)
                                   actttgtcACTGGTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCtggccctgtgcgcgggcaagcccgccaggcctgcgcggtcactgtag                                                                                                        <  2:607088‑M1/147‑48 (MQ=255)
                                   aCATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTgg                                                                                                        <  1:342698/150‑1 (MQ=255)
                                         ccACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCtt                                                                                                  <  2:565046/150‑1 (MQ=255)
                                              ttttCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAg                                                                                             <  1:401326/150‑1 (MQ=255)
                                                ttCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCg                                                                                           <  2:264150/150‑1 (MQ=255)
                                                 tcatcaGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCgg                                                                                          <  1:245031/150‑1 (MQ=255)
                                                    ttaGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGccc                                                                                       >  1:28325/3‑150 (MQ=255)
                                                             tttCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCtggccctgtgcgcgggcaagcccgccaggcctgcgc                                                                                                                   >  1:338668‑M1/1‑77 (MQ=255)
                                                             tttCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCtggccctgtgcgcgggcaagcccgccaggcctgcgc                                                                                                                   <  2:338668‑M1/113‑37 (MQ=255)
                                                             tttCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgc                                                                                                                   <  1:302816‑M2/113‑5 (MQ=255)
                                                             tttCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgc                                                                                                                   >  2:302816‑M2/1‑109 (MQ=255)
                                                                  tCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGc                                                                         >  1:227261/1‑150 (MQ=255)
                                                                    ggTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCaa                                                                       <  2:28325/150‑1 (MQ=255)
                                                                          ggCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGtt                                                                 <  2:830709/150‑1 (MQ=255)
                                                                            cTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTgg                                                               <  2:227261/150‑1 (MQ=255)
                                                                             tGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTggg                                                              <  2:773206/150‑1 (MQ=255)
                                                                              ggggCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTgggg                                                             <  1:592688/150‑1 (MQ=255)
                                                                                 gCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGc                                                                                         <  1:558381/119‑1 (MQ=255)
                                                                                 gCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGc                                                                                         >  2:558381/1‑119 (MQ=255)
                                                                                  cTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAAt                                                         <  1:596817/150‑1 (MQ=255)
                                                                                       gtTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGACTTGCCCACCCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGc                                                    >  2:418889/1‑150 (MQ=255)
                                                                                          gctgctGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTa                                                                            >  2:657237/1‑123 (MQ=255)
                                                                                          gctgctGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGACTGGCGGGCTTGCCCGCGCACATGGCCAGCGCTCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTa                                                                            <  1:657237/123‑1 (MQ=255)
                                                                                            tgctgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGc                                               >  1:838719/1‑150 (MQ=255)
                                                                                               tgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTgggg                                                             >  2:629135/1‑133 (MQ=255)
                                                                                               tgATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTgggg                                                             <  1:629135/133‑1 (MQ=255)
                                                                                                    gCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCAGGCTAAGCTCgccgcc                                       >  2:594906/1‑150 (MQ=255)
                                                                                                      ctctTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCctc                                     >  2:568052/1‑150 (MQ=255)
                                                                                                                  aaGCCGCTCCTACAGTGACCGCGCAGGCCTGGCTGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAAc                         >  2:508227/1‑150 (MQ=255)
                                                                                                                           cTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgcggtcactgtaggagcggcttgcccgcgaagaggccatcagcagcaacacaagccccagcctgaaccgaggaaatagacctgatgaaccagtgacaaagt                <  2:607088‑M2/150‑104 (MQ=255)
                                                                                                                            tACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATAcc               <  1:568052/150‑1 (MQ=255)
                                                                                                                             aCAGTGACCGCGCAGGCAACGCGCGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCg              >  2:438411/1‑150 (MQ=255)
                                                                                                                                    ccGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTcgccgc                                        >  2:443364/1‑117 (MQ=255)
                                                                                                                                    ccGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTcgccgc                                        <  1:443364/117‑1 (MQ=255)
                                                                                                                                      gcgctggccctgtgcgcgggcaagcccgccaggcctgcgcggtcactgtaggagcggcttgcccgcgaagaggccatcagcagcaacacaagccccagcctgaaccgaggaaa                                          <  1:302816‑M1/113‑110 (MQ=255)
                                                                                                                                      gcgctggccctgtgcgcgggcaagcccgccaggcctgcgcggtcactgtaggagcggcttgcccgcgaagaggccatcagcagcaacacaagccccagcctgaaccgaggaaa                                          >  2:302816‑M1/1‑4 (MQ=255)
                                                                                                                                      gcgcAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgcggtcactgtaggagcggcttgcccgcgaagaggccatcagcagcaacacaagccccagcctgaaccgaggaaa                                          >  1:338668‑M2/1‑36 (MQ=255)
                                                                                                                                      gcgcAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAgcgcggtcactgtaggagcggcttgcccgcgaagaggccatcagcagcaacacaagccccagcctgaaccgaggaaa                                          <  2:338668‑M2/113‑78 (MQ=255)
                                                                                                                                       cgctggccctgtgcgcgggcaagcccgccaggcctgcgcggtcactgtaggagcggcttgcccgcgaagaggccatcagcagcaacacaagccccagcctgaaccgaggaaatagacctgatgaaaaagtggcaatgtattgtctgtggc    <  2:596817‑M1/150‑148 (MQ=255)
                                                                                                                                         ctggccctgtgcgcgggcaagcccgccaggcctgcgcggtcactgtaggagcggcttgcccgcgaagaggccatcagcagcaacacaagccccagcctgaaccgaggaaatagacctgatgaaaaagtggcaatgtattgtctgtggcct  <  2:592688‑M1/150‑150 (MQ=255)
                                                                                                                                         |                                                                                                                                                     
CCTTCGGCTTCGTCGTAGATCAGGCCACAGACAATACATTGCCACTTTTTCATCAGGTCTATTTCCTCGGTTCAGGCTGGGGCTTGTGTTGCTGCTGATGGCCTCTTCGCGGGCAAGCCGCTCCTACAGTGACCGCGCAGGCCTGGCGGGCTTGCCCGCGCACAGGGCCAGCGCGCAATGCGTGGGGGCTTTGTAGCGGCCCTGCCGCCAGTATGCAAGCAGTTGGGGCAATCATGCTAAGCTCGCCGCCTCTCTGGCTGTAACAAAGCATACCGACGTGTCGTACG  >  NC_002947/6062891‑6063177

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