New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 218516241 (1.190)3 (0.090) 3/262 4.1 7.2% coding (140/789 nt) thiM hydoxyethylthiazole kinase
?NC_000913 = 2185172 38 (1.170)coding (130/789 nt) thiM hydoxyethylthiazole kinase
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

AGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCA  >  NC_000913/2185016‑2185300
                                                                                                                                                  |                                                                                                                                          
aGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGAt                                                                                                                                          <  4:126060/149‑1 (MQ=255)
          agagCTTTTTGCTTGCTCAACTGCAGCACGCATCGCCTGAGCGAGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCATCGAgcg                                                                                                                                <  7:106115/149‑1 (MQ=255)
               tttttGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAg                                                                                                                            >  3:40985/1‑148 (MQ=255)
                  ttgcttgcTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCaa                                                                                                                         <  7:90238/148‑1 (MQ=255)
                  ttgcttgcTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAg                                                                                                                        >  5:109544/1‑149 (MQ=255)
                         ctCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTgg                                                                                                                 >  1:32973/1‑149 (MQ=255)
                           cAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCg                                                                                                               >  4:38253/1‑149 (MQ=255)
                            aaCGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGCTGCACCGAGCGCCAGCAAGGTATTGGCgg                                                                                                              >  7:122440/1‑149 (MQ=255)
                            aaCGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGc                                                                                                                <  6:133820/147‑1 (MQ=255)
                                     aCGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGttt                                                                                                     <  3:38253/149‑1 (MQ=255)
                                       gCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGc                                                                                                   <  2:32973/149‑1 (MQ=255)
                                         aTCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTAtt                                                                                                                   <  3:89637/131‑1 (MQ=255)
                                         aTCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTAtt                                                                                                                   >  4:89637/1‑131 (MQ=255)
                                            gCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGcacca                                                                                               <  6:109544/148‑1 (MQ=255)
                                                   cgcgTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTg                                                                                       <  4:40985/149‑1 (MQ=255)
                                                          cTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGc                                                                                <  2:125215/149‑1 (MQ=255)
                                                               tCAGTGTGCAAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACAGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTg                                                                            >  6:82786/1‑148 (MQ=255)
                                                                     tgCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAaga                                                                     <  5:82786/149‑1 (MQ=255)
                                                                       ccAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAgg                                                                   >  4:56749/1‑149 (MQ=255)
                                                                          aCGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAgggg                                                                 >  4:59948/1‑148 (MQ=255)
                                                                               aaTCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGt                                                            <  8:122440/148‑1 (MQ=255)
                                                                                   accaaGGCACTGGCGATAGCCGCAAACTGACTGGCCGCTTCGGTTTCGTTAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGt                                                        <  8:72064/146‑1 (MQ=255)
                                                                                     caaGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGa                                                      >  8:29608/1‑148 (MQ=255)
                                                                                             tGGCGATAGCCGCAAGCTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtg                                              >  5:30685/1‑148 (MQ=255)
                                                                                              ggCGATAGCCGCAAACTGCCTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGgtgt                                             >  8:130987/1‑148 (MQ=255)
                                                                                                      ccGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGt                                                                                                       >  6:144815/1‑82 (MQ=255)
                                                                                                      ccGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGt                                                                                                       <  5:144815/82‑1 (MQ=255)
                                                                                                        gCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTg                                                                            <  4:29485/107‑1 (MQ=255)
                                                                                                        gCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTg                                                                            >  3:29485/1‑107 (MQ=255)
                                                                                                         cAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTgcg                                   >  6:11973/1‑147 (MQ=255)
                                                                                                                 ctgGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGGGTAACGCGTGCGCAGATTgcgc                         <  3:59948/149‑1 (MQ=255)
                                                                                                                    gCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTgcgc                         <  5:14125/146‑1 (MQ=255)
                                                                                                                    gCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTgcgc                         >  6:14125/1‑146 (MQ=255)
                                                                                                                    gCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGa                      <  7:29608/149‑1 (MQ=255)
                                                                                                                       tctTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAg                                                 <  1:30695/119‑1 (MQ=255)
                                                                                                                       tctTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAg                                                 >  2:30695/1‑119 (MQ=255)
                                                                                                                            ggTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCcagc               <  6:30685/148‑1 (MQ=255)
                                                                                                                                   aTAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGTCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTgcg                                   <  3:63367/121‑1 (MQ=255)
                                                                                                                                   aTAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTgcg                                   >  4:63367/1‑121 (MQ=255)
                                                                                                                                   aTAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGAc       >  6:105056/1‑149 (MQ=255)
                                                                                                                                         aTCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCa  >  1:158840/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                          
AGGGTGTTTGAGAGCTTTTTGCTTGCTCAACGGCAGCACGCATCGCCTGAGCGCGTGGCTGCGTCAGTGTGCCAACGTTAATCAACAAGGCACTGGCGATAGCCGCAAACTGACTGGCCTCTTCGGTTTCGATAACCATCGCTGGCGATGCACCGAGCGCCAGCAAGGTATTGGCGGTAAAGGTTTGCACCACATCATTGGTCATGCAGTGCACAAGAGGGGAATGTTGGTGAAAAAGGTGTAACGCGTGCGCAGATTGCGCTGAACCCAGCAGGTCGACTTGCA  >  NC_000913/2185016‑2185300

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.