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)51 (0.970) 42/144 0.2 97.9% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 1 (0.020)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 2774196NA (NA)38 (0.720) 32/144 0.7 97.2% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 1 (0.020)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

GGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2774125‑2774196
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524522‑4524596
                                                                                                                                                   
GGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGC                                                                         <  1:1454531/76‑1
 GCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCA                                                                        <  2:632545/76‑1
     CGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCC                                                                    >  1:1293834/1‑76
     CGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCC                                                                    >  2:928527/1‑76
       TCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCG                                                                  <  1:856183/76‑1
             CGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCT                                                            <  1:218873/76‑1
             CGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCT                                                            >  1:353792/1‑76
                   ACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAAT                                                      <  2:310535/76‑1
                       AGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGG                                                  >  2:1610925/1‑76
                             CGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGA                                            <  1:682464/76‑1
                                ATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCT                                         >  1:67572/1‑76
                                ATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCT                                         >  2:380478/1‑76
                                  TCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTC                                       >  1:39612/1‑76
                                        GATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTT                                 <  2:437958/76‑1
                                         ATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTC                                >  1:675471/1‑76
                                         ATGGTTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTC                                <  2:1141832/76‑1
                                            GTTCAACCCGTAAGTTAGCGCTTATGGGCCGCAACCGCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGT                             <  2:675453/76‑1
                                             TTCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTA                            >  2:71234/1‑76
                                              TCAACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAT                           >  2:318368/1‑76
                                                AACCCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCA                         <  2:454370/76‑1
                                                   CCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATA                      <  1:918449/76‑1
                                                   CCTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATA                      <  1:1697136/76‑1
                                                    CTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAA                     <  1:1380291/76‑1
                                                    CTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAA                     <  1:1097381/76‑1
                                                    CTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAA                     >  1:523738/1‑76
                                                    CTTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCAGATGCGACCTTCAAGGTTCAGTATCAATAA                     <  2:165318/76‑1
                                                     TTAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAC                    >  2:845415/1‑76
                                                      TAAGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACC                   <  2:1371432/76‑1
                                                        AGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTA                 >  2:1471233/1‑76
                                                        AGTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTA                 <  1:970475/76‑1
                                                         GTTAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAC                >  2:1490307/1‑76
                                                           TAGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCC              >  1:1614814/1‑76
                                                            AGCGCTTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCA             >  2:1406075/1‑76
                                                                 TTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTC        >  2:1652131/1‑76
                                                                 TTATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTC        >  2:1293599/1‑76
                                                                  TATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCA       <  2:1460109/76‑1
                                                                  TATGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCA       >  2:80799/1‑76
                                                                   ATGGGCCGCAACCCCGGGTGCNGCTAATGCGGATGCGAGCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAG      <  2:287139/76‑1
                                                                       GCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  <  1:839717/76‑1
                                                                                                                                                   
GGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGG‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑  >  REL606/2774125‑2774196
‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑‑gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524522‑4524596

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