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)65 (1.150) 43/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)50 (0.890) 41/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)

Only 100 of 125 total aligned reads displayed.

CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524447‑4524596
                                                                           |                                                                          
cgaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggc                                                                            >  2:412949‑M2/76‑76 (MQ=255)
cATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGc                                                                            >  1:286796/1‑76 (MQ=255)
 gaggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggcc                                                                           >  1:501822‑M2/75‑76 (MQ=255)
  tttAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCg                                                                          >  1:445571/1‑76 (MQ=255)
  aggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCg                                                                          <  2:1003619‑M2/3‑1 (MQ=255)
   ggcatcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGc                                                                         <  1:1478383‑M2/4‑1 (MQ=255)
       ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                     >  2:1158095‑M2/1‑74 (MQ=255)
       ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                     >  1:450880‑M2/1‑74 (MQ=255)
       ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                     <  2:1697537‑M2/76‑3 (MQ=255)
       tcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcc                                                                     <  1:1566612‑M2/8‑1 (MQ=255)
        tCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                    <  1:141693‑M2/76‑4 (MQ=255)
        cgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAccc                                                                    >  2:235654‑M2/68‑76 (MQ=255)
         gtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAcccc                                                                   >  1:770901‑M2/67‑76 (MQ=255)
          aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                  >  1:242039‑M2/1‑71 (MQ=255)
          aGAAACAACCCTGAATAACGGAACCAATACCATCCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                  >  1:603875‑M2/1‑71 (MQ=255)
           cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                 >  1:215817‑M2/65‑76 (MQ=255)
           cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                 <  1:673136‑M2/12‑1 (MQ=255)
            aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                >  1:770591‑M2/1‑69 (MQ=255)
            aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                >  1:1005908‑M2/1‑69 (MQ=255)
             aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                               >  1:1257594‑M2/1‑68 (MQ=255)
              acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                              >  2:303070‑M2/1‑67 (MQ=255)
              acaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagc                                                              >  2:526480‑M2/1‑67 (MQ=255)
              aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                              >  2:1291713‑M2/62‑76 (MQ=255)
                cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                            >  2:1045836‑M2/60‑76 (MQ=255)
                 acCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgct                                                           <  2:807434‑M2/76‑13 (MQ=255)
                     tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                       >  2:700133‑M2/1‑60 (MQ=255)
                     tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                       >  2:1803833‑M2/1‑60 (MQ=255)
                     caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                       <  1:344704‑M2/22‑1 (MQ=255)
                      gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                      >  1:14582‑M2/1‑59 (MQ=255)
                      gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                      >  1:1627099‑M2/1‑59 (MQ=255)
                       ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                     <  1:51435‑M2/24‑1 (MQ=255)
                       aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                     <  2:682446‑M2/76‑19 (MQ=255)
                        cgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGc                                                    <  1:964574‑M2/25‑1 (MQ=255)
                        aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                    >  1:1080213‑M2/1‑57 (MQ=255)
                         tAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaag                                                   >  1:438970‑M2/1‑56 (MQ=255)
                          aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                  <  1:554227‑M2/76‑22 (MQ=255)
                          aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                  <  2:1489326‑M2/76‑22 (MQ=255)
                          aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                  <  2:617127‑M2/76‑22 (MQ=255)
                          aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                  <  1:1205102‑M2/76‑22 (MQ=255)
                          aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                  <  1:154343‑M2/76‑22 (MQ=255)
                          aaCGGAACCAATACCATTCCGTTCCAGGCGCGGTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                  <  1:1597947‑M2/76‑22 (MQ=255)
                          aaCGGAACCAATACCATCCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAActcataagcgctaacttaagg                                                  >  1:911960‑M2/1‑55 (MQ=255)
                            cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                <  2:51871‑M2/76‑24 (MQ=255)
                              gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                              <  1:213932‑M2/31‑1 (MQ=255)
                              gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                              <  1:1679738‑M2/31‑1 (MQ=255)
                                cgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGa                                            >  2:204264‑M2/44‑76 (MQ=255)
                                   attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                         >  2:92230‑M2/41‑76 (MQ=255)
                                   aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggta                                         >  1:796707‑M2/1‑46 (MQ=255)
                                      aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                      >  2:634570‑M2/1‑43 (MQ=255)
                                      aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                      >  1:709442‑M2/1‑43 (MQ=255)
                                       ttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCaa                                     <  1:953587‑M2/40‑1 (MQ=255)
                                        tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                    >  1:954475‑M2/36‑76 (MQ=255)
                                         aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                   <  1:1287741‑M2/76‑37 (MQ=255)
                                           tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                 >  2:121382‑M2/1‑38 (MQ=255)
                                           tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtgggattacgcct                                 <  1:1185800‑M2/76‑39 (MQ=255)
                                           gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                 <  1:1230911‑M2/44‑1 (MQ=255)
                                            atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                <  2:1481512‑M2/45‑1 (MQ=255)
                                             tggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                               >  1:1762721‑M2/31‑76 (MQ=255)
                                               ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgata                             <  1:1750553‑M2/76‑43 (MQ=255)
                                               gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                             <  1:449245‑M2/48‑1 (MQ=255)
                                               gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                             <  1:124436‑M2/48‑1 (MQ=255)
                                               gttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                             <  1:1219893‑M2/48‑1 (MQ=255)
                                                tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                            <  1:1402967‑M2/76‑44 (MQ=255)
                                                tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                            >  2:163391‑M2/1‑33 (MQ=255)
                                                 tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                           <  2:1676888‑M2/50‑1 (MQ=255)
                                                 tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                           >  2:858317‑M2/27‑76 (MQ=255)
                                                 tcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTAt                           <  2:446017‑M2/50‑1 (MQ=255)
                                                 ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                           <  2:318690‑M2/76‑45 (MQ=255)
                                                 ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                           >  1:1733268‑M2/1‑32 (MQ=255)
                                                   aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                         <  1:682851‑M2/76‑47 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                      >  1:1001926‑M2/1‑27 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                      <  2:1587862‑M2/76‑50 (MQ=255)
                                                      cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                      >  1:1525332‑M2/1‑27 (MQ=255)
                                                       gcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaa                     >  1:752613‑M2/1‑26 (MQ=255)
                                                         taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                   <  1:579626‑M2/58‑1 (MQ=255)
                                                          nagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                  >  1:544602‑M2/18‑76 (MQ=255)
                                                          aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                  >  2:1226858‑M2/18‑76 (MQ=255)
                                                           tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                 >  2:1801223‑M2/1‑22 (MQ=255)
                                                           agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                 <  1:528049‑M2/60‑1 (MQ=255)
                                                            aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                <  1:52166‑M2/76‑56 (MQ=255)
                                                              tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg              >  2:1049645‑M2/1‑19 (MQ=255)
                                                              tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg              <  1:1605127‑M2/76‑58 (MQ=255)
                                                              tttGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccg              >  1:1571429‑M2/1‑19 (MQ=255)
                                                               agcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCa             <  2:1309319‑M2/64‑1 (MQ=255)
                                                                tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgat            <  1:1790889‑M2/76‑60 (MQ=255)
                                                                tGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgan            >  1:1003618‑M2/1‑17 (MQ=255)
                                                                 gCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatg           >  2:1752823‑M2/1‑16 (MQ=255)
                                                                   cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt         >  1:18873‑M2/9‑76 (MQ=255)
                                                                     tatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCa       <  1:683092‑M2/70‑1 (MQ=255)
                                                                      cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta      >  2:608171‑M2/1‑11 (MQ=255)
                                                                      cGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaatta      >  2:1604242‑M2/1‑11 (MQ=255)
                                                                       tgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg     >  2:1781061‑M2/5‑76 (MQ=255)
                                                                        gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg    <  1:381821‑M2/73‑1 (MQ=255)
                                                                        gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg    <  1:371092‑M2/73‑1 (MQ=255)
                                                                         ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc   <  2:723692‑M2/76‑69 (MQ=255)
                                                                         ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc   <  2:886138‑M2/76‑69 (MQ=255)
                                                                         ggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactc   <  1:1259393‑M2/76‑69 (MQ=255)
                                                                         ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa   <  2:1374709‑M2/74‑1 (MQ=255)
                                                                         ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa   <  2:7597‑M2/74‑1 (MQ=255)
                                                                          gCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactct  >  2:742150‑M2/1‑7 (MQ=255)
                                                                           |                                                                          
CATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGAC  >  REL606/4524447‑4524596

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.