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.120) 39/144 0.1 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)46 (1.080) 39/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)

GAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTG  >  REL606/2774123‑2774270
                                                                         |                                                                          
gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtt                                                                          <  2:1369919/76‑1 (MQ=32)
 aGGCATCGTCTTANCGAGGCACCGAGGCNTCNGATTCTTCAGATGGTTGAACCCTTAAGTTAGCGCTTATGGGttt                                                                         <  2:1279342/76‑1 (MQ=11)
   gCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTa                                                                       <  1:1092863/76‑1 (MQ=38)
   gCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTa                                                                       <  2:1044970/76‑1 (MQ=38)
    cATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAc                                                                      >  2:38361/1‑76 (MQ=39)
     aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                     >  2:782741/1‑76 (MQ=40)
         tCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcg                                                                 >  2:1157520/1‑76 (MQ=255)
             aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtg                                                             <  2:632533/76‑1 (MQ=255)
              aCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtgg                                                            <  1:334174/76‑1 (MQ=255)
               cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggt                                                           >  1:333991/1‑76 (MQ=255)
               cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggt                                                           <  2:327930/76‑1 (MQ=255)
                gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAGCCCTTAAGTTAGCGCTTATTGGTTTTACCCCCGCgtggtg                                                          <  2:1164258/76‑1 (MQ=255)
                gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                          <  2:291663/76‑1 (MQ=255)
                  ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGc                                                        <  1:427018/76‑1 (MQ=255)
                   gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCt                                                       >  2:1164997/1‑76 (MQ=255)
                   gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCt                                                       >  2:1389124/1‑76 (MQ=255)
                     aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGGGGCTAc                                                     >  1:374761/1‑76 (MQ=255)
                      ccGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTAct                                                    >  2:449197/1‑76 (MQ=255)
                       cGAGGCGTCGCATTCTTCANATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActt                                                   >  1:1321745/1‑76 (MQ=255)
                       cGAGGCGTCGCATTCTTCAGATGGTTCNACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActt                                                   >  1:1424376/1‑76 (MQ=255)
                       cGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActt                                                   >  1:313575/1‑76 (MQ=255)
                            cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCg                                              <  1:556544/76‑1 (MQ=255)
                              tCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAt                                            <  1:1000165/76‑1 (MQ=255)
                                   ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                       <  1:555919/76‑1 (MQ=255)
                                   ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                       <  1:955802/76‑1 (MQ=255)
                                    tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                      >  2:654647/1‑76 (MQ=255)
                                       tcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAGTCGTAAGCg                                   <  2:489729/76‑1 (MQ=255)
                                         aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgag                                 <  1:722353/76‑1 (MQ=255)
                                             ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                             <  2:458357/76‑1 (MQ=255)
                                                 cAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAGCCGCCt                         <  1:87885/76‑1 (MQ=255)
                                                   aCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTNTTCGATTCGTAAGCGAGAACTGCCTTa                       <  1:1417166/76‑1 (MQ=255)
                                                      cTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAAc                    <  1:531487/76‑1 (MQ=255)
                                                        tAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTc                  >  2:1085530/1‑76 (MQ=255)
                                                           gTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCAt               >  1:424970/1‑76 (MQ=255)
                                                            ttAGCGCTTATGGGTTTTACCCCCGCGTGATGGCTACTTCTTCGATTCGTAAGCGAGNACTGCCTTAAACTCCATa              <  1:1355202/76‑1 (MQ=25)
                                                             tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc             >  1:483745/1‑76 (MQ=255)
                                                             tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc             <  1:387369/76‑1 (MQ=255)
                                                                   ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt       >  2:252482/1‑76 (MQ=255)
                                                                   ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt       >  1:1014694/1‑76 (MQ=255)
                                                                    tATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtc      >  2:1014511/1‑76 (MQ=255)
                                                                    tATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtc      >  2:125322/1‑76 (MQ=255)
                                                                     aTGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtct     <  1:377792/76‑1 (MQ=255)
                                                                       gggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAATCGAGAACTGCCTTAAACTCCATACCACGTTctct   <  2:15617/76‑1 (MQ=255)
                                                                       gggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTctct   >  2:813134/1‑76 (MQ=255)
                                                                        ggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTg  <  2:148745/76‑1 (MQ=255)
                                                                        ggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCANACCACGTtctacg  >  1:177116/1‑73 (MQ=255)
                                                                         |                                                                          
GAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTG  >  REL606/2774123‑2774270

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