New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1307027 =7 (0.270)3 (0.130) 3/114 2.2 22.2% coding (124/330 nt) yciU UPF0263 family protein
?NC_000913 1307124 = 15 (0.670)coding (27/330 nt) yciU UPF0263 family protein

GATCGAGATTGTCGGCAGCCAGTTCGAGAAAAATATCGTAGGCCTGTTCAAGCGTTTCATCTTCAGTCAGGCGATTGTTTAGATCCATATCCATAGTTACTACCTGTTTAACCTCTGTTGGCGACGT  >  NC_000913/1307057‑1307183
                                                                   |                                                           
gatatttttctcgagctggctgccgacactctcgatccagctgatgttctgctgttcaatcttcagtc                                                             >  2:654046‑M2/68‑68 (MQ=255)
          gTCGGCCGCCAGTTCGAGAAAAATATCGTAGGCCTGTTCAAGCGTTTCATCTTCAGTCAGGCGATTGt                                                   <  2:735754/68‑1 (MQ=255)
          gTCGGCAGCCAGTTCGAGAAAAATATCGTAGGCCTGTTCAAGCGTTTCATCTTCAGTCAGGCGATTGt                                                   <  2:16450/68‑1 (MQ=255)
            cGGCAGCCAGTTCGAGAAAAATATCGTAGGCCGGTTAAAGCGTGTCATATTCAGTCAGGCGAATGttt                                                 >  1:387279/1‑68 (MQ=255)
                aGCCAGTTCGAGTAAATAATCTTAGGCCTGTTCAAGCGGTTCATCTTCAGGCAGGCTATTGTTTAGat                                             <  2:343242/68‑1 (MQ=255)
                         gagaAAAATATCATAGGCCTGTTCAAGCGTTTCATCTTCAGTCAGGCGATTGTTTAGatccatatcca                                    <  2:727855/68‑1 (MQ=255)
                             aaaaTATCGTAGGCCTGTTCAAGCGTTTCATCTTCAGTCAGGCGATTGTTTAGATCCATATCCATAGt                                <  2:310786/68‑1 (MQ=255)
                             aaaaTATCGTAGGCCTGTTCAAGCGCTCCATCTTCAGTCAGGCGATCGCCTAGACCCATATCCATAgc                                <  2:232831/68‑2 (MQ=255)
                               tcgatccagctgatgttctgctgttcaatcttcagtCAGGCGATTGTTTAGATCCATATCCATAGTta                              <  1:654046‑M2/32‑1 (MQ=255)
                                cgatccagctgatgttctgctgttcaatcttcagtCAGGCGATTGTTTAGATCCATATCCATAGTtac                             >  2:51019‑M2/36‑68 (MQ=255)
                                  atCGTAGGCCTGTTAAAGCGTTTCATCTTCAGTCAGGCGATTGTTTAGATCCATATCCATAGTtacta                           >  1:939502/1‑68 (MQ=255)
                                     gTAGGCCTGTTCAAGCGTTTCATCTTCAGTCAGGCg                                                        >  1:288114/1‑36 (MQ=255)
                                     gTAGGCCTGTTCAAGCGTTTCATCTTCAGTCAGGCg                                                        <  2:288114/36‑1 (MQ=255)
                                          ccTGTTCAAGCGTTTCATCTTCAGTCAGGCGATTGTTTAGATCCATATCCATAGTTACTACCTGTTTa                   >  1:334408/1‑68 (MQ=255)
                                          ccTGTTCAAGCGTTTCATCTTCAGTCAGGCGATTGTTTAGATCCATATCCATAGTTACTACCTGTTTa                   >  2:788837/1‑68 (MQ=255)
                                                cAAGCGTTTCATCTTCAGTCAGGCGATTGTGTAGATCCATATCCATAGTTACTACCTGTTTAACctct             <  2:82926/68‑1 (MQ=255)
                                                  aGCGTTTCATCTTCAGTCAGGCGATTGTTTAGATCCATATCCATAGTTACTACCTGTTTAACCTCTGt           >  1:200502/1‑68 (MQ=255)
                                                    cGTTTCATCTTCAGTCAGGCGATTGTTTAGATCCATATCCATAGTTACTACCTGTTTAACCTTGTTgg        >  1:931191/1‑68 (MQ=255)
                                                           tcttcagtCAGGCGATTGTTTAGATCCATATCCATAGTTACTACCTGTTTAACCTCTGTTGGCGACGt  <  2:604982‑M2/60‑1 (MQ=255)
                                                                   |                                                           
GATCGAGATTGTCGGCAGCCAGTTCGAGAAAAATATCGTAGGCCTGTTCAAGCGTTTCATCTTCAGTCAGGCGATTGTTTAGATCCATATCCATAGTTACTACCTGTTTAACCTCTGTTGGCGACGT  >  NC_000913/1307057‑1307183

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.