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)296 (1.260) 117/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)249 (1.060) 107/144 0.1 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 284 total aligned reads displayed.

CGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772780‑2772924
                                                                          |                                                                      
   aaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                                                                    >  2:6655896/1‑76 (MQ=255)
   aaTTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGAAACACCCTGAACGGTGGGGAAAAACCCCAt                                                                    >  2:1991967/1‑76 (MQ=255)
     ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGGACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                  >  2:4155838/1‑76 (MQ=255)
     ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                  <  2:781059/76‑1 (MQ=255)
     ttttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATaa                                                                  <  1:2831747/76‑1 (MQ=255)
       ttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgc                                                                <  1:2936450/76‑1 (MQ=255)
       ttGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAACCCATAAgc                                                                >  1:4143863/1‑76 (MQ=255)
        tGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAgcg                                                               <  2:1032117/76‑1 (MQ=255)
           ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                            >  2:6250833/1‑76 (MQ=255)
           ttCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                                                            <  1:920653/76‑1 (MQ=255)
             cATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAAc                                                          >  1:361771/1‑76 (MQ=255)
              aTCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACt                                                         >  1:2089373/1‑76 (MQ=255)
              aTCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACt                                                         <  1:7632448/76‑1 (MQ=255)
               tCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                                                        >  1:3475943/1‑76 (MQ=255)
                cGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa                                                       >  2:748956/1‑76 (MQ=255)
                 gTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                      <  1:6740344/76‑1 (MQ=255)
                 gTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                      >  2:6932066/1‑76 (MQ=255)
                 gTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTaa                                                      >  1:1100088/1‑76 (MQ=255)
                   gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    >  1:1031982/1‑76 (MQ=255)
                   gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    <  1:7201186/76‑1 (MQ=255)
                   gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    <  1:5900223/76‑1 (MQ=255)
                   gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    <  1:4870307/76‑1 (MQ=255)
                   gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    <  2:7156889/76‑1 (MQ=255)
                   gatgaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg                                                    <  1:1008380/76‑1 (MQ=255)
                      gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                 >  1:5771812/1‑76 (MQ=255)
                      gaCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGtt                                                 <  1:4061034/76‑1 (MQ=255)
                       aCCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtg                                                <  2:5145361/76‑1 (MQ=255)
                        ccGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGCAAAAACGCCATAAGCGCTAACTTAAGGGgtgt                                               >  2:1394818/1‑76 (MQ=37)
                          ggTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                             >  2:2803926/1‑76 (MQ=255)
                          ggTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTgg                                             >  2:1587920/1‑76 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          <  1:3731724/76‑1 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  1:2713958/1‑76 (MQ=255)
                             ttATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAt                                          >  2:7037214/1‑76 (MQ=255)
                              tATCCTCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:310077/76‑1 (MQ=255)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAATCCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         <  1:3930538/76‑1 (MQ=255)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         >  2:2907878/1‑76 (MQ=255)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         >  2:5681069/1‑76 (MQ=255)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         >  2:3969581/1‑76 (MQ=255)
                              tATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt                                         >  2:1126163/1‑76 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAGGCGCTAACTTAAGGGTTGTGGTATTa                                        <  1:5154073/76‑1 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTCACTTAAGGGTTGTGGTATTa                                        >  1:687876/1‑76 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTGAGGGTTGTGGTATTa                                        >  1:2014123/1‑76 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAttc                                        >  1:3253359/1‑75 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        <  1:2324659/76‑1 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        >  2:2352445/1‑76 (MQ=255)
                               aTCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAAACCCATAAGCGCTAACTTAAGGGTTGTGGTATTa                                        >  1:6017062/1‑76 (MQ=255)
                                tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       >  2:7271103/1‑76 (MQ=255)
                                tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  1:1140534/76‑1 (MQ=255)
                                tCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc                                       <  1:2491008/76‑1 (MQ=255)
                                 ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGGAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      >  2:3927875/1‑76 (MQ=255)
                                 ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      >  1:1668443/1‑76 (MQ=255)
                                 ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      <  1:5412162/76‑1 (MQ=255)
                                 ccccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg                                      <  1:1696854/76‑1 (MQ=255)
                                  cccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc                                     <  2:6888120/76‑1 (MQ=255)
                                   ccGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGcc                                    >  2:5882892/1‑76 (MQ=255)
                                    cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                   >  1:455316/1‑76 (MQ=255)
                                    cGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCt                                   >  2:2135184/1‑76 (MQ=255)
                                     gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  <  2:6531902/76‑1 (MQ=255)
                                     gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  <  2:6434953/76‑1 (MQ=255)
                                     gCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTg                                  <  1:480799/76‑1 (MQ=255)
                                       tGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGat                                >  2:502497/1‑76 (MQ=255)
                                        ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                               <  1:5086904/76‑1 (MQ=255)
                                        ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                               >  1:2617650/1‑76 (MQ=255)
                                        ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                               >  1:3323262/1‑76 (MQ=255)
                                        ggCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGata                               >  1:7420950/1‑76 (MQ=255)
                                          cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             >  2:437595/1‑76 (MQ=255)
                                          cgcgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATg                             <  2:3286937/76‑1 (MQ=255)
                                            cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAt                           >  1:5665348/1‑76 (MQ=255)
                                             ggggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGAtt                          >  1:2441249/1‑76 (MQ=255)
                                               ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        >  2:3741781/1‑76 (MQ=255)
                                               ggAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTa                        <  1:4676129/76‑1 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       <  1:344126/76‑1 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       <  2:4239035/76‑1 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       <  1:477442/76‑1 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       >  1:2210296/1‑76 (MQ=255)
                                                gAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTaa                       >  1:5531431/1‑76 (MQ=255)
                                                  acacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACg                     <  2:2669515/76‑1 (MQ=255)
                                                   cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGt                    >  2:152745/1‑76 (MQ=255)
                                                       cTGAACGGTGGGGAAAAGCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                <  1:7335449/76‑1 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                >  1:6463727/1‑76 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                <  2:6786960/76‑1 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                >  1:2503921/1‑76 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                >  2:5566597/1‑76 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                <  1:701090/76‑1 (MQ=255)
                                                       cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCg                >  1:1431843/1‑76 (MQ=255)
                                                         gAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGAt              >  2:4766431/1‑76 (MQ=255)
                                                          aaCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATg             >  1:427989/1‑76 (MQ=255)
                                                           aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa            <  2:4426704/76‑1 (MQ=255)
                                                           aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGa            <  1:176753/76‑1 (MQ=255)
                                                            cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGaa           >  2:2205323/1‑76 (MQ=255)
                                                             ggTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAt          <  1:5212684/76‑1 (MQ=255)
                                                              gTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAAtt         <  2:5964690/76‑1 (MQ=255)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        >  1:1874022/1‑76 (MQ=40)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTa        >  1:4045054/1‑76 (MQ=40)
                                                               tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATCTAACGTGCCGATGAATTa        >  1:1199867/1‑76 (MQ=255)
                                                                ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc       >  2:5171163/1‑76 (MQ=39)
                                                                ggggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAc       >  2:4911143/1‑76 (MQ=39)
                                                                 gggAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTAct      >  1:7837357/1‑76 (MQ=38)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACtctcc  >  1:7926/1‑75 (MQ=14)
                                                                     aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCa  <  1:6157646/76‑1 (MQ=31)
                                                                          |                                                                      
CGAAATTTTGATTCATCGTGATGACCGGTTTATCCCCGCTGGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGCCTGATATGATTTAACGTGCCGATGAATTACTCTCA  >  REL606/2772780‑2772924

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