New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? CP000730 1236098 =70 (1.010)3 (0.050) 3/282 2.2 4.9% coding (254/642 nt) yloS thiamine diphosphokinase
?CP000730 1236143 = 50 (0.770)coding (299/642 nt) yloS thiamine diphosphokinase
Rejected: Frequency below/above cutoff threshold.

TTCAGTCAGTAAAGAAGAACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCCTAAAGGTCAACATACAGTTGAAAAAGATAAG  >  CP000730/1235994‑1236288
                                                                                                                                                     |                                                                                                                                                  
aattttattttgttgatcgataacttctatatgaacatcatgtttataatatgcttttttcaataataactgaattgccccaaaaaagtgatctaatcgtccgcCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGc                                                                                                                                                   >  1:788395‑M1/105‑151 (MQ=255)
  cagtcaatAAATAAAAACGCCAACTTCTAACAGAACAGTTATAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTATAATTTATGGTGCAACAGGCgg                                                                                                                                                 <  2:667675/151‑1 (MQ=255)
   agtcagtAAAGAAGAACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGa                                                                                                                                                <  1:264773/151‑1 (MQ=255)
     tattttgttgatcgataacttctatatgaacatcatgtttataatatgcttttttcaataataactgaattgccccaaaaaagtgatctaatcgtccgcCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACg                                                                                                                                              >  2:434594‑M1/100‑151 (MQ=255)
        gtAAAGAAGAACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGAtt                                                                                                                                           >  1:1048302/1‑151 (MQ=255)
         tAAAGAAGAACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACACTTTATGGTGCAACAGGCGGACGATTa                                                                                                                                          >  2:70703/1‑151 (MQ=255)
         tAAAGAAGAACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTa                                                                                                                                          <  2:648461/151‑1 (MQ=255)
          aaaGAAGAACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAg                                                                                                                                         >  2:336298/1‑151 (MQ=255)
          aaaGAAGAACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAg                                                                                                                                         >  1:599910/1‑151 (MQ=255)
           aagaagAACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGa                                                                                                                                        >  2:910434/1‑151 (MQ=255)
            agaagaACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGAt                                                                                                                                       <  1:15653/151‑1 (MQ=255)
               agagCGCTAACTTCTAACAGATCATTTACAAATCAAACCAGTTCTGGCTGAATAAGCTGATACGGATTTAGCTTTAGTGGTTGATAAAGCTGTTGCACTTGGATGTGATAGTATTACGATTTATGGTGCAACAGGCGGACGATTAGATCAc                                                                                                                                    <  2:277811/151‑1 (MQ=255)
               agaaCGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGGCGATTAGCTCAc                                                                                                                                    >  1:459222/1‑151 (MQ=255)
               agaaCGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCAc                                                                                                                                    >  1:927463/1‑151 (MQ=255)
                gaaCGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACt                                                                                                                                   >  2:486126/1‑151 (MQ=255)
                         cTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTgggg                                                                                                                          >  1:557281/1‑151 (MQ=255)
                          ttCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGc                                                                                                                         <  1:413329/151‑1 (MQ=255)
                          ttCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGc                                                                                                                         >  1:570679/1‑151 (MQ=255)
                           tCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCa                                                                                                                        >  2:727491/1‑151 (MQ=255)
                              aacagaacagTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAAtt                                                                                                                     >  1:804422/1‑151 (MQ=255)
                               acagaacagTTACAAATCAAACCAGTTCAAGCTGTAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTc                                                                                                                    <  1:981467/151‑1 (MQ=255)
                                  gaacagTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGt                                                                                                                 <  1:1023730/151‑1 (MQ=255)
                                         tACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTattatt                                                                                                          <  1:957575/151‑1 (MQ=255)
                                          aCAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTg                                                                                                         <  1:954512/151‑1 (MQ=255)
                                            aaaTCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGaa                                                                                                      <  2:585740/151‑1 (MQ=255)
                                                 aaaCCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAg                                                                                                 <  1:578834/151‑1 (MQ=255)
                                                 aaaCCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAg                                                                                                 >  1:608461/1‑151 (MQ=255)
                                                   aCCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCa                                                                                               <  2:664652/151‑1 (MQ=255)
                                                          cAAGATGAAAAAGCGAATACGGATTTAACTTAAGCAGATAATAAAAATGTTGAAATTGGATTTGAGAGTAGTACAATTTATGGTGCAACAGGCAGACGATTAGATCATTTTTTTGTGGCAATTCAAAAAAAAAAAAAAAAAGCATATtata                                                                                        <  2:10864/151‑1 (MQ=255)
                                                           aaGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATaa                                                                                       <  1:48123/151‑1 (MQ=255)
                                                               tGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACat                                                                                   >  1:556557/1‑151 (MQ=255)
                                                                      gCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTc                                                                            >  1:355409/1‑151 (MQ=255)
                                                                        tGATACGGATTTAGCTTTAGCGGTTGATACAGCTCTCTTCCAGCGTTTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCat                                                                          >  1:515429/1‑151 (MQ=255)
                                                                          aTACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCatat                                                                        <  1:902236/151‑1 (MQ=255)
                                                                               gATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCGTATTATATTCATGATGTTCATATAGAAg                                                                   <  2:508335/151‑1 (MQ=255)
                                                                                aTTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGt                                                                  >  1:884875/1‑151 (MQ=255)
                                                                                  ttagctttagcGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTa                                                                <  2:1048302/151‑1 (MQ=255)
                                                                                       tttagcGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGAt                                                           <  1:1012435/151‑1 (MQ=255)
                                                                                        ttagcGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATGAACATCATGTTtataatatgcttttttcaa                                                          <  1:941384/151‑16 (MQ=255)
                                                                                        ttagcGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATc                                                          <  1:741599/151‑1 (MQ=255)
                                                                                         tagcGGTTGATAAAGCTGTTGCACTTGGTTTTGATATTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTAAAGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATca                                                         <  2:967151/151‑1 (MQ=255)
                                                                                            cGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATcaaca                                                      >  2:625190/1‑151 (MQ=255)
                                                                                                tgctAAAGCTGTTCAACCTGGTTTGGTTAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACaaaat                                                  <  1:152376/148‑1 (MQ=255)
                                                                                                 cgtccgcctgttgcaccataaattgtaatactatcaaatccaagtgcaacagGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACaaaata                                                 >  2:434594‑M2/53‑151 (MQ=255)
                                                                                                  aTAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACaaaataa                                                <  1:899972/151‑1 (MQ=255)
                                                                                                      gcctgttgcaccataaattgtaatactatcaaatccaagtgcaacagGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAAtt                                            >  1:788395‑M2/48‑151 (MQ=255)
                                                                                                         TGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGaa                                         <  2:653537‑M1/151‑1 (MQ=255)
                                                                                                               ACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTAATAAAAAAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTg                                   >  2:58427‑M1/1‑151 (MQ=255)
                                                                                                                 ttGGATTTGATAGTATTGTGATTAATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGcc                                 <  1:9268/151‑1 (MQ=255)
                                                                                                                 ttGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGcc                                 <  2:339913/151‑1 (MQ=255)
                                                                                                                          aTAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCCTAAAGGTCa                        <  1:408761/151‑1 (MQ=255)
                                                                                                                            aGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTAAAAAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCCTAAAGGTCAAc                      <  2:99743/151‑1 (MQ=255)
                                                                                                                                ttACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCCTAAAGGTCAACATa                   >  1:547846/1‑151 (MQ=255)
                                                                                                                                     aatccaagtgcaacagGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCCTAAAGGTCAACATACAGTTg             <  1:434594‑M2/135‑1 (MQ=255)
                                                                                                                                      tttATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCCTAAAGGTCAACATACAGTTGa            <  1:361538/151‑1 (MQ=255)
                                                                                                                                           ggTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCCTAAAGGTCAACATACAGTTGAAAAAg       >  1:466250/1‑151 (MQ=255)
                                                                                                                                           ggTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCATAAAGGTCAACATACAGTTGAATAAg       >  2:77044/1‑151 (MQ=255)
                                                                                                                                                aacagGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCCTAAAGGTCAACATACAGTTGAAAAAGATAag  <  1:21767‑M2/146‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TTCAGTCAGTAAAGAAGAACGCCAACTTCTAACAGAACAGTTACAAATCAAACCAGTTCAAGCTGAAAAAGCTGATACGGATTTAGCTTTAGCGGTTGATAAAGCTGTTGCACTTGGATTTGATAGTATTACAATTTATGGTGCAACAGGCGGACGATTAGATCACTTTTTTGGGGCAATTCAGTTATTATTGAAAAAAGCATATTATAAACATGATGTTCATATAGAAGTTATCGATCAACAAAATAAAATTGAATTATTGCCTAAAGGTCAACATACAGTTGAAAAAGATAAG  >  CP000730/1235994‑1236288

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