New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 559411926 (0.500)3 (0.070) 3/204 NT 12.2% intergenic (‑66/+79) PP_4921/thiC NCS1 family transporter/phosphomethylpyrimidine synthase
?NC_002947 = 5594156 23 (0.570)intergenic (‑103/+42) PP_4921/thiC NCS1 family transporter/phosphomethylpyrimidine synthase

GGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCGTAC  >  NC_002947/5594008‑5594303
                                                                                                                                                    |                                                                                                                                                   
ggATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgcttgc                                                                                                                                                     <  7:205208/149‑1 (MQ=255)
 gATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGGGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgcttgc                                                                                                                                                     <  1:221364/148‑1 (MQ=255)
    cctgactaatggcacaatgatttctacgatttggttataactccaatcgcagtgaggaataaatttggggccctaatgccttcctgtgggagcgggcacGCCCGCGAACACCGGCAAAGCCGGTGGCATCCACCGCAGTGGATGCTTc                                                                                                                                                  <  1:304541‑M2/49‑1 (MQ=255)
       ctgactaatggcacaatgatttctacgatttggttataactccaatcgcagtgaggaataaatttggggccccgcgccccttgtgggagcgggcacGCCCGCGAACACCGGCAAAGCCGGTGGCATCCACCGCAGTGGATGCTTCGCgg                                                                                                                                              <  5:77697‑M2/53‑1 (MQ=255)
        actaatggcacaatgatttctacgatttggttataactccaatcgcagtgaggaataaatttggggccctaatgccttcctgtgggagcgggcacGCCCGCGAACACCGGCAAAGCCGGTGGCATCCACCGCAGTGGATGCTTCGCgg                                                                                                                                              <  6:254837‑M2/53‑1 (MQ=255)
                     gagagaACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctccc                                                                                                                                <  2:17769‑M2/149‑22 (MQ=255)
                                          tcatggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGtt                                                                                                           <  5:261245/149‑1 (MQ=255)
                                                tcatgTCCGCAGGAGTGATGTTAAGTGGCATACTCCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACt                                                                                                     >  1:44183/1‑149 (MQ=255)
                                                      ccGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCCGTTTAGACTTtgttg                                                                                               >  8:178563/1‑147 (MQ=255)
                                                       cGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGt                                                                                              <  2:44183/149‑1 (MQ=255)
                                                                gATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTTGCCGCGCACACCGGCAAAGCCGGGGCCCATCCCCCCAGTGCTTGGTTGGCGGGGAAGCCCCCGCCCCCCGGGGTAGGTGGGGGGGGTTTGTCTTTTTTGTACAgctggg                                                                                     >  8:288495/1‑148 (MQ=12)
                                                                                   caccCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTa                                                                                                         <  1:282784/110‑1 (MQ=255)
                                                                                                 gggCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCAc                                                    >  3:298947/1‑149 (MQ=255)
                                                                                                   gCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGc                                                  <  7:288495/149‑1 (MQ=255)
                                                                                                     ttGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgac                                                                                    >  6:63638‑M2/1‑48 (MQ=255)
                                                                                                     ttGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgac                                                                                    <  5:63638‑M2/113‑66 (MQ=255)
                                                                                                       gCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGc                                              >  2:26162/1‑149 (MQ=255)
                                                                                                        cccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGcc                                             >  2:327075/1‑149 (MQ=255)
                                                                                                        cccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGcc                                             >  7:11412/1‑149 (MQ=255)
                                                                                                        cccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGcc                                             >  1:227516/1‑149 (MQ=255)
                                                                                                         ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatc                                            <  1:52045‑M2/149‑106 (MQ=255)
                                                                                                             tggcacCGGCTTTGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGc                                                                                  <  7:97750/104‑1 (MQ=37)
                                                                                                             tggcacCGGCTTTGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGc                                                                                  >  8:97750/4‑107 (MQ=37)
                                                                                                             gAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCt                                        >  5:108862/1‑149 (MQ=255)
                                                                                                                  ccGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCt                                                                                >  4:119954/1‑104 (MQ=255)
                                                                                                                  ccGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCt                                                                                <  3:119954/104‑1 (MQ=255)
                                                                                                                   cGGCAAAGCCGGTGCCATGAACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCg                                  >  4:156901/1‑149 (MQ=255)
                                                                                                                              gaagcaagCACTGCGGTGCATGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACg                       >  1:83956/8‑149 (MQ=255)
                                                                                                                               tGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACgg                      <  1:327075/149‑1 (MQ=255)
                                                                                                                                   aTGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTc                  <  1:26162/149‑1 (MQ=255)
                                                                                                                                          gCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCAAATCCACGGCCTCGATCTTg           >  8:199768/1‑149 (MQ=255)
                                                                                                                                           cggTGCATGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCtt            <  2:83956/145‑1 (MQ=255)
                                                                                                                                                   gcttCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCGTAc  >  5:241094/1‑149 (MQ=255)
                                                                                                                                                   gcttCGCGGGCAAGCCCGCTCCCACAGGGTAAGATGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCGTAc  >  5:103621/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTGCTCGACCGTCACATCCACGGCCTCGATCTTGGCAGCGTAC  >  NC_002947/5594008‑5594303

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