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)48 (1.210) 40/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)37 (0.940) 29/144 0.3 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

TTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTC  >  REL606/2772785‑2772923
                                                                     |                                                                     
ttttGATTCATCGTGATAACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                 <  1:145316/76‑1 (MQ=255)
 tttGATTCATCGTGATGACCGGTTTATCCCCGCCGGCGCGGGGGACACCCCGGACGGGGGGGGAAAAACCCCTAAg                                                                >  1:990245/1‑76 (MQ=21)
        cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                         <  2:315421/76‑1 (MQ=255)
         aTCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACt                                                        <  2:973363/76‑1 (MQ=255)
           cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                      <  2:320651/76‑1 (MQ=255)
           cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                      >  2:1309967/1‑76 (MQ=255)
                tgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt                                                 <  2:1030498/76‑1 (MQ=255)
                 gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                >  2:1001624/1‑76 (MQ=255)
                 gaCCCGATTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAACAACCCCATAAGCGCTAACTTAAGGGtt                                                <  1:1302311/76‑1 (MQ=38)
                       tttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                          >  2:622733/1‑76 (MQ=255)
                       tttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                          <  1:757575/76‑1 (MQ=255)
                        ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                         >  2:81959/1‑76 (MQ=255)
                        ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGNTGTGGTAt                                         >  1:1155433/1‑76 (MQ=255)
                             cccGGTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                    <  2:555778/76‑1 (MQ=255)
                               cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                  <  1:710657/76‑1 (MQ=255)
                                  gggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                               <  1:881428/75‑1 (MQ=255)
                                   ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                              <  1:923560/76‑1 (MQ=255)
                                       tggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                          >  2:1073440/2‑76 (MQ=255)
                                         gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                        <  2:69288/76‑1 (MQ=255)
                                          ggAACACCCTGAACGGTGGGGAAAATCCNCATAAGCGCACACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                       <  1:112874/76‑1 (MQ=18)
                                          ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                       <  1:834267/76‑1 (MQ=255)
                                          ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                       <  1:996382/76‑1 (MQ=255)
                                          ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGGTGTGGTATTACGCCTGATATGATTTa                       <  1:834269/76‑1 (MQ=255)
                                             acacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                    <  2:747121/76‑1 (MQ=255)
                                             acacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                    >  2:133224/1‑76 (MQ=255)
                                             acacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                    >  2:499708/1‑76 (MQ=255)
                                               acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTACCGTg                  <  2:1299185/76‑1 (MQ=255)
                                               acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                  <  2:129979/76‑1 (MQ=255)
                                                 ccTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGcc                >  2:688575/1‑76 (MQ=255)
                                                  cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg               <  2:1028883/76‑1 (MQ=255)
                                                  cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg               <  2:1042304/76‑1 (MQ=255)
                                                  cTGAACCGTTGGGGAAAACCCCATAAGCCCTAGCATAACGGTTGTCGTATTACGCCTGATATGATTTAACGTGCCg               <  2:1028896/76‑1 (MQ=11)
                                                   tGAACGGTGGTGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGa              >  1:1218333/1‑76 (MQ=255)
                                                   tGAACGGTGGGGAAAAACCCCATNAGCGCNAACTTCAGGGTTGTNGTATTAAGCCTGATATGATTTAACGTGCCGa              >  2:172515/1‑76 (MQ=17)
                                                    gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt             >  1:268063/1‑76 (MQ=255)
                                                    gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt             >  1:672126/1‑76 (MQ=255)
                                                     aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg            <  1:460478/76‑1 (MQ=255)
                                                     aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg            >  2:1026302/1‑76 (MQ=255)
                                                       cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa          <  1:576093/76‑1 (MQ=255)
                                                               aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc  >  1:547854/1‑76 (MQ=32)
                                                               aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc  <  2:1030698/76‑1 (MQ=32)
                                                                     |                                                                     
TTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTC  >  REL606/2772785‑2772923

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