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)43 (1.100) 35/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)34 (0.870) 29/144 0.4 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

GGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTGA  >  REL606/2774125‑2774271
                                                                       |                                                                           
ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                         <  1:411821/76‑1 (MQ=35)
   aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                      >  2:545426/1‑76 (MQ=40)
     cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCATCCCTTATGTTAGCGCTTNTGGGTTTTAcccc                                                                    >  1:161871/1‑76 (MQ=255)
           aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtg                                                              <  2:1294585/76‑1 (MQ=255)
              gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                           <  1:518984/76‑1 (MQ=255)
               aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGtggtgg                                                          <  2:398261/76‑1 (MQ=255)
               aGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGtggtgg                                                          >  1:720897/1‑76 (MQ=255)
                 gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCt                                                        <  2:535268/76‑1 (MQ=255)
                 gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCt                                                        <  2:260092/76‑1 (MQ=255)
                   aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTAc                                                      >  1:728700/1‑76 (MQ=255)
                        ggCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtctt                                                 <  2:1007755/76‑1 (MQ=255)
                         gCGTCGCATTCTTNAGATGGTTCAACCCNTANGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACttcttc                                                <  2:1166155/76‑1 (MQ=255)
                           gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGa                                              >  1:876396/1‑76 (MQ=255)
                             cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAtt                                            >  1:1073347/1‑76 (MQ=255)
                              gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAtta                                           <  2:1296591/76‑2 (MQ=255)
                              gCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTc                                           >  1:465279/1‑76 (MQ=255)
                               cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGNTGGCTACTTCTTCGATTCg                                          <  1:1210995/76‑1 (MQ=255)
                               cATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCg                                          >  1:650881/1‑76 (MQ=255)
                                aTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGt                                         >  1:506562/1‑76 (MQ=255)
                                 ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                        <  1:541891/76‑1 (MQ=255)
                                     tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCg                                    <  1:291678/76‑1 (MQ=255)
                                           ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                              >  2:847020/1‑76 (MQ=255)
                                           ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                              <  2:461130/76‑1 (MQ=255)
                                              tCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGcc                           >  2:660238/1‑76 (MQ=255)
                                               cAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCt                          >  2:331133/1‑76 (MQ=255)
                                                aaCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTTGCACTCGAGAACTGCttt                         >  2:105431/1‑76 (MQ=25)
                                                aaCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCtt                         >  2:105420/1‑76 (MQ=255)
                                                  cccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTaa                       >  1:700623/1‑76 (MQ=255)
                                                       aaGTTAGCGCTTATGGGTTTTACCCCCGCGTGTTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTcc                  <  2:316912/76‑1 (MQ=255)
                                                       aaGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGCTTCGTCAGCGAGAACTGCCTTAAACTcc                  >  1:586683/1‑76 (MQ=255)
                                                       aaGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTcc                  <  1:947412/76‑1 (MQ=255)
                                                       aaGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCCTAAACGAGGACTGCCTTAAAATcc                  >  1:864956/1‑76 (MQ=37)
                                                            aGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAcc             <  1:604468/76‑1 (MQ=255)
                                                            aGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCATAAACTCCATAcc             >  1:1110777/1‑76 (MQ=255)
                                                                 ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCNTACCACGtt        <  1:24194/76‑1 (MQ=255)
                                                                 ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt        >  2:863820/1‑76 (MQ=255)
                                                                       gTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTTAACTCCATACCACGTTCTCTGa  <  1:9544/76‑1 (MQ=255)
                                                                       gTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTGa  <  2:977471/76‑1 (MQ=255)
                                                                       |                                                                           
GGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTGA  >  REL606/2774125‑2774271

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.