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)83 (1.350) 58/144 0.1 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)61 (0.990) 48/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)

CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTG  >  REL606/2774122‑2774270
                                                                          |                                                                          
cgAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGt                                                                           <  1:1826496/76‑1 (MQ=31)
  aGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGttt                                                                         >  2:1201396/1‑76 (MQ=34)
   ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                        <  2:1461970/76‑1 (MQ=35)
     cATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAc                                                                      <  2:1691902/76‑1 (MQ=39)
     cATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCACATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAc                                                                      >  1:279205/1‑76 (MQ=255)
      aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                     <  2:1958495/76‑1 (MQ=40)
          tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTGCCCCcg                                                                 >  1:1318149/1‑76 (MQ=255)
            ttAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgcg                                                               <  2:1531546/76‑1 (MQ=255)
              aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtg                                                             >  2:1562276/1‑76 (MQ=255)
                 gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                          >  1:95722/1‑76 (MQ=255)
                 gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                          >  1:1689560/1‑76 (MQ=255)
                 gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                          >  1:1689572/1‑76 (MQ=255)
                     cACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTa                                                      <  1:1422458/76‑1 (MQ=255)
                       ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTAct                                                    <  1:480583/76‑1 (MQ=255)
                        cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActt                                                   <  1:1939746/76‑1 (MQ=255)
                        cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCAGCGTGGTGGCTActt                                                   >  2:1202527/1‑76 (MQ=255)
                          aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtct                                                 >  1:808919/1‑76 (MQ=255)
                          aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtct                                                 <  1:1073124/76‑1 (MQ=255)
                            gCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACttcttc                                               <  1:642084/76‑1 (MQ=255)
                             cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCg                                              <  1:1445069/76‑1 (MQ=255)
                               tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAt                                            >  1:1235497/1‑76 (MQ=255)
                               tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAt                                            <  1:58506/76‑1 (MQ=255)
                               tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAt                                            >  2:154409/1‑76 (MQ=255)
                                cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAtt                                           >  2:1483175/1‑76 (MQ=255)
                                cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAtt                                           >  1:1115682/1‑76 (MQ=255)
                                  cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCg                                         >  2:1204747/1‑76 (MQ=255)
                                  cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCg                                         >  2:1485509/1‑76 (MQ=255)
                                  cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCg                                         >  2:233123/1‑76 (MQ=255)
                                   aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGt                                        >  1:1841010/1‑76 (MQ=255)
                                   aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGt                                        >  2:1336183/1‑76 (MQ=255)
                                    ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                       <  1:90151/76‑1 (MQ=255)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                      >  2:964633/1‑76 (MQ=255)
                                     tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCAGCGTGGTGGCTACTTCTTCGATTCGTaa                                      >  2:1248217/1‑76 (MQ=255)
                                      cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAg                                     <  2:1371954/76‑1 (MQ=255)
                                      cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAg                                     <  2:1494495/76‑1 (MQ=255)
                                      cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAg                                     <  2:1710628/76‑1 (MQ=255)
                                       ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGc                                    >  2:1967452/1‑76 (MQ=255)
                                        tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCg                                   <  2:1851291/76‑1 (MQ=255)
                                        tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCg                                   >  1:1592687/1‑76 (MQ=255)
                                           gATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgaga                                <  2:1448231/76‑1 (MQ=255)
                                           gATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgaga                                >  2:1518935/1‑76 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                               <  1:1638094/76‑1 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                               <  1:1638110/76‑1 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                               >  1:449151/1‑76 (MQ=255)
                                            aTGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGaa                               <  2:808919/76‑1 (MQ=255)
                                             tGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAAc                              <  2:193516/76‑1 (MQ=255)
                                             nggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAAc                              >  1:180863/2‑76 (MQ=255)
                                               gTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTg                            <  2:1110982/76‑1 (MQ=255)
                                               gTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTg                            >  1:490130/1‑76 (MQ=255)
                                                ttCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGc                           >  1:371461/1‑76 (MQ=255)
                                                 tCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGcc                          >  2:659901/1‑76 (MQ=255)
                                                   aaCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCtt                        >  1:2021027/1‑76 (MQ=255)
                                                      ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaaa                     <  1:78043/76‑1 (MQ=255)
                                                      ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaaa                     >  2:1585179/1‑76 (MQ=255)
                                                      ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaaa                     >  2:1153746/1‑76 (MQ=255)
                                                       cTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAc                    <  2:1012159/76‑1 (MQ=255)
                                                        ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACt                   >  1:1046560/1‑76 (MQ=255)
                                                        ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACt                   >  2:1953940/1‑76 (MQ=255)
                                                           aGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCa                >  1:407223/1‑76 (MQ=255)
                                                           aGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCa                >  1:486799/1‑76 (MQ=255)
                                                           aGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCa                >  1:650742/1‑76 (MQ=255)
                                                           aGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCa                >  2:877775/1‑76 (MQ=255)
                                                            gTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCAt               <  2:695570/76‑1 (MQ=255)
                                                             ttAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATa              <  2:1269958/76‑1 (MQ=255)
                                                             ttAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATa              >  1:1901226/1‑76 (MQ=255)
                                                               aGCGCTCATGGGTTTTACCCCCGCGTGGTGGCTACTTTTTCGATTTGTAAGCGAGAACTGCCTTAAACTCCATAcc            >  1:1184140/1‑76 (MQ=255)
                                                                gcgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCa           <  2:1675347/76‑1 (MQ=255)
                                                                gcgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCa           <  2:427262/76‑1 (MQ=255)
                                                                   cTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGt        >  1:1168527/1‑76 (MQ=255)
                                                                    ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt       >  2:977974/1‑76 (MQ=255)
                                                                    ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt       >  2:1028593/1‑76 (MQ=255)
                                                                     tATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtc      <  2:803860/76‑1 (MQ=255)
                                                                      aTGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtct     >  2:1564152/1‑76 (MQ=255)
                                                                      aTGGGTTTTACCCCCGAGTGGTGGCTACTTCTTCGTTTCGTAAGAGAGAACTGCCTTAAACTCCATACAACGTtct     >  2:1564166/1‑76 (MQ=255)
                                                                        gggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTctct   <  2:95722/76‑1 (MQ=255)
                                                                         ggTTTTACCGCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTg  <  2:1728468/76‑1 (MQ=255)
                                                                         ggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTg  >  2:1274779/1‑76 (MQ=255)
                                                                          |                                                                          
CGAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTG  >  REL606/2774122‑2774270

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