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)47 (1.150) 36/144 0.2 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)32 (0.780) 26/144 0.7 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

ACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2774119‑2774196
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCA  >  REL606/4524522‑4524591
                                                                                                                                                    
ACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATG                                                                          <  1:1053472/76‑1
ACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATG                                                                          >  2:1146556/1‑76
 CGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGG                                                                         <  1:740595/76‑1
  GCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG                                                                        <  2:1238872/76‑1
  GCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG                                                                        >  1:403517/1‑76
        CATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAA                                                                  >  1:1261159/1‑76
         ATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCGGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAAC                                                                 <  2:1253017/76‑1
         ATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGTTTATCGGCCGCAAC                                                                 >  2:81512/1‑76
         ATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAAC                                                                 >  2:899950/1‑76
              CTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGG                                                            >  2:393043/1‑76
                TAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGT                                                          <  2:691654/76‑1
                 AACGAGGCACCGAGNCGTCGCATTCTTCAGATGGTTCAACCCTTAAGGTAGCGCTTATGGGCCGCACCTCNGGGTG                                                         <  2:1340299/76‑1
                     AGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACACTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGC                                                     >  2:1065080/1‑76
                       GCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTA                                                   <  2:471468/76‑1
                       GCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTA                                                   <  2:865408/76‑1
                         ACCGAGGCGTCGCATTCTTCAGTTGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAAT                                                 >  2:471946/1‑76
                            GAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCG                                              >  2:422858/1‑76
                              GGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGA                                            <  2:1185714/76‑1
                               GCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGAT                                           <  2:450479/76‑1
                               GCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGAT                                           <  2:813726/76‑1
                               GCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGAT                                           <  1:1278554/76‑1
                                CGTCGCATTCTTCAGGTGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATG                                          <  1:1109026/76‑1
                                     CATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGNATGCGACC                                     <  1:1289746/76‑1
                                       TTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTT                                   <  1:157806/76‑1
                                         CTTCAGATGGTTCAACCCTTAAGTTAGCGCTTTTGGGCCGCAACCCCGGGTGCTGCTAATGCGGGTGCGACCTTCA                                 <  1:79372/76‑1
                                           TCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG                               <  1:562237/76‑1
                                            CAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGG                              <  1:1275587/76‑1
                                             tGATGGTTCAACCATTAAGTTAGCGCTTATGGGCCGCATCCCCGGGTGCTGCTAATGCGGGCGCTACCTTCACGGT                             >  2:17078/2‑76
                                                     CAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATC                     <  2:1197086/76‑1
                                                     CAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATC                     <  1:727125/76‑1
                                                        CCCTTAAGCTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGTGGATGCGACTTTCATGGTTCAGTATCAAT                  >  2:20974/1‑76
                                                         CCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATA                 <  2:696308/76‑1
                                                             AAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCT             <  2:932776/76‑1
                                                              AGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTA            <  2:679142/76‑1
                                                              AGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTA            >  1:417687/1‑76
                                                                      CTTATGAGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTT    >  2:741443/1‑76
                                                                        TATGGGCCGCCACCCCGAGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATATCCTACCCAGGTTCA  >  1:112438/1‑76
                                                                                                                                                    
ACGCGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2774119‑2774196
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCA  >  REL606/4524522‑4524591

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