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)41 (1.130) 31/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)39 (1.070) 32/144 0.2 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:1052793/1‑76 (MQ=32)
gAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtt                                                                          >  1:1165748/1‑76 (MQ=32)
 aGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGttt                                                                         >  1:751130/1‑76 (MQ=34)
 aGGCATCGTCTTAACGAGGCACCGACGCGTCGCATTCTTCNNATGGTTCAACNNTAAAGTTAGCGCTTATGGGttt                                                                         >  1:1084039/1‑76 (MQ=11)
   gCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTa                                                                       <  2:433991/76‑1 (MQ=38)
    cATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAc                                                                      <  2:496033/76‑1 (MQ=39)
     aTCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcc                                                                     >  1:955302/1‑76 (MQ=40)
       cGTCTTAACGAGGCACCGANGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcccc                                                                   >  1:1128022/1‑76 (MQ=255)
       cGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAcccc                                                                   <  2:676463/76‑1 (MQ=255)
        gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATTGGTTTTAccccc                                                                  <  2:730577/76‑1 (MQ=255)
        gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAccccc                                                                  >  2:781514/1‑76 (MQ=255)
        gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAccccc                                                                  <  2:716095/76‑1 (MQ=255)
        gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAccccc                                                                  <  2:951382/76‑1 (MQ=255)
        gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAccccc                                                                  <  2:328539/76‑1 (MQ=255)
        gTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAccccc                                                                  >  1:711880/1‑76 (MQ=255)
        gTCTTAACGAGGCACCGAGGCGTCGCATTCATCAGATGGTTCAACCCTTAAGTTAGCGCTTANGGGTTTTAccccc                                                                  >  2:1099480/1‑76 (MQ=14)
             aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtg                                                             >  2:314163/1‑76 (MQ=255)
             aaCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtg                                                             <  1:322363/76‑1 (MQ=255)
               cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggt                                                           >  2:291245/1‑76 (MQ=255)
               cGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggt                                                           <  2:458766/76‑1 (MQ=255)
                gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                          <  2:986086/76‑1 (MQ=255)
                gAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggtg                                                          <  2:894230/76‑1 (MQ=255)
                 aGGCACCGCAGCGTCGCATTCATCANATGGTTCAACCCTTNAGTTAGCTCTNANGGGTTTTACTCCCGCGtggtgg                                                         >  2:20226/1‑76 (MQ=17)
                   gCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTTGCt                                                       <  2:844797/76‑1 (MQ=255)
                    cACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTa                                                      <  2:694627/76‑1 (MQ=255)
                     aCCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTAc                                                     <  2:681469/76‑1 (MQ=255)
                        gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtc                                                  <  2:136187/76‑1 (MQ=255)
                        gAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtc                                                  >  2:586146/1‑76 (MQ=255)
                         aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtct                                                 <  2:790825/76‑1 (MQ=255)
                            cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCg                                              >  1:637578/1‑76 (MQ=255)
                               cGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGAtt                                           >  1:638085/1‑76 (MQ=255)
                                  gttcttcAGATGGNTCAACCCTTACGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGt                                        <  1:157468/75‑1 (MQ=255)
                                  aTTCTTCAGATGGTTCAACCCTTAAATTAGCTCTTATGGCNNTTACCCCCGCNNGGTGGCTACTTCTTCGATTCGt                                        >  1:1093126/1‑76 (MQ=25)
                                   ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                       >  1:1045126/1‑76 (MQ=255)
                                    tcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                      >  1:602578/1‑76 (MQ=255)
                                      ttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGc                                    <  1:512413/76‑1 (MQ=255)
                                         aGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCgag                                 >  2:1067512/1‑76 (MQ=255)
                                             ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                             <  2:889011/76‑1 (MQ=255)
                                                 cAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCt                         >  2:292603/1‑76 (MQ=255)
                                                  aaCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCtt                        <  1:404206/76‑1 (MQ=255)
                                                   aCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTa                       <  2:900977/76‑1 (MQ=255)
                                                       ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACt                   >  2:447796/1‑76 (MQ=255)
                                                         aaGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAACGAGAACTGCCTTAAACTcc                 >  1:1166595/1‑76 (MQ=255)
                                                          aGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCa                <  2:232259/76‑1 (MQ=255)
                                                            ttAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATa              <  1:1179362/76‑1 (MQ=255)
                                                             tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGGTTCGTAAGCGAGGACGGGCGTAAACTCCAtcc             >  2:931985/1‑74 (MQ=21)
                                                             tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc             >  2:931991/1‑76 (MQ=255)
                                                               gcgcTTATNGGTTTGACCCNCGCGTGGTGGCTACTTCTTCGANTCGTAAGCGAGAACTGCCTTAAACTCCATACCa           <  2:1010676/76‑1 (MQ=255)
                                                                cgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCAc          <  1:190779/76‑1 (MQ=255)
                                                                  cTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGt        >  1:371851/1‑76 (MQ=255)
                                                                    tATGGGTTTTACCCCCGCTTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtc      <  1:321564/76‑1 (MQ=255)
                                                                       gggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTctct   <  1:995382/76‑1 (MQ=255)
                                                                        ggTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTg  >  2:451069/1‑76 (MQ=255)
                                                                         |                                                                          
GAGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCTCTG  >  REL606/2774123‑2774270

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