New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 4160055179 (1.920)3 (0.030) 3/184 NT 1.8% coding (736/1401 nt) sthA pyridine nucleotide transhydrogenase, soluble
?NC_000913 = 4160081 159 (1.850)coding (710/1401 nt) sthA pyridine nucleotide transhydrogenase, soluble
Rejected: Frequency below/above cutoff threshold.

TTGACCTTCAGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4159904‑4160055
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgtcacaACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATACCGATTCGCTGGCGTTACAGAACATTGGGC  <  NC_000913/4160081‑4159938
                                                                                                                                                                                                                                                                                                        
TTGACCTTCAGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGT                                                                                                                                                      <  4:392382/148‑1
      TTCAGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAAC                                                                                                                                                >  3:667128/1‑148
                         CAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCAC                                                                                                                                                   <  2:1046297/126‑1
                         CAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCAC                                                                                                                                                   >  1:1046297/1‑126
                                       ATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGAT                                                                                                               <  4:667128/148‑1
                                                    CCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGT                                                                                                                                                      <  3:685470/96‑1
                                                    CCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGT                                                                                                                                                      >  4:685470/1‑96
                                                                                        AGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACG                                                                                                                 >  1:271097/1‑97
                                                                                        AGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACG                                                                                                                 <  2:271097/97‑1
                                                                                                                                                TCGTCACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATACCGATTCGCTGGCGTTACAGAACATT      <  1:244540/148‑1
                                                                                                                                                TCGTCACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATACCGATTCGCTGGCGTTACAGAACATT      <  2:612288/148‑1
                                                                                                                                                TCGTCACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATACCGATTCGCTGGCGTTAC              >  3:680032/1‑140
                                                                                                                                                TCGTCACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATACCGATTCGCTGGCGTTAC              <  4:680032/140‑1
                                                                                                                                                   TCACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGT                                             >  4:872135/1‑106
                                                                                                                                                   TCACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGT                                             <  3:872135/106‑1
                                                                                                                                                    CACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATACCGATTCGCTGGCGTTACAGAACATTGGGC  >  1:62561/1‑148
                                                                                                                                                    CACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATACCGATTCGCTGGCGTTACAGAACATTGG    <  3:160386/146‑1
                                                                                                                                                    CACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATACCGATTCGCTGGCGTTACAGAACATTGG    >  4:160386/1‑146
                                                                                                                                                    CACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATAC                               <  3:1144889/119‑1
                                                                                                                                                    CACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATAC                               >  4:1144889/1‑119
                                                                                                                                                    CACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAAC                                                                             >  2:256667/1‑73
                                                                                                                                                    CACAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAAC                                                                             <  1:256667/73‑1
                                                                                                                                                      CAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGA                                                                                                        >  2:74295/1‑44
                                                                                                                                                      CAACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGA                                                                                                        <  1:74295/44‑1
                                                                                                                                                       AACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAA                                                                              <  2:457939/69‑1
                                                                                                                                                       AACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAA                                                                              >  2:132338/1‑69
                                                                                                                                                       AACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAA                                                                              >  1:457939/1‑69
                                                                                                                                                       AACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAA                                                                              <  1:132338/69‑1
                                                                                                                                                                                                                                                                                                        
TTGACCTTCAGCTGTCCGCGGCTGTCAGTTTCTAGCCCAATGTTCTGTAACGCCAGCGAATCGGTATTACCGGTGCGACCGTTGGCATAGAGCAGGCAGTCAGCTTTCAGTTTTTTACCCGACTTCAGATGCATGATCACACCATCGTCACA‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  NC_000913/4159904‑4160055
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tcgtcacaACGAAGAGTACGAGAAGATCGAAGGCTGTGACGATGGTGTGATCATGCATCTGAAGTCGGGTAAAAAACTGAAAGCTGACTGCCTGCTCTATGCCAACGGTCGCACCGGTAATACCGATTCGCTGGCGTTACAGAACATTGGGC  <  NC_000913/4160081‑4159938

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