New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 583970453 (0.940)3 (0.060) 3/248 NT 6.0% intergenic (+81/+22) yheT/sseA putative hydrolase/3‑mercaptopyruvate sulfurtransferase
?NC_002947 = 5839733 44 (0.810)coding (848/855 nt) sseA 3‑mercaptopyruvate sulfurtransferase
Rejected: Coverage evenness skew score above cutoff.

GCCTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAGGTTATGGCAAGCGGTTACCC  >  NC_002947/5839559‑5839850
                                                                                                                                                 |                                                                                                                                                  
gCCTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCgg                                                                                                                                                 >  6:287864/1‑149 (MQ=255)
  cTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCgggg                                                                                                                                               <  3:232399/149‑1 (MQ=255)
  cTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCgggg                                                                                                                                               <  5:109569/149‑1 (MQ=255)
     cGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCggggg                                                                                                                                              <  2:255152/147‑1 (MQ=255)
     cGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGt                                                                                                                                             <  7:19231/148‑1 (MQ=255)
       aaaTCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTc                                                                                                                                           >  5:246046/1‑148 (MQ=255)
        aaTCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCaa                                                                                                                                         >  7:41196/1‑149 (MQ=255)
              gggTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGt                                                                                                                                   <  6:246046/149‑1 (MQ=255)
                     actTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTgt                                                                                                                               <  4:229029/146‑2 (MQ=255)
                     actTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTgt                                                                                                                               >  3:229029/1‑145 (MQ=255)
                      ctTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTa                                                                                                                            >  2:251138/1‑148 (MQ=255)
                      ctTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTAc                                                                                                                           <  3:237940/149‑1 (MQ=255)
                      ctTGGAACGACGGATTCCGCAGTGGCGGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGTACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCtt                                                                                                                             <  3:276208/147‑1 (MQ=255)
                           aaCGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGcc                                                                                                                      >  1:14276/1‑149 (MQ=255)
                              gacgGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGc                                                                                                                               >  4:19023/1‑137 (MQ=255)
                              gacgGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGc                                                                                                                               <  3:19023/137‑1 (MQ=255)
                               acgTATTCCGGAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTTCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCTCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTg                                                                                                                  <  7:178128/149‑1 (MQ=255)
                                              gCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTg                                                                                                                  <  1:60661/134‑1 (MQ=255)
                                              gCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTg                                                                                                                  >  2:60661/1‑134 (MQ=255)
                                                tggtggACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCgg                                                                                                 <  1:155556/149‑1 (MQ=255)
                                                  gtggACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACt                                                                                                                          <  2:274186/122‑1 (MQ=255)
                                                  gtggACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACt                                                                                                                          >  1:274186/1‑122 (MQ=255)
                                                            gctagcCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCca                                                                                     >  2:169943/1‑149 (MQ=255)
                                                                gcCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCcactca                                                                                 <  5:287864/149‑1 (MQ=255)
                                                                  ccGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCcactcact                                                                               >  1:59675/1‑149 (MQ=255)
                                                                    gcgAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTcc                                                                             <  8:41196/149‑1 (MQ=255)
                                                                        aTTGGTAATGGGTGCGCGCCGCCTGTGGCTGCGGGCATGCCGGCGAACACGGGCTTTGCCGGTGCCACCCACCGGGGGGGTTCAACCGGGTATGCTTAGTCGCTGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCAc                                                                           <  1:251138/147‑1 (MQ=255)
                                                                             tAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCcac                                                                                    <  8:234267/133‑1 (MQ=255)
                                                                             tAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCg                                                                    >  5:151651/1‑149 (MQ=255)
                                                                             tAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACACACCGCGGGGGTTCAAACCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCcac                                                                                    >  7:234267/1‑133 (MQ=255)
                                                                              aaaGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCg                                                                    >  5:100163/1‑148 (MQ=255)
                                                                                          cccccTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCGCTCACTCCACGACCCCGCATACAGCTTGcc                                                       >  3:298752/1‑149 (MQ=255)
                                                                                          cccccTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGcc                                                       >  6:74757/1‑149 (MQ=255)
                                                                                             ccTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAg                                                    >  7:219855/1‑149 (MQ=255)
                                                                                              cTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGc                                                   <  6:151651/149‑1 (MQ=255)
                                                                                                 tgGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCgg                                                 <  4:232721/148‑1 (MQ=255)
                                                                                                 tgGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCg                                                  <  8:219855/147‑1 (MQ=255)
                                                                                                    gAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTAc                                             >  4:324638/1‑149 (MQ=255)
                                                                                                    gAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCAACGACCCCGCATACAGCTTGCCCAGCGGGTAc                                             >  3:31734/1‑149 (MQ=255)
                                                                                                          gCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgcc                                       >  4:13523/1‑149 (MQ=255)
                                                                                                            aTGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgccag                                     <  2:59675/149‑1 (MQ=255)
                                                                                                                 cgcgAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgccagcgc                                  >  1:346471/1‑147 (MQ=255)
                                                                                                                 cgcgAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgccagcgc                                  >  6:7424/1‑147 (MQ=255)
                                                                                                                  gcgAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTACACGACCCCGCATACAGCTTGCCCAGCGGGTACCCcgccagcgccag                               >  4:272641/1‑149 (MQ=255)
                                                                                                                          cGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCAcc                                                                                                            <  6:345266/64‑1 (MQ=255)
                                                                                                                          cGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCaaa                         >  4:89254/1‑147 (MQ=255)
                                                                                                                            gCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAgg                     >  4:309964/1‑149 (MQ=255)
                                                                                                                               ttGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAGGTTa                  <  5:7424/149‑1 (MQ=255)
                                                                                                                                 gCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAGGTTATg                >  1:168458/1‑149 (MQ=255)
                                                                                                                                  ccGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAGGTTATg                <  2:346471/148‑1 (MQ=255)
                                                                                                                                  ccGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAGGTTATg                <  2:14276/148‑1 (MQ=255)
                                                                                                                                   ttaacctacccaccGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAGGTTATGGc              <  1:177624/142‑1 (MQ=255)
                                                                                                                                     gTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAGGTTATGGCaa            <  5:74757/149‑1 (MQ=255)
                                                                                                                                       gCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCACATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAGGTTATGGCAAGc          >  8:82302/1‑149 (MQ=255)
                                                                                                                                          acccaccGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCg                                                                    <  6:235782/88‑1 (MQ=255)
                                                                                                                                          acccaccGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCg                                                                    >  5:235782/1‑88 (MQ=255)
                                                                                                                                               ccGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCATCGCAAACAGGTTATGGCAAGCGGATAccc  <  4:31734/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCCTGCGAAATCCGGGGTACTACTTGGAACGACGGATTCCGCAGTGGCTGGTGGACGGTAGCTAGCCCGCGAATTGGTAAAGGCTGCGCGCCCCCTGTGGGAGCGGGCATGCCCGCGAACACCGGCATTGCCGGTGCCACCCACCGCGGGGGTTCAACCCGGTCTGCTTACTCGCCGGTGGCCACCGCATGCTGCGGGTTGTTGATCCACTCACTCCACGACCCCGCATACAGCTTGCCCAGCGGGTACCCCGCCAGCGCCAGCGCAAACAGGTTATGGCAAGCGGTTACCC  >  NC_002947/5839559‑5839850

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