Predicted mutation
evidence seq id position mutation annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp coding (494‑499/549 nt) fimA → major type 1 subunit fimbrin (pilin)

New junction evidence
  seq id position reads (cov) reads (cov) score skew freq annotation gene product
* ? REL606 2772854 =NA (NA)27 (1.300) 26/144 0.0 100% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 0 (0.000)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 2774196NA (NA)23 (1.110) 20/144 0.2 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

GTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2772931‑2772854
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTC  <  REL606/4524527‑4524458
                                                                                                                                                    
GTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATG                                                                          >  2:591825/1‑76
GTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATG                                                                          >  1:672597/1‑76
 TTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGG                                                                         >  2:360576/1‑76
 TTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGG                                                                         >  1:487299/1‑76
  TATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG                                                                        >  1:233206/1‑76
  TATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG                                                                        >  1:415032/1‑76
  TATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG                                                                        <  1:456855/76‑1
        GAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGG                                                                  >  2:310035/1‑76
               ATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGT                                                           <  1:531216/76‑1
                TTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTT                                                          <  1:501427/76‑1
                 TCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTG                                                         >  2:574757/1‑76
                               AATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCC                                           >  1:382691/1‑76
                                  CATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGG                                        >  1:482307/1‑76
                                   ATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGA                                       >  1:249692/1‑76
                                   ATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGAAAAATAACGCGCCTGGA                                       >  1:534494/1‑76
                                    TATCATGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAA                                      >  2:191954/1‑76
                                        AGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGA                                  >  2:125174/1‑76
                                              AATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTA                            >  1:162226/1‑76
                                               ATACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTAT                           >  2:3609/1‑76
                                                TACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATT                          <  1:107100/76‑1
                                                 ACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTG                         >  1:142007/1‑76
                                                 ACCACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTG                         <  1:545118/76‑1
                                                    ACAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTT                      <  1:549532/76‑1
                                                     CAACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAAAGGTATTGGTTC                     >  1:439230/1‑76
                                                      AACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCC                    <  2:301886/76‑1
                                                      AACCCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCC                    >  2:34090/1‑76
                                                         CCTTAAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTT                 <  2:381200/76‑1
                                                             AAGTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTC             >  1:409759/1‑76
                                                               GTTAGCGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAG           <  2:26368/76‑1
                                                                    aGCTGATGGGTTGCGGCCCCGGTTGCGAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTG      <  1:433929/75‑1
                                                                    CGCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTG      >  1:280345/1‑76
                                                                     GCTTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGT     >  2:522621/1‑76
                                                                       TTATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTT   >  1:311111/1‑76
                                                                        TATGGGTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTC  >  2:360047/1‑76
                                                                                                                                                    
GTTATCGTGAGAGTAATTCATCGGCACGTTAAATCATATCAGGCGTAATACCACAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  <  REL606/2772931‑2772854
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggTTGCGGCCCCGGTTGCAAAATAACGCGCCTGGAACGGAATGGTATTGGTTCCGTTATTCAGGGTTGTTTC  <  REL606/4524527‑4524458

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