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)239 (1.070) 105/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)254 (1.140) 113/144 0.0 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

Only 100 of 256 total aligned reads displayed.

GCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772779‑2772924
                                                                           |                                                                      
gCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAcc                                                                        <  2:3321760/76‑1 (MQ=255)
 cGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAtcc                                                                       >  1:4603763/1‑76 (MQ=255)
 cGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAccc                                                                       >  1:7473710/1‑76 (MQ=255)
 cGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAccc                                                                       <  1:84301/76‑1 (MQ=255)
 cGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAccc                                                                       <  2:3098727/76‑1 (MQ=255)
  gAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAcccc                                                                      <  2:4764471/76‑1 (MQ=255)
     aTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATa                                                                   <  1:374691/76‑1 (MQ=255)
      ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                  >  1:5591374/1‑76 (MQ=255)
      ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                  <  1:1831944/76‑1 (MQ=255)
      ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                  >  1:3613727/1‑76 (MQ=255)
       tttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAg                                                                 >  1:6115442/1‑76 (MQ=255)
        ttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgc                                                                <  2:516119/76‑1 (MQ=255)
        ttGATTCATCGTGATGACCGGTTTATCCCCGCTCGCGCGGGGGACACCCTGCACGGGGGGAAAAAAACCCATAAgc                                                                >  1:2453508/1‑76 (MQ=25)
         tGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgcg                                                               >  2:4375239/1‑76 (MQ=255)
            ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                            >  2:803089/1‑76 (MQ=255)
             tCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGCACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                                                           >  1:5960779/1‑76 (MQ=255)
             tCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                                                           <  1:5400202/76‑1 (MQ=255)
             tCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCGGAACGGTGGGGAAAAACCCCATAAGCGCTaa                                                           <  2:5809666/76‑1 (MQ=255)
              cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGGAAAACCCCCTAAGCGCTAAc                                                          >  1:1267499/1‑76 (MQ=255)
              cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCCTAAGCGCTAAc                                                          >  1:1158381/1‑76 (MQ=255)
                tCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        <  2:464276/76‑1 (MQ=255)
                 cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                       <  1:670522/76‑1 (MQ=255)
                 cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                       <  1:5015672/76‑1 (MQ=255)
                 cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACCTa                                                       >  1:5185073/1‑76 (MQ=255)
                  gTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCGCATAAGCGCTAACTTag                                                      >  2:186949/1‑75 (MQ=255)
                   tgatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAg                                                     <  2:3620589/76‑1 (MQ=255)
                    gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    <  1:2560490/76‑1 (MQ=255)
                    gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    >  2:6655253/1‑76 (MQ=255)
                     atgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCCTAAGCGCTAACTTAAggg                                                   >  2:1963869/1‑76 (MQ=255)
                      tgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt                                                  <  1:110036/76‑1 (MQ=255)
                      tgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt                                                  >  2:1008015/1‑76 (MQ=255)
                       gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                 >  2:6873996/1‑76 (MQ=255)
                       gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                 <  2:4512011/76‑1 (MQ=255)
                        aCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtg                                                >  2:1241875/1‑76 (MQ=255)
                          cGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg                                              <  2:7418769/76‑1 (MQ=255)
                          cGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg                                              <  1:3367857/76‑1 (MQ=255)
                          cGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg                                              <  1:5993086/76‑1 (MQ=255)
                           ggTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                             <  1:6388794/76‑1 (MQ=255)
                            gTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                            <  2:960775/76‑1 (MQ=255)
                            gTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt                                            >  1:4026082/1‑76 (MQ=255)
                             tttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTa                                           <  2:595177/76‑1 (MQ=255)
                             tttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTATGGTa                                           <  1:1380856/76‑1 (MQ=255)
                              ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAGGGGTTGTGGTAt                                          >  1:7144983/1‑76 (MQ=255)
                              ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAATGGTTGGGGCAt                                          >  2:1754462/1‑76 (MQ=38)
                              ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          <  1:2264263/76‑1 (MQ=255)
                              ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          <  1:1960332/76‑1 (MQ=255)
                              ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  2:6826653/1‑76 (MQ=255)
                               tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         >  2:4172669/1‑76 (MQ=255)
                               tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:2762856/76‑1 (MQ=255)
                               tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  2:4896109/76‑1 (MQ=255)
                               tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:5629734/76‑1 (MQ=255)
                                aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  1:7380952/76‑1 (MQ=255)
                                aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  2:170452/76‑1 (MQ=255)
                                 tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  1:3234985/76‑1 (MQ=255)
                                 tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  2:2302680/76‑1 (MQ=255)
                                 accccGCTCGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  1:1694944/75‑1 (MQ=255)
                                  ccccGCTGGCGCGGGGAACACCCTGAACGGTTGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      >  1:1473574/1‑76 (MQ=255)
                                  ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTaag                                      >  1:6927676/1‑74 (MQ=255)
                                  ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAct                                      <  1:2335154/76‑2 (MQ=255)
                                  ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      >  1:3092721/1‑76 (MQ=255)
                                  ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      <  2:5642622/76‑1 (MQ=255)
                                  ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      <  2:7473398/76‑1 (MQ=255)
                                   cccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                     >  2:2592617/1‑76 (MQ=255)
                                      gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  <  1:4227938/76‑1 (MQ=255)
                                       cTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGa                                 >  2:1552057/1‑76 (MQ=255)
                                        tGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                >  1:6231011/1‑76 (MQ=255)
                                         ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGatt                               >  1:5655332/1‑75 (MQ=255)
                                         ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                               <  1:2946516/76‑1 (MQ=255)
                                         ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                               <  1:4981248/76‑1 (MQ=255)
                                         ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                               >  2:7317819/1‑76 (MQ=255)
                                           cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             >  2:6985521/1‑76 (MQ=255)
                                            gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                            <  1:2678653/76‑1 (MQ=255)
                                            gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGa                            <  1:3183337/76‑1 (MQ=255)
                                             cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                           >  2:7130028/1‑76 (MQ=255)
                                              ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          <  2:6333827/76‑1 (MQ=255)
                                               gggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAttt                         >  1:6015345/1‑76 (MQ=255)
                                                ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        >  2:2455123/1‑76 (MQ=255)
                                                 gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       >  2:5443406/1‑76 (MQ=255)
                                                 gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       >  2:1918050/1‑76 (MQ=255)
                                                     acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  2:3133191/76‑1 (MQ=255)
                                                     acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  2:2572188/76‑1 (MQ=255)
                                                     acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  2:2238825/76‑1 (MQ=255)
                                                     acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTg                   <  2:2102680/76‑1 (MQ=255)
                                                      cccTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGc                  >  1:2399322/1‑76 (MQ=255)
                                                        cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                >  1:7017491/1‑76 (MQ=255)
                                                          gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt              >  1:6603362/1‑76 (MQ=255)
                                                          gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt              <  2:3748933/76‑1 (MQ=255)
                                                           aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg             <  1:934138/76‑1 (MQ=255)
                                                             cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa           >  2:658911/1‑76 (MQ=255)
                                                                tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        <  1:474038/76‑1 (MQ=40)
                                                                 ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc       >  2:121259/1‑76 (MQ=39)
                                                                 ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc       <  1:1000815/76‑1 (MQ=39)
                                                                  gggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct      >  1:87472/1‑76 (MQ=38)
                                                                  gggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct      <  1:4232365/76‑1 (MQ=38)
                                                                   ggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActc     >  1:6054406/1‑76 (MQ=35)
                                                                    gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  1:3716110/76‑1 (MQ=34)
                                                                    gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    <  1:3469523/76‑1 (MQ=34)
                                                                    gAAAAACACCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTActct    >  2:6176704/1‑76 (MQ=21)
                                                                     aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctc   <  1:3881488/76‑1 (MQ=32)
                                                                      aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  >  2:6787312/1‑76 (MQ=31)
                                                                           |                                                                      
GCGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772779‑2772924

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