New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 279163 =NA (NA)17 (0.040) 12/194 25.7 3.0% noncoding (768/768 nt) IS1 repeat region
?NC_000913 = 3623367 552 (1.460)coding (4176/4236 nt) rhsB rhs element protein RhsB
Rejected: Coverage evenness skew score above cutoff.
Rejected: Frequency below/above cutoff threshold.

ACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/279262‑279163
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGACTTGATATTGTTATCTTTACTTTTAACTATTTTCTCATCTTTACC  <  NC_000913/3623367‑3623269
                                                                                                                                                                                                       
ACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCC                                                                                                    >  1:585638/1‑101
ACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCC                                                                                                    <  1:586154/101‑1
ACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCC                                                                                                    <  1:586155/101‑1
ACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCC                                                                                                    >  1:585407/1‑101
                                     CATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCG                                                               <  1:586156/101‑1
                                             AGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAA                                                       <  1:586157/101‑1
                                                     GGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGA                                               >  2:6807473/1‑101
                                                      GCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGAC                                              >  2:6807203/1‑101
                                                      GCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCC                                                                                                    <  2:3611480/47‑1
                                                       CATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCCCG                                                                                                 >  2:3611481/1‑49
                                                          TATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGACTTGA                                          <  1:6723133/101‑1
                                                           ATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGACTTGAT                                         >  2:6806898/1‑101
                                                                   ATAAAACACTATCAATAAGTTGGAGTCATTACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGACTTGATATTGTTAT                                 >  1:6807202/1‑101
                                                                                    AGTTGGAGTCATTACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGACTTGATATTGTTATCTTTACTTTTAACTATT                >  2:6807340/1‑101
                                                                                       TGGAGTCATTACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGACTTGATATTGTTATCTTTACTTTTAACTATTTTC             >  1:6807220/1‑101
                                                                                                TACCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGACTTGATATTGTTATCTTTACTTTTAACTATTTTCTCATCTTTA    >  1:6806959/1‑101
                                                                                                  CCCCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGACTTGATATTGTTATCTTTACTTTTAACTATTTTCTCATCTTTACC  >  2:6806933/1‑101
                                                                                                                                                                                                       
ACGGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCAGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/279262‑279163
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑CCCGTGGTTATCTATTGCTTTAAAATTACCACTGGGCGTACCTAAATAAAATGACTTGATATTGTTATCTTTACTTTTAACTATTTTCTCATCTTTACC  <  NC_000913/3623367‑3623269

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.