New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 87106216 (0.820)4 (0.230) 3/118 1.6 20.9% coding (96/921 nt) gsiC glutathione ABC transporter permease
?NC_000913 = 871080 16 (0.920)coding (114/921 nt) gsiC glutathione ABC transporter permease

ATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGG  >  NC_000913/871018‑871133
                                                              |                                                     
aTCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgg                                                  >  2:131685‑M2/1‑63 (MQ=255)
  cGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggc                                                >  1:270526‑M2/1‑61 (MQ=255)
   gTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGcc                                               >  1:245495/1‑68 (MQ=255)
     ctcGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgtcgggcagc                                             >  2:292424‑M2/1‑58 (MQ=255)
     ctcGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggcagc                                             >  1:37768‑M2/1‑58 (MQ=255)
           gCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTg                                       >  1:358126/1‑68 (MQ=255)
                tgtTTTTATTTGTCCATATGCTGCCCTGCGATCCGGCGCGATTGATTgccgggcagcatatggacaaa                                  >  1:516662‑M2/1‑47 (MQ=255)
                    tttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAgg                              <  1:633164/68‑1 (MQ=255)
                     ttATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGt                             <  1:549204/68‑1 (MQ=255)
                      tATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGATGCTGATGCGCAGGtt                            <  1:392223/68‑1 (MQ=255)
                      tATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGtt                            <  1:269939/68‑1 (MQ=255)
                      tATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGtt                            <  1:164539/68‑1 (MQ=255)
                      tATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGtt                            <  1:149467/68‑1 (MQ=255)
                        tttGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTgccgggcagcatatggacaaataaaaaca                          >  1:429025‑M2/1‑39 (MQ=255)
                         ttGTCCATATGCTGTCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGGtata                         <  1:18082/68‑1 (MQ=255)
                           gTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGa                       <  1:217998/68‑1 (MQ=255)
                               atatGCTGCCCGGCGAGCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGGTATAGAACTg                   <  1:141806/68‑1 (MQ=255)
                                  tgctgcCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTg                >  1:22196/1‑68 (MQ=255)
                                   gctgcACGGCGATCCGGCGAGATTGATTGCCGGGCCCGAAGCTGATGCGCTGGTTATAGAACTGGTGc               >  2:254517/1‑68 (MQ=255)
                                                ccGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCGGGTTATAGAACTGGTGCGTCAGCAGCTggg  >  1:279585/1‑68 (MQ=255)
                                                ccGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTggg  <  1:583211/68‑1 (MQ=255)
                                                ccGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGGCAGCAGCTggg  <  1:269228/68‑1 (MQ=255)
                                                              |                                                     
ATCGTCTCGGTGCTGGTGTTTTTATTTGTCCATATGCTGCCCGGCGATCCGGCGCGATTGATTGCCGGGCCCGAAGCTGATGCGCAGGTTATAGAACTGGTGCGTCAGCAGCTGGG  >  NC_000913/871018‑871133

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