Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC NC_000913 4,542,065 IS1 (+) +8 bp 8.2% coding (29‑36/597 nt) fimE → regulator for fimA

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? NC_000913 3583428 =NA (NA)12 (0.090) 12/176 NT 8.5% noncoding (1/768 nt) IS1 repeat region
?NC_000913 = 4542072 117 (0.920)coding (36/597 nt) fimE regulator for fimA
* ? NC_000913 = 3584195NA (NA)9 (0.060) 9/176 NT 6.5% noncoding (768/768 nt) IS1 repeat region
?NC_000913 4542065 = 117 (0.920)coding (29/597 nt) fimE regulator for fimA

ACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3583515‑3583428
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGT  <  NC_000913/4542072‑4541976
                                                                                                                                                                                         
ACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCAACTTCTTTAC                                                                                        <  2:3632666/99‑1
                    AAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTA                                                                    <  2:1796139/99‑1
                            GCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAAT                                                            >  2:1246353/1‑99
                                 TGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAG                                                       >  1:652090/1‑99
                                                  ACCATCATACACTAAATCAGTAAGTTGGCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATA                                      <  1:3626557/99‑1
                                                   CCATCATACACTAAATCAGTAAGTTGGCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATAT                                     >  2:1101354/1‑99
                                                          ACACTAAATCAGTAAGTTGGCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATATATTTGAA                              <  2:3635721/99‑1
                                                            ACTAAATCAGTAAGTTGGCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATATATTTGAATG                            >  1:3558533/1‑99
                                                              TAAATCAGTAAGTTGGCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATATATTTGAATGAA                          <  2:2434774/99‑1
                                                                             GCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATATATTTGAATGAATACCTATAGGAAACC           <  2:1880573/99‑1
                                                                             GCAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATATATTTGAATGAATACCTATAGGAAACC           >  2:3519313/1‑99
                                                                              CAGCATCACCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATATATTTGAATGAATACCTATAGGAAACCT          >  1:1427907/1‑99
                                                                                      CCAACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGT  <  1:374894/99‑1
                                                                                                                                                                                         
ACAGGAGGGACAGCTGATAGAAACAGAAGCCACTGGAGCACCTCAAAAACACCATCATACACTAAATCAGTAAGTTGGCAGCATCACC‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  NC_000913/3583515‑3583428
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑AACTTCTTTACCGGTAAGATAACGACGTTTACTCACAATAGTAGTACTCCTAACTGAGATATATTTGAATGAATACCTATAGGAAACCTCAATCGGT  <  NC_000913/4542072‑4541976

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.