Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp 22.4% 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 15386 =NA (NA)6 (0.110) 6/64 NT 17.6% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 28 (0.490)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 16728NA (NA)9 (0.170) 8/64 NT 27.3% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 24 (0.420)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562
                                   |                                   
tggtttaacccttaagttagcgcttatgggCCGCaa                                     <  1:7201532‑M2/6‑1 (MQ=255)
cGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCaa                                     <  1:1021727/36‑1 (MQ=255)
 tgtCCAGTCGCGTTATTTGGCACCCGGGGCCGCAAc                                    <  1:1368918/34‑1 (MQ=25)
 gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                    >  1:6150680‑M2/1‑35 (MQ=255)
 gTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAc                                    <  1:4228150‑M2/36‑2 (MQ=255)
  ttCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                   <  1:4662307‑M2/36‑3 (MQ=255)
  gttcaacccttaagttagcgcttatgggCCGCAAcc                                   <  1:7332263‑M2/8‑1 (MQ=255)
   tCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                  >  1:1498893‑M2/1‑33 (MQ=255)
    tcaacccttaagttagcgcttatgggCCGCAAcccc                                 >  1:4148069‑M2/27‑36 (MQ=255)
      aGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCgg                               <  1:6355333/36‑1 (MQ=255)
       ggCGCGTTATTTTGCAACCGGGGCCGCAACCCCggg                              <  1:1898923/36‑1 (MQ=255)
       ggCGCGTTATTTTGCAACCGGGGCCGCAACCCCggg                              <  1:3240003/36‑1 (MQ=255)
        gcgcGTTATTTTGCAACCGGGGCCGCAACCCCGGGt                             >  1:860755/1‑36 (MQ=255)
        gcgcGTTATTTTGCAACCGGGGCCGCAACCCCGGGt                             <  1:563167/36‑1 (MQ=255)
        gcgcGTTATTTTGCAACCGGGGCCGCAACCCCGGGt                             <  1:442249/36‑1 (MQ=255)
        cccttaagttagcgcttatgggCCGCAACCCCGGGt                             <  1:5014619‑M2/14‑1 (MQ=255)
          gcgTTATTTTGCAACCGGGGCCGCAACTCCGGGtgc                           >  1:5426700/1‑36 (MQ=255)
          gcgTTATTTTGCAACCGGGGCCGCAACCCCGGGtgc                           <  1:1348133/36‑1 (MQ=255)
          gcgTTATTTTGCAACCGGGGCCGCAACCCCGGGtgc                           >  1:5517369/1‑36 (MQ=255)
           cgTTATTTTGCAACCGGGGCCGCAACCCCGGGtgct                          <  1:6836219/36‑1 (MQ=255)
            gTTATTTTGCACCCGGGGCCGCAACCCCGGGtgctg                         >  1:4701878/1‑36 (MQ=255)
              tATTTTGCAACCGGCGCCGCAAcccataagcgctaa                       >  1:4254583‑M2/1‑22 (MQ=255)
               aTTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTa                      >  1:3876791/1‑36 (MQ=255)
               aTTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTa                      <  1:1877201/36‑1 (MQ=255)
               aTTTTGCAACCGGGGCCCCAAcccataagcgctaac                      <  1:144910‑M2/36‑16 (MQ=255)
                 tttGCAACCGGGGCCGCAAcccataagcgctcactt                    <  1:2686897‑M2/36‑18 (MQ=255)
                 tttGCAACCGGGGCCGCAACCCCGGGTGCTGCTAAt                    <  1:4766822/36‑1 (MQ=255)
                 tttGCAACCGGGGCCGCAACCCCGGGTGCTGCTAAt                    >  1:4669705/1‑36 (MQ=255)
                     cAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCgg                <  1:4727434/36‑1 (MQ=255)
                      cttatgggCCGCAACCCCGGGTGCTGCTAATGCGGa               <  1:1372558‑M2/28‑1 (MQ=255)
                      aaCCGGGGCCGCAAcccataagcgctaacttaaggg               <  1:696160‑M2/36‑23 (MQ=255)
                      aaCCGGGGCCGCAATCCCGGGAGCTGCTAATGCGGa               >  1:1124410/1‑36 (MQ=37)
                      aaCCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGa               <  1:1521566/36‑1 (MQ=255)
                       aCCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGAt              >  1:6245911/1‑36 (MQ=255)
                        cctgggCCGCAACCCCGGGTGCTGCTAATGCGGATg             <  1:6440325‑M2/30‑1 (MQ=255)
                         cGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGc            >  1:4191695/1‑36 (MQ=255)
                          ggggCCGCAACCCCGGGTGCTGCTAATGCGGATGCg           <  1:2544800/36‑1 (MQ=255)
                             gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc        >  1:1499155‑M2/2‑36 (MQ=255)
                             gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAcc        >  1:3995582‑M2/2‑36 (MQ=255)
                              ccGCAACCCCGGGTGCTGCTAATGCGGATGCGACCt       <  1:7087568/36‑1 (MQ=255)
                               cGCAACCCCGGGTGCTGCTAATGCGGATGCGACCtt      <  1:2926567/36‑1 (MQ=255)
                                gCAAcccataagcgctaacttaagggttgtggtatt     <  1:1569491‑M2/36‑33 (MQ=255)
                                gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc     >  1:1729861/1‑36 (MQ=255)
                                gCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc     <  1:3791635/36‑1 (MQ=255)
                                 cAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa    <  1:4258775/36‑1 (MQ=255)
                                 cAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCa    <  1:6101522/36‑1 (MQ=255)
                                  aacccataagcgctaacttaagggttgtggtattac   <  1:6830039‑M2/36‑35 (MQ=255)
                                   aCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAg  <  1:3621869/36‑1 (MQ=255)
                                   |                                   
CGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAG  >  REL606/4524492‑4524562

Base quality scores:  ATCG  < 3 ≤  ATCG  < 15 ≤  ATCG  < 23 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG