New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 = 3653230NA (NA)27 (0.950) 21/276 0.3 90.0% noncoding (1/1195 nt) IS5 repeat region
?NC_000913 = 3766089 3 (0.110)coding (3907/4134 nt) rhsA rhs element protein RhsA

TTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGC  >  NC_000913/3653102‑3653367
                                                                                                                                |                                                                                                                                         
ttCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGacac                                                                                                                                <  2:441118/140‑1 (MQ=33)
        tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAtttttt                                                                                                                        >  2:80051/1‑140 (MQ=40)
        tGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTAtttttt                                                                                                                        <  2:58505/140‑1 (MQ=40)
                cTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATaaa                                                                                                                >  1:402358/1‑140 (MQ=255)
                   aCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGaa                                                                                                             <  2:21107/140‑1 (MQ=255)
                   aCTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGaa                                                                                                             <  1:385345/140‑1 (MQ=255)
                          aaTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTa                                                                                                      >  1:92693/1‑140 (MQ=255)
                          aaTTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTa                                                                                                      >  1:257280/1‑140 (MQ=255)
                            ttCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAcc                                                                                                    >  1:85184/1‑140 (MQ=255)
                              cACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTAcccc                                                                                                  >  1:15533/1‑140 (MQ=255)
                                 tGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaa                                                                                               >  1:201415/1‑140 (MQ=255)
                                  gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                              <  2:15533/140‑1 (MQ=255)
                                  gTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCaaa                                                                                              <  1:102389/140‑1 (MQ=255)
                                   tCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAAt                                                                                             <  1:144295/140‑1 (MQ=255)
                                       cGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATtata                                                                                         <  2:236286/140‑1 (MQ=255)
                                            gTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAg                                                                                    <  2:60881/140‑1 (MQ=255)
                                              aaGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGAt                                                                                  <  1:224254/140‑1 (MQ=255)
                                                    atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACaa                                                                            >  2:374460/1‑140 (MQ=255)
                                                    atgaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACaa                                                                            >  2:66995/1‑140 (MQ=255)
                                                      gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAata                                                                          >  2:387750/1‑140 (MQ=255)
                                                      gaCTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAata                                                                          >  1:108648/1‑140 (MQ=255)
                                                             gatgaACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGac                                                                   >  2:201193/1‑140 (MQ=255)
                                                                      tGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATc                                                          <  2:92693/140‑1 (MQ=255)
                                                                         tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGaa                                                       <  1:66995/140‑1 (MQ=255)
                                                                         tCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGaa                                                       <  1:374460/140‑1 (MQ=255)
                                                                                  ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTAttt                                              >  2:288602/1‑140 (MQ=255)
                                                                                  ccAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTAttt                                              >  1:181025/1‑140 (MQ=255)
                                                                                      aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCaa                                          <  2:321878/140‑1 (MQ=255)
                                                                                      aTGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCaa                                          <  2:257280/140‑1 (MQ=255)
                                                                                                  aTCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCa                              >  2:90721/1‑140 (MQ=255)
                                                                                                   tctcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAt                             <  2:269514/140‑1 (MQ=255)
                                                                                                   tctcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAt                             <  2:179574/140‑1 (MQ=255)
                                                                                                    ctcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAtt                            >  1:137340/1‑140 (MQ=255)
                                                                                                    ctcATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCAtt                            >  1:132494/1‑140 (MQ=255)
                                                                                                              ggACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAAcc                  >  1:338512/1‑140 (MQ=255)
                                                                                                                              tccttAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGc  >  2:341839/1‑140 (MQ=255)
                                                                                                                                |                                                                                                                                         
TTCTGGTCTGACGGCGCTTACTGCTGAATTCACTGTCGGCGAAGGTAAGTTGATGACTCATGATGAACCCTGTTCTATGGCTCCAGATGACAAACATGATCTCATATCAGGGACTTGTTCGCACCTTCCTTAGGTGACACTATTTTTTTACATAAAGAACACCCTACCCCCAAATTATATAAAGATGGACAATATCTTGACAAACAAATCGAAGGTTTATTTTCAACATTATTAACCATTAATGACAACCTTTTACGAGCAAAAGC  >  NC_000913/3653102‑3653367

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