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)63 (1.370) 47/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)47 (1.030) 41/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 113 total aligned reads displayed.

GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524453‑4524602
                                                                          |                                                                           
atcgtcttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAAc                                                                            >  1:822553‑M2/70‑76 (MQ=255)
 ttCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                           <  2:2464‑M2/76‑3 (MQ=255)
    aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccat                                                                        >  2:808832‑M2/1‑71 (MQ=255)
     cttaacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCgg                                                                       >  2:1169673‑M2/65‑76 (MQ=255)
      aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccatan                                                                      >  1:295482‑M2/1‑69 (MQ=255)
      aaaCAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataa                                                                      >  2:124390‑M2/1‑69 (MQ=255)
       aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                     >  1:37228‑M2/1‑68 (MQ=255)
       aacaacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataag                                                                     >  2:1360743‑M2/1‑68 (MQ=255)
        aacgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtg                                                                    <  2:166837‑M2/15‑1 (MQ=255)
          cgaggcaccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGtgct                                                                  >  2:1198188‑M2/60‑76 (MQ=255)
            cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgcta                                                                >  1:268595‑M2/1‑63 (MQ=255)
             ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                               >  1:1033722‑M2/1‑62 (MQ=255)
             ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                               >  2:1242317‑M2/1‑62 (MQ=255)
             ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                               <  2:1289273‑M2/76‑15 (MQ=255)
             ccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaa                                                               >  1:1055520‑M2/1‑62 (MQ=255)
               tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             <  1:1403945‑M2/76‑17 (MQ=255)
               tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaact                                                             <  2:849834‑M2/76‑17 (MQ=255)
               caccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTaa                                                             >  1:280965‑M2/55‑76 (MQ=255)
                gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactt                                                            <  1:974809‑M2/76‑18 (MQ=255)
                gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactn                                                            >  1:933730‑M2/1‑59 (MQ=255)
                 ccgaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATg                                                           <  1:1280118‑M2/24‑1 (MQ=255)
                 aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           <  2:1271427‑M2/76‑19 (MQ=255)
                 aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           >  2:1211783‑M2/1‑58 (MQ=255)
                 aaTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaactta                                                           >  2:436073‑M2/1‑58 (MQ=255)
                  aTAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaa                                                          <  2:390315‑M2/76‑20 (MQ=255)
                   gaggcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCg                                                         >  2:585544‑M2/51‑76 (MQ=255)
                   gaggcgccgcattcttcatacggttcaactcttaagttagctcttatggtCCGCAACCCCGGGTGCTGCTAATGCg                                                         <  2:604684‑M2/26‑1 (MQ=255)
                    aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                        <  2:741731‑M2/76‑22 (MQ=255)
                    aaCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagg                                                        <  2:691032‑M2/76‑22 (MQ=255)
                      gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                      <  1:534607‑M2/29‑1 (MQ=255)
                      gcgtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                      <  1:367699‑M2/29‑1 (MQ=255)
                      cGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggt                                                      >  1:1179056‑M2/1‑53 (MQ=255)
                       ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     <  2:826108‑M2/76‑25 (MQ=255)
                       ggAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggtt                                                     >  2:680350‑M2/1‑52 (MQ=255)
                        gtcgcattcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGc                                                    <  2:1503443‑M2/31‑1 (MQ=255)
                        gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                                    >  1:384907‑M2/1‑51 (MQ=255)
                        gAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttg                                                    >  1:1020641‑M2/1‑51 (MQ=255)
                         aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                                   <  2:352312‑M2/76‑27 (MQ=255)
                         aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgt                                                   >  2:1214664‑M2/1‑50 (MQ=255)
                          aCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtg                                                  >  1:1492495‑M2/1‑49 (MQ=255)
                             attcttcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt                                               >  1:355421‑M2/41‑76 (MQ=255)
                               tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                             <  2:1258905‑M2/76‑33 (MQ=255)
                               tACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatt                                             >  2:1033537‑M2/1‑44 (MQ=255)
                                ctgcagatggttcaacccgtaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa                                            <  2:79222‑M2/39‑1 (MQ=255)
                                aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtatta                                            <  2:1407232‑M2/76‑34 (MQ=255)
                                 ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  2:302862‑M2/1‑42 (MQ=255)
                                 ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattac                                           >  1:305812‑M2/1‑42 (MQ=255)
                                  tcagatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg                                          >  2:1240453‑M2/36‑76 (MQ=255)
                                  cATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacg                                          >  1:826108‑M2/1‑41 (MQ=255)
                                   aTTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgc                                         <  1:1089903‑M2/76‑37 (MQ=255)
                                     tCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcct                                       <  1:20240‑M2/76‑39 (MQ=255)
                                     gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       <  2:1430984‑M2/44‑1 (MQ=255)
                                     gatggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGtt                                       <  1:897845‑M2/44‑1 (MQ=255)
                                      ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      >  2:1345761‑M2/1‑37 (MQ=255)
                                      ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctg                                      >  1:701115‑M2/1‑37 (MQ=255)
                                      atggttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                      >  1:397586‑M2/32‑76 (MQ=255)
                                          ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                                  >  1:1515580‑M2/28‑76 (MQ=255)
                                          ttcaacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTa                                  >  2:1380013‑M2/28‑76 (MQ=255)
                                          tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatat                                  >  1:647616‑M2/1‑33 (MQ=255)
                                          cccaGGCGCTTTCTTCTGCAGCCGGGGCCGCAAcccataagcgctaccttaagggttgtggtattacgcctgatat                                  >  2:1028547‑M2/2‑33 (MQ=255)
                                           ccAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatg                                 <  1:990472‑M2/76‑45 (MQ=255)
                                            caacccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                                <  1:904089‑M2/51‑1 (MQ=255)
                                             aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                               <  1:1082095‑M2/76‑47 (MQ=255)
                                             aGGCGCGTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgat                               >  1:78645‑M2/1‑30 (MQ=255)
                                              acccttaagtgagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCaa                              <  1:1375497‑M2/53‑1 (MQ=255)
                                                cgcgTTATTTTGCAACCGGGGCCGTAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                            >  2:13619‑M2/1‑27 (MQ=255)
                                                cgcgTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgattta                            >  2:177348‑M2/1‑27 (MQ=255)
                                                ccttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATa                            <  2:660530‑M2/55‑1 (MQ=255)
                                                  ttaagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                          >  2:475491‑M2/20‑76 (MQ=255)
                                                   taagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAcc                         <  1:209147‑M2/58‑1 (MQ=255)
                                                   gTTATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacg                         >  1:1095485‑M2/1‑24 (MQ=255)
                                                    ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                        >  1:1124868‑M2/1‑23 (MQ=255)
                                                    ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                        >  2:536245‑M2/1‑23 (MQ=255)
                                                    ttATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgt                        >  2:345975‑M2/1‑23 (MQ=255)
                                                    aagttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCt                        >  2:483804‑M2/18‑76 (MQ=255)
                                                     tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       <  1:436185‑M2/76‑55 (MQ=255)
                                                     tATTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtg                       >  2:647471‑M2/1‑22 (MQ=255)
                                                     agttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTa                       >  1:184541‑M2/17‑76 (MQ=255)
                                                      gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                      >  2:162066‑M2/16‑76 (MQ=255)
                                                      gttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAc                      <  1:939150‑M2/61‑1 (MQ=255)
                                                      aTTTTGCAACCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgc                      <  2:1380721‑M2/76‑56 (MQ=255)
                                                       ttagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAcc                     >  1:1240270‑M2/15‑76 (MQ=255)
                                                        tagcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc                    >  1:73775‑M2/14‑76 (MQ=255)
                                                          gcgcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAg                  <  2:859830‑M2/65‑1 (MQ=255)
                                                            gcttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt                >  1:107018‑M2/10‑76 (MQ=255)
                                                            cAACCGGGGCCGCAAcccataagcgctaccttaagggttgtgggattacgcctgatatgatttaacgtgccgatga                <  1:1397969‑M2/76‑62 (MQ=255)
                                                             cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGtt               <  2:1266547‑M2/68‑1 (MQ=255)
                                                             aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaa               >  1:908915‑M2/1‑14 (MQ=255)
                                                             aaCCGGGGCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgggccgatgaa               <  2:1146364‑M2/76‑63 (MQ=255)
                                                              ttatgggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTc              <  2:1298064‑M2/69‑1 (MQ=255)
                                                                 ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac           <  1:698945‑M2/76‑67 (MQ=255)
                                                                 ggggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattac           >  1:67802‑M2/1‑10 (MQ=255)
                                                                  gggCCGCAAcccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattact          >  1:925353‑M2/1‑9 (MQ=255)
                                                                  gggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg          <  2:104753‑M2/73‑1 (MQ=255)
                                                                   ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa         <  2:1347274‑M2/74‑1 (MQ=255)
                                                                   ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGa         <  2:619229‑M2/74‑1 (MQ=255)
                                                                     ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACg       >  2:1196495/1‑76 (MQ=255)
                                                                      cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGt      <  2:42383/76‑1 (MQ=255)
                                                                       gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTc     <  1:288578/76‑1 (MQ=255)
                                                                          acccataagcgctaacttaagggttgtggtattacgcctgatatgatttaacgtgccgatgaattactctcacgat  >  2:838247‑M2/1‑1 (MQ=255)
                                                                          |                                                                           
GTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATT  >  REL606/4524453‑4524602

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.