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)54 (1.110) 40/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)41 (0.850) 35/144 0.4 100% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 0 (0.000)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

AGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCT  >  REL606/2774124‑2774267
                                                                        |                                                                       
aGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGttt                                                                      >  2:139/1‑76 (MQ=34)
aGGCATCGTCTTAAAGAGGCACCGGGGCATCNCATTCTTCTGATGGNTCAACNCTTAAGTTAGCGCTTATGGGttt                                                                      <  2:202485/76‑1 (MQ=11)
 ggCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGtttt                                                                     <  2:717644/76‑1 (MQ=35)
       gTCTTAACGAGTCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTAccccc                                                               <  1:694729/76‑1 (MQ=255)
         cTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCcgc                                                             >  1:541066/1‑76 (MQ=255)
           tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgt                                                           >  2:818754/1‑76 (MQ=255)
           tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgt                                                           <  1:1409766/76‑1 (MQ=255)
           tAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgt                                                           >  2:705884/1‑76 (MQ=255)
             aCGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtgg                                                         <  2:1401482/76‑1 (MQ=255)
              cGAGGCACCGAGGCGTCGCNTTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCgtggt                                                        >  1:1510911/1‑76 (MQ=255)
                 ggCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGc                                                     <  1:1255330/76‑1 (MQ=255)
                        aGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtct                                              <  1:296365/76‑1 (MQ=255)
                         ggCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTActtctt                                             <  1:1305225/76‑1 (MQ=255)
                           cGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCg                                           >  1:1070230/1‑76 (MQ=255)
                            gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCgg                                          >  1:741368/1‑75 (MQ=255)
                            gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGa                                          >  1:741339/1‑76 (MQ=255)
                            gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGa                                          >  1:958018/1‑76 (MQ=255)
                            gTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGa                                          >  2:1066621/1‑76 (MQ=255)
                                cATTCTGCNGCTGGGTCAANCCTTAAGTTAGCGCTTATGGGTNTTACCCCCGCGTGGTGGCTACTTCTTCGATTCg                                      <  2:1347981/76‑1 (MQ=25)
                                  ttcttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTa                                    <  1:973756/76‑1 (MQ=255)
                                   tcttnAGATGGTTCAACCCTTAAGTTAGGGCCTATGGGTNTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                   <  1:193108/76‑1 (MQ=37)
                                   tcttcAGATGGTTNAACCCTTAAGTTAGNGCNGATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTaa                                   <  2:1453108/76‑1 (MQ=255)
                                    cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAg                                  <  1:646390/76‑1 (MQ=255)
                                    cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAg                                  <  1:646375/76‑1 (MQ=255)
                                    cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAg                                  >  2:629994/1‑76 (MQ=255)
                                    cttcAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGGGGCTACTTCTTTGGTTCGTAAg                                  >  1:1269834/1‑76 (MQ=25)
                                       cAGATGGTTCAACCCTTAAGTTAGCGCTTATGGNTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCga                               >  2:1376159/1‑76 (MQ=255)
                                            ggTTCAACCCTTACGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                          >  2:1157159/1‑76 (MQ=255)
                                            ggTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACt                          <  1:1102864/76‑1 (MQ=255)
                                             gTTCAACCCTTAAGTTAGCGCTTATGGGGTTTACCCCCGCGTGGGGGCTACATCTTCGAATCGTGAGCAAGGACtc                         >  2:628838/1‑75 (MQ=18)
                                              ttCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGc                        >  2:967610/1‑76 (MQ=255)
                                                 aaCNCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGNTACTTCTTCGATTCGTAAGCGAGAACTGCCtt                     <  2:1510911/76‑1 (MQ=255)
                                                 aaCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCtt                     <  2:1159266/76‑1 (MQ=255)
                                                 aaCCATTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCtt                     >  2:963101/1‑76 (MQ=255)
                                                  aCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTa                    <  1:1489253/76‑1 (MQ=255)
                                                  aCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTa                    >  2:1080635/1‑76 (MQ=255)
                                                  aCCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTa                    <  1:275580/76‑1 (MQ=255)
                                                    ccTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAATCGAGAACTGCCTTaaa                  <  1:410202/76‑1 (MQ=255)
                                                      ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTTTTCGCTTCGTCAGCGAGAACTGCCTTAAACt                >  2:1180956/1‑76 (MQ=25)
                                                      ttAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACt                >  1:354739/1‑76 (MQ=255)
                                                       tAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTc               >  2:1228170/1‑76 (MQ=255)
                                                         aGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGTCCTCTTCTTCGANTCATAAGCGAGAACTGCCTTAAACTCCa             >  1:139628/1‑76 (MQ=37)
                                                         aGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCa             <  1:769924/76‑1 (MQ=255)
                                                          gTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCAt            <  1:240173/76‑1 (MQ=255)
                                                          attAGCGCTTATGNGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCAt            <  1:204527/75‑1 (MQ=255)
                                                           ttAGCGCGTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATa           <  1:773096/76‑1 (MQ=255)
                                                            tAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATAc          <  2:597330/76‑1 (MQ=255)
                                                             aGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTNAGCGAGAACTGCCTTAAACTCCATAcc         <  1:1605189/76‑1 (MQ=255)
                                                               cgcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCAc       >  1:909337/1‑76 (MQ=255)
                                                                gcTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGCACTGCCTTACACTCCATACCACg      >  1:1394367/1‑76 (MQ=255)
                                                                  ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt    >  1:254570/1‑76 (MQ=255)
                                                                  ttATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGtt    >  2:1410616/1‑76 (MQ=255)
                                                                   tATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtc   >  1:948402/1‑76 (MQ=255)
                                                                    aTGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTtct  >  1:644113/1‑76 (MQ=255)
                                                                        |                                                                       
AGGCATCGTCTTAACGAGGCACCGAGGCGTCGCATTCTTCAGATGGTTCAACCCTTAAGTTAGCGCTTATGGGTTTTACCCCCGCGTGGTGGCTACTTCTTCGATTCGTAAGCGAGAACTGCCTTAAACTCCATACCACGTTCT  >  REL606/2774124‑2774267

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