New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? AP009048 = 359307947 (0.800)3 (0.060) 3/90 NT 7.3% coding (185/627 nt) dsbA periplasmic protein disulfide isomerase I
?AP009048 = 3593101 35 (0.680)coding (163/627 nt) dsbA periplasmic protein disulfide isomerase I

AGGTCACCACCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/3593004‑3593079
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatcagTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGTGAAGATGACTAAAT  <  AP009048/3593101‑3593032
                                                                                                                                                  
AGGTCACCACCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATT                                                                             <  1:5040255/71‑1
 GGTCACCACCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTA                                                                            <  1:2713652/71‑1
 GGTCACCACCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTA                                                                            <  1:3827041/71‑1
  GTCACCACCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTAT                                                                           >  1:1312381/1‑71
     ACCACCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCA                                                                         <  1:2348495/70‑1
            ATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATT                                                                             <  1:3028945/59‑1
                   TGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGTTTGAAGAAGTTCT                                                          >  1:1152078/1‑71
                                    ATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATG                                         >  1:4716061/1‑71
                                                          TTTTCTTCACATTATCAGTTTGAAGAAGTTCTGCATATTTC                                                 >  1:4903778/1‑41
                                                                       ATCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGTGAAG            <  1:4616743/65‑1
                                                                        TCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGTGAAGTTGACgt     >  1:2278596/1‑69
                                                                        TCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGTGAAGATGACTA     <  1:3417963/71‑1
                                                                        TCAGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTG                           <  1:3224933/49‑1
                                                                          AGTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGTGAAGATGACTAA    <  1:1346480/70‑1
                                                                           GTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGTGAAGATGACTAAAT  <  1:1500085/71‑1
                                                                                                                                                  
AGGTCACCACCCATGAAGTTGACGTGGTATTTAGTCATCTTCACGCCTTCCGGCAGTTTTTTCTTCACATTATCAG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  AP009048/3593004‑3593079
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑tatcagTTTGAAGAAGTTCTGCATATTTCTGATAATGTGAAGAAAAAACTGCCGGAAGGCGTGAAGATGACTAAAT  <  AP009048/3593101‑3593032

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