New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_002947 = 559417023 (0.750)6 (0.250) 5/214 2.6 21.7% intergenic (‑66/+79) PP_4921/thiC NCS1 family transporter/phosphomethylpyrimidine synthase
?NC_002947 = 5594207 25 (1.020)intergenic (‑103/+42) PP_4921/thiC NCS1 family transporter/phosphomethylpyrimidine synthase

AAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTG  >  NC_002947/5594026‑5594314
                                                                                                                                                |                                                                                                                                                
aaaaCAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaag                                                                                                                                             >  1:53850‑M1/1‑145 (MQ=255)
           cgAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGt                                                                                                                                  <  2:305760/150‑1 (MQ=255)
                  cGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGACGGTGCCATGc                                                                                                                           >  1:407458/1‑150 (MQ=255)
                              gctgGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgctt                                                                                                               <  2:168445/150‑1 (MQ=255)
                                tgGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCAAGCACTGCGGtgcatgg                                                                                                             <  1:79508/150‑2 (MQ=38)
                                      gggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcg                                                                                                       >  2:678945‑M2/1‑144 (MQ=255)
                                       ggATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgg                                                                                                      <  2:53850‑M1/150‑45 (MQ=255)
                                        gATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcggg                                                                                                     <  1:727542‑M2/150‑9 (MQ=255)
                                                ggTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccg                                                                                             >  2:259716‑M2/1‑134 (MQ=255)
                                                     ggAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctccc                                                                                        <  1:562174‑M2/150‑22 (MQ=255)
                                                     ggAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTACCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgg                                                                                                      <  1:465702‑M1/136‑45 (MQ=255)
                                                     ggAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTACCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgg                                                                                                      >  2:465702‑M1/1‑92 (MQ=255)
                                                                gctggGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagg                                                                                   <  2:131605‑M2/144‑27 (MQ=255)
                                                                gctggGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagg                                                                                   >  1:131605‑M2/1‑118 (MQ=255)
                                                                gctggGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAgg                                                                             >  1:471474/1‑150 (MQ=255)
                                                                     gTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgctt                                                                                                               >  2:329657/1‑111 (MQ=255)
                                                                     gTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGtgctt                                                                                                               <  1:329657/111‑1 (MQ=255)
                                                                           tcatggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaagcccg                                                                                             <  2:102127‑M1/123‑54 (MQ=255)
                                                                           tcatggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaagcccg                                                                                             >  1:102127‑M1/1‑70 (MQ=255)
                                                                             atggtcatgTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAg                                                                <  2:727542/150‑1 (MQ=255)
                                                                                         gCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGtgctgc                                                                       >  1:633071/1‑131 (MQ=255)
                                                                                         gCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGtgctgc                                                                       <  2:633071/131‑1 (MQ=255)
                                                                                          cAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGAGGACAACCGCGCACCCACAGAGTAAGGTGCTGCGGTATAGACTTTGTGGTACa                                                   >  2:449236/1‑150 (MQ=37)
                                                                                                        aGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTcc                                     <  1:493170/150‑1 (MQ=255)
                                                                                                          tGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTg                                   <  2:562174/150‑1 (MQ=255)
                                                                                                              aTACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGtt                                                                   <  2:386109/114‑1 (MQ=255)
                                                                                                              aTACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGtt                                                                   >  1:386109/1‑114 (MQ=255)
                                                                                                                     accCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGAGTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCt                        >  1:3979/1‑150 (MQ=255)
                                                                                                                       ccTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagc                            >  1:131605‑M1/1‑26 (MQ=255)
                                                                                                                       ccTGTGGGAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagc                            <  2:131605‑M1/144‑119 (MQ=255)
                                                                                                                       ccTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCg                      <  1:747347/150‑1 (MQ=255)
                                                                                                                       ccTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGGTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCg                      >  2:717498/1‑150 (MQ=255)
                                                                                                                            gggAGCGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcc                 <  1:562174‑M1/150‑130 (MQ=255)
                                                                                                                              gAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCt               >  1:41433/1‑150 (MQ=255)
                                                                                                                                 cGGGCTTGCCCGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgacc            >  2:259716‑M1/1‑16 (MQ=255)
                                                                                                                                 cGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatga                                       <  2:102127‑M2/123‑71 (MQ=255)
                                                                                                                                 cGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatga                                       >  1:102127‑M2/1‑53 (MQ=255)
                                                                                                                                    gCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCa         >  1:342433/1‑150 (MQ=255)
                                                                                                                                      ttGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATg       >  2:125077/1‑150 (MQ=255)
                                                                                                                                        gCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGcc     <  2:471474/150‑1 (MQ=255)
                                                                                                                                         cccGCGAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatc    <  1:727542‑M1/150‑143 (MQ=255)
                                                                                                                                          ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggtaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcc                 <  1:465702‑M2/136‑93 (MQ=255)
                                                                                                                                          ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggtaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcc                 >  2:465702‑M2/1‑44 (MQ=255)
                                                                                                                                          ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatcc   <  2:53850‑M2/150‑107 (MQ=255)
                                                                                                                                          ccGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGGGGCGGTTTAGACTTTGTGGTACAGCTGACTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCtt   >  1:239345/1‑150 (MQ=255)
                                                                                                                                           cgcgAAgcaagcactgcggtgcatggcaccggctttgccggtgttcgcgggcaagcccgctcccacagggtggtgtatgccacttaacatcactcctgcggacatgaccatgacatcacccagccagttctctcccgaccatccaatccc  >  2:678945‑M1/1‑6 (MQ=255)
                                                                                                                                                |                                                                                                                                                
AAAACAGGTCGCGAGCACCGAAAATGCGCTGCTGGATGGGGATTGGATGGTCGGGAGAGAACTGGCTGGGTGATGTCATGGTCATGTCCGCAGGAGTGATGTTAAGTGGCATACACCACCCTGTGGGAGCGGGCTTGCCCGCGAACACCGGCAAAGCCGGTGCCATGCACCGCAGTGCTTGCTTCGCGGGCAAGCCCGCTCCCACAGGGTAAGGTGCTGCGGTTTAGACTTTGTGGTACAGCTGGCTGCCTTCCTGGCGGAACCGCTCGGCCTGCTCACGCATGCCTTG  >  NC_002947/5594026‑5594314

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.