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.070) 35/64 0.1 94.7% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 3 (0.220)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 2774196NA (NA)45 (0.890) 31/64 0.2 95.7% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 2 (0.160)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

GGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAT  >  REL606/2772820‑2772884
                                  |                              
ggCGCTGGGAACACCCTGAACGGTGGGGAAAAAccc                               <  2:2281835/36‑1 (MQ=255)
ggCGCGGGGAACACCCTGAACGGTGGGGAAAAAccc                               <  1:1417601/36‑1 (MQ=255)
 gcgcGGGGAACACCCTGAACGGTGGGGAAAAAcccc                              <  2:1258128/36‑1 (MQ=255)
 gcgcGGGGAACACCCTGAACGGTGGGGAAAAAcccc                              <  1:652743/36‑1 (MQ=255)
   gcgGGGGACGCCCTGAACGGTGGGGAAAAACCCCAt                            <  1:1412791/36‑1 (MQ=37)
   gcgGGGAACACCCTGAACGGTGGGGAAAAACCCCAt                            <  1:1226148/36‑1 (MQ=255)
    cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATa                           <  2:1841034/36‑1 (MQ=255)
    cgGGGAACACCCTGAACGGTGGGGAAAAACCCCATa                           <  2:2257541/36‑1 (MQ=255)
     ggggAACACCCTGAACGGTGGGGAAAAACCCCATaa                          >  2:1464967/1‑36 (MQ=255)
      gggAACACCCTGAACGGTGGGGAAAAACCCCATAAg                         <  1:2334768/36‑1 (MQ=255)
       ggAACACCCTGAACGGTGGGGAAAAACCCCATAAgc                        >  2:2820683/1‑36 (MQ=255)
       ggAACACCCTGAACGGTGGGGAAAAACCCCATAAgc                        >  1:469934/1‑36 (MQ=255)
       ggAACACCCTGAACGGCGGGGAAAAACCCCACaacc                        >  1:1877803/1‑34 (MQ=25)
         aaCACCCTGAACGGTGGGGAAAAACCCCATAAgcgc                      <  1:2059401/36‑1 (MQ=255)
          acatCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                     <  2:1103436/36‑1 (MQ=255)
          acacCCTGAACGGTGGGGAAAAACCCCATAAGCGCt                     >  2:2625517/1‑36 (MQ=255)
           cacCCTGAACGGTGGGGAAAAACCCCATAAGCGCTa                    <  2:3585447/36‑1 (MQ=255)
            acCCTGAACGGTGGGGAAAAACCCCATAAGCGCTaa                   >  2:3533059/1‑36 (MQ=255)
             cccTGAACGGTGGGGAAACACCCCATAAGCGCTAAc                  >  2:3350500/1‑36 (MQ=255)
               cagAACAGTGAGGAAAAACCCCATAAGCGCTAACtt                <  1:722775/34‑1 (MQ=18)
               cTGAACGGTGGGTAAAAACCCCATAAGCGCTAACtt                <  2:789720/36‑1 (MQ=255)
               cTGAACGGTGGGGAAAAACCCCATAAGCGCTAACtt                <  2:2491391/36‑1 (MQ=255)
                tGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTa               >  1:3048963/1‑36 (MQ=255)
                   aCGGTGGGGAAAAACCCCATAAGCGCTAACTTAAgg            >  2:1470753/1‑36 (MQ=255)
                    cGGTGGGGAAAAACCCCATAAGCGCTAACTTAAggg           >  2:868949/1‑36 (MQ=255)
                     ggTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGt          >  1:2967724/1‑36 (MQ=255)
                      gTGGAGAAAAACCCCATAAGCGCTAACTTAAGGGtt         <  1:637595/36‑1 (MQ=255)
                       tGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTtg        <  1:2124588/36‑1 (MQ=255)
                         gggAGAGACGCCATAAGCGCTAACTTAAGGGTtgtg      <  1:1769417/36‑1 (MQ=11)
                         gggAAAAACCCCATAAGCGCTAACTTAAGGGTtgtg      <  2:1895844/36‑1 (MQ=255)
                          ggAAAAACCCCATAAGCGCNAACTTAAGGGTTGTgg     >  1:1020958/1‑36 (MQ=255)
                           gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt    <  2:254619/36‑1 (MQ=38)
                           gAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGt    <  1:1058710/36‑1 (MQ=38)
                            aaaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTa   <  2:1512418/36‑1 (MQ=35)
                             aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTAt  >  1:3511048/1‑36 (MQ=32)
                             aaaaCCCCATAAGCGCTAACTTAAGGGTTGTGGTAt  >  1:1296371/1‑36 (MQ=32)
                                  |                              
GGCGCGGGGAACACCCTGAACGGTGGGGAAAAACCCCATAAGCGCTAACTTAAGGGTTGTGGTAT  >  REL606/2772820‑2772884

Base quality scores:  ATCG  < 3 ≤  ATCG  < 13 ≤  ATCG  < 22 ≤  ATCG  < 30 ≤  ATCG  < 32 ≤  ATCG