New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1411899 =52 (1.410)7 (0.240) 7/204 2.7 19.1% coding (50/936 nt) ttcA tRNA s(2)C32 thioltransferase, iron sulfur cluster protein
?NC_000913 = 1434958 18 (0.610)pseudogene (51/51 nt) ttcC pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related
Rejected: Frequency below/above cutoff threshold.

GGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTG  >  NC_000913/1434818‑1435103
                                                                                                                                            |                                                                                                                                                 
ggCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcag                                                                                                                                           >  2:77389‑M2/1‑141 (MQ=255)
   aaaaCGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggtt                                                                                                                                        <  6:97705‑M2/149‑12 (MQ=255)
            cAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgtt                                                                                                                                >  3:124309‑M2/1‑129 (MQ=255)
                      aaaaaaggggcataaaaAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTAc                                                                                                                                                   <  5:15367/119‑1 (MQ=255)
                      aaaaaaggggcaaaaaaAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTAc                                                                                                                                                   >  6:15367/1‑119 (MQ=255)
                               gcaaaaaaAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttctttctttgtaatttgttga                                                                                                            <  1:178956‑M2/149‑40 (MQ=255)
                               gcaaaaaaAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGa                                                                                                            >  1:157420/1‑149 (MQ=255)
                                    aaaaGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttctttctttgtaatttgttgattttc                                                                                                       >  2:178956‑M2/1‑105 (MQ=255)
                                           tGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTAATGTTGTATTCTTCTTTCTTACAAATTTCTTGATTTTTCGACAtt                                                                                                >  1:139127/1‑149 (MQ=255)
                                                           aTTGTTGTCGCTGATGATTTTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttctttctttgtaatttgttgattttcttgcattatttcagttctctgg                                                                                 <  8:121965‑M2/148‑67 (MQ=255)
                                                                 gtCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTta                                                                           >  6:146652/1‑148 (MQ=255)
                                                                    gCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttcttt                                                                                                                            <  6:134832‑M2/96‑24 (MQ=255)
                                                                    gCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttcttt                                                                                                                            >  5:134832‑M2/1‑73 (MQ=255)
                                                                    gCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttctt                                                                                                                             <  8:121938‑M2/95‑23 (MQ=255)
                                                                    gCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttctt                                                                                                                             >  7:121938‑M2/1‑73 (MQ=255)
                                                                    gCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTatat                                                                        <  5:146652/148‑1 (MQ=255)
                                                                     cTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATgg                                                                      <  5:97705/149‑1 (MQ=255)
                                                                       gatgatATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                                                                             >  4:111369/1‑140 (MQ=255)
                                                                       gatgatATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGt                                                                             <  3:111369/140‑1 (MQ=255)
                                                                           atttTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttctttctttgtaatttgttgattttcttgcattatttcagttctctggtactaaatggggcaaa                                                                 >  7:121965‑M2/1‑66 (MQ=255)
                                                                           atatTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCag                                                                <  2:157420/149‑1 (MQ=255)
                                                                            tatTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCag                                                                >  6:39123/1‑148 (MQ=255)
                                                                            tatTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCACTTACttgttcaggttgtattgttcttt                                                                                                                            <  1:105217‑M2/88‑24 (MQ=255)
                                                                            tatTTAACAAATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttcttt                                                                                                                            >  2:105217‑M2/1‑65 (MQ=255)
                                                                                tAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCagaaga                                                            <  4:58642/148‑1 (MQ=255)
                                                                                 aaCACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttctttctttgtaatttgttgattttcttgcattatttcagttctctggtactaaatggggcaaattgggg                                                           >  7:54276‑M2/1‑60 (MQ=255)
                                                                                            aaaTAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACttgttcaggttgtattgttctttctttgtaatttgttgattttcttgcattatttcagttctctggtactaaatggggcaaattgggggcaaactttgc                                                <  7:87395‑M2/148‑100 (MQ=255)
                                                                                               tAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCgg                                            <  5:39123/149‑1 (MQ=255)
                                                                                                    aaaaTACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGc                                       >  8:33295/1‑149 (MQ=255)
                                                                                                                  tCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGaaa                          >  5:57646/1‑148 (MQ=255)
                                                                                                                       gagTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGat                    <  7:33295/149‑1 (MQ=255)
                                                                                                                       gagTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGa                     <  2:139127/148‑1 (MQ=255)
                                                                                                                           tgtgATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGcc                                      >  8:119300/1‑127 (MQ=255)
                                                                                                                           tgtgATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGcc                                      <  7:119300/127‑1 (MQ=255)
                                                                                                                            gtgATTTTTGCCCTTACttgttcaggttgtattgttctttctttgtaatttgttgattttcttgcattatttcagttctctggtactaaatggggcaaattgggggcaaactttgcaactacgataaccgcgcattcaacatggctat                >  8:87395‑M2/1‑17 (MQ=255)
                                                                                                                                         ttACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCgtg  >  1:100260/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATGATATTTAACACATTGAAAAATAAGTAAAATACTTATGAGTCAGAGAGTTGTGATTTTTGCCCTTACTTGTTCAGGTTGTATTGTTCTTTCTTACTAATTTCTTGATTTTGCGACATTTAAAAGCGACTCAATTCGTTATATGGCATCAGAAGAGTATGCGTCATGCCGGAACGCCCAGCATAAGAAATCTGATATAAAAAACTGTGGCGTG  >  NC_000913/1434818‑1435103

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