New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 38884311 (0.010)32 (0.390) 22/96 2.1 98.5% intergenic (+239/‑4) mnmE/tnaC 5‑carboxymethylaminomethyluridine‑tRNA synthase GTPase subunit/tnaAB operon leader peptide
?NC_000913 3906599 = 0 (0.000)intergenic (‑32/+254) bglG/phoU transcriptional antiterminator BglG/negative regulator of the pho regulon

TTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAG  >  NC_000913/3906555‑3906648
                                            |                                                 
tctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACtt                                              >  1:4493296‑M2/45‑50 (MQ=255)
tctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACtt                                              >  1:4855047‑M2/45‑50 (MQ=255)
tctgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACtt                                              >  1:538666‑M2/45‑50 (MQ=255)
  tgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACtt                                              <  1:5276437‑M2/6‑1 (MQ=255)
  tgagtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCa                                            >  1:3215640‑M2/43‑50 (MQ=255)
    agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCa                                            <  1:2642385‑M2/8‑1 (MQ=255)
    agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCAcc                                          <  1:3980679‑M2/10‑1 (MQ=255)
    agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCAcc                                          <  1:2802856‑M2/10‑1 (MQ=255)
    agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCAc                                           <  1:4438928‑M2/9‑1 (MQ=255)
    agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCAc                                           <  1:5144535‑M2/9‑1 (MQ=255)
    agtgtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCAc                                           <  1:6021024‑M2/9‑1 (MQ=255)
       gtaataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAg                                        >  1:5070692‑M2/38‑49 (MQ=255)
         aataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTAt                                     >  1:1541704‑M2/36‑50 (MQ=255)
         aataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTAt                                     >  1:6284988‑M2/36‑50 (MQ=255)
          ataatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTAtt                                    >  1:5443886‑M2/35‑50 (MQ=255)
           taatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtc                                   >  1:4612308‑M2/34‑50 (MQ=255)
           taatgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtc                                   >  1:1390341‑M2/34‑50 (MQ=255)
             atgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtctc                                 >  1:3399547‑M2/32‑50 (MQ=255)
             atgtagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATtct                                  <  1:4873928‑M2/18‑1 (MQ=255)
                tagcctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTgg                              >  1:4526682‑M2/29‑50 (MQ=255)
                   cctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTa                           >  1:6357883‑M2/26‑50 (MQ=255)
                    ctcgtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTAt                          >  1:4754394‑M2/25‑50 (MQ=255)
                       gtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATg                         <  1:2556130‑M2/27‑1 (MQ=255)
                       gtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGt                        <  1:5239942‑M2/28‑1 (MQ=255)
                       gtgtcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGt                        <  1:3030215‑M2/28‑1 (MQ=255)
                          tcttgcgaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAgg                    <  1:3279115‑M2/32‑1 (MQ=255)
                                gaggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGc              >  1:2366512‑M2/13‑50 (MQ=255)
                                  ggataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCt            <  1:239083‑M2/40‑1 (MQ=255)
                                    ataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGc          <  1:2336481‑M2/42‑1 (MQ=255)
                                    ataagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGc          <  1:853724‑M2/42‑1 (MQ=255)
                                     taagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCg         <  1:3181115‑M2/43‑1 (MQ=255)
                                     taagtgcCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCg         <  1:1036936‑M2/43‑1 (MQ=255)
                                            cGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAg  >  1:5546144/1‑50 (MQ=255)
                                            cGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAg  <  1:1834560/50‑1 (MQ=255)
                                            |                                                 
TTTGCATGTTCATAGCAAGGACCTTTTTTATAAACAAAAAAACCCGACTTCACCAGTATTCTCTGGTTATGTCAGGTTTTGCCTGCGAATGCAG  >  NC_000913/3906555‑3906648

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