Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 4,542,238 IS1 (–) +8 bp 5.0% coding (202‑209/597 nt) fimE → regulator for fimA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 1978503 =NA (NA)9 (0.080) 9/176 NT 6.2% noncoding (768/768 nt) IS1 repeat region
?NC_000913 = 4542245 123 (1.170)coding (209/597 nt) fimE regulator for fimA
* ? NC_000913 3583428 =NA (NA)4 (0.040) 4/174 NT 2.9% noncoding (1/768 nt) IS1 repeat region
?NC_000913 4542238 = 122 (1.180)coding (202/597 nt) fimE regulator for fimA
Rejected: Frequency below/above cutoff threshold.

GGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1978600‑1978503
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAAAATCCGTTCTTCAGTCGGCGAATATTTATTCTACCTTCATTAAGGTCAAGGTCCTGATAATGCAGATCAAGCAGTTCACTAATACGCAT  <  NC_000913/4542245‑4542154
                                                                                                                                                                                              
GGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCG                                                                                             >  2:1343004/1‑99
 GAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGA                                                                                            <  1:133628/99‑1
 GAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGA                                                                                            <  1:2215973/99‑1
 GAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGA                                                                                            >  1:2560444/1‑99
   AGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAAA                                                                                          <  2:385987/99‑1
     TCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAAAAT                                                                                        >  1:779981/1‑99
        CTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAAAATCCG                                                                                     <  1:2211287/99‑1
                     AATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAAAATCCGTTCTTCAGTCGGC                                                                        <  2:829141/99‑1
                                 TGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAAAATCCGTTCTTCAGTCGGCGAATATTTATTC                                                            <  2:2940344/99‑1
                                                       TTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACCGAAAATCCGTTCTTCAGTCGGCGAATATTTATTCTACCTTCATTAAGGTCAAGGTC                                      <  1:881383/99‑1
                                                                         CTATCAATAAGTTGGAGTCATTACCGAAAATCCGTTCTTCAGTCGGCGAATATTTATTCTACCTTCATTAAGGTCAAGGTCCTGATAATGCAGATCAAG                    >  2:281596/1‑99
                                                                                  AGTTGGAGTCATTACCGAAAATCCGTTCTTCAGTCGGCGAATATTTATTCTACCTTCATTAAGGTCAAGGTCCTGATAATGCAGATCAAGCAGTTCACT           >  1:2052533/1‑99
                                                                                     TGGAGTCATTACCGAAAATCCGTTCTTCAGTCGGCGAATATTTATTCTACCTTCATTAAGGTCAAGGTCCTGATAATGCAGATCAAGCAGTTCACTAAT        <  1:1712139/99‑1
                                                                                           CATTACCGAAAATCCGTTCTTCAGTCGGCGAATATTTATTCTACCTTCATTAAGGTCAAGGTCCTGATAATGCAGATCAAGCAGTTCACTAATACGCAT  <  2:1120156/99‑1
                                                                                                                                                                                              
GGAAGTCGCTGTCGTTCTCAAAATCGGTGGAGCTGCATGACAAAGTCATCGGGCATTATCTGAACATAAAACACTATCAATAAGTTGGAGTCATTACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/1978600‑1978503
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑GAAAATCCGTTCTTCAGTCGGCGAATATTTATTCTACCTTCATTAAGGTCAAGGTCCTGATAATGCAGATCAAGCAGTTCACTAATACGCAT  <  NC_000913/4542245‑4542154

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