Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp 84.3% 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)40 (0.710) 29/64 NT 78.4% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 11 (0.180)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 16728NA (NA)46 (0.810) 31/64 NT 90.2% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 5 (0.080)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

Only 100 of 112 total aligned reads displayed.

GATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGC  >  REL606/15351‑15421
                                   |                                   
gttccaggcgcgttattttgcaaccggggccgcaac                                     <  1:3566822‑M1/1‑1 (MQ=255)
gttccaggcgcgttattttgcaaccggggccgcaac                                     <  1:3306910‑M1/1‑1 (MQ=255)
  tccaggcgcgttattttgcaaccggggccgcaaccc                                   <  1:3676004‑M1/3‑1 (MQ=255)
  tccaggcgcgttattttgcaaccggggccgcaaccc                                   >  1:4394492‑M1/34‑36 (MQ=255)
  tccaggcgcgttattttgcaaccggggccgcaaccc                                   >  1:4537539‑M1/34‑36 (MQ=255)
  tatGCGAGCACCTGAAAAGTGGCGGGGATCACTccc                                   <  1:6765983/36‑1 (MQ=255)
   ccaggcgcgttattttgcaaccggggccgcaaCCCa                                  <  1:3420464‑M1/4‑1 (MQ=255)
   atGCGAGCACCTGAAAAGTGGCGGGGGTCACTtccc                                  >  1:3654043/1‑35 (MQ=37)
    tGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAt                                 <  1:5016016/36‑1 (MQ=255)
    caggcgcgttattttgcaaccggggccgcaaCCCAt                                 >  1:3685538‑M1/32‑36 (MQ=255)
    caggcgcgttattttgcaaccggggccgcaaCCCAt                                 >  1:8252873‑M1/32‑36 (MQ=255)
     gCGAGCACCTGAAAAGTGGCGGGGGTCACGCCCCTa                                >  1:201404/1‑36 (MQ=25)
     aggcgcgttattttgcaaccggggccgcaaCCCATa                                >  1:460962‑M1/31‑36 (MQ=255)
     aggcgcgttattttgcaaccggggccgcaaCCCATa                                >  1:3571826‑M1/31‑36 (MQ=255)
      ggcgcgttattttgcaaccggggccgcaaCCCATaa                               >  1:5405410‑M1/30‑36 (MQ=255)
      ggcgcgttattttgcaaccggggccgcaaCCCATaa                               >  1:4524361‑M1/30‑36 (MQ=255)
      cGAGCACCTGAAAAGTGGCGGGGGTCACTCCCCTaa                               >  1:6625362/1‑36 (MQ=37)
      cGAGCACCTGAAAAGTGGCGGGGATCACTCCCATaa                               <  1:1211579/36‑1 (MQ=255)
       gAGCACCTGAAAAGTGGCGGGGGTCACTCCCCTAAg                              >  1:4448177/1‑36 (MQ=37)
       gAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAg                              <  1:7876528/36‑1 (MQ=255)
       gAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAg                              <  1:3767140/36‑1 (MQ=255)
        cgcgttattttgcaaccggggccgcaaCCCATAAgc                             >  1:6764757‑M1/28‑36 (MQ=255)
        aGCACCTGAAAAGTGGCGGGGGTCACCCCCCTaaac                             >  1:2686632/1‑34 (MQ=21)
        aGCACCTGAAAAGTGGCGGGGATCACTCCCATAAgc                             <  1:2646543/36‑1 (MQ=255)
          cACCTGAAAAGTGGCGGGGATCACTCCCATAAgcgc                           <  1:3505873/36‑1 (MQ=255)
           gttattttgcaaccggggccgcaaCCCATaagcgat                          <  1:1616650‑M1/12‑3 (MQ=255)
           aCCTGAAAAGTGGCGGGGATCACTCCCCTAAGCGCt                          >  1:6642108/1‑36 (MQ=255)
           aCCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCt                          <  1:3892770/36‑1 (MQ=255)
           aCCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCt                          <  1:3379686/36‑1 (MQ=255)
            ttattttgcaaccggggccgcaaCCCATAAGCGCTa                         >  1:4806779‑M1/24‑36 (MQ=255)
             tattttgcaaccggggccgcaaCCCATAAGCGCTaa                        >  1:4282991‑M1/23‑36 (MQ=255)
             cTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTaa                        <  1:4796802/36‑1 (MQ=255)
              tGAAAAGTGGCGGGGGTCACTCCCCTAAACGCTAAc                       >  1:7288992/1‑36 (MQ=37)
              tGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAAc                       <  1:7669167/36‑1 (MQ=255)
               ttttgcaaccggggccgcaaCCCATAAGCGCTAACt                      <  1:4002364‑M1/16‑1 (MQ=255)
               gAAAAGTGGCGGGGGTCACTCCCCTAAGCGCTAACt                      >  1:6405546/1‑36 (MQ=21)
               gAAAAGTGGCGGGGGTCACTCCCATAAGCGCTAACt                      >  1:278124/1‑36 (MQ=255)
                ttttcaaccggggccgcaaCCCATAAGCGCTAACtt                     >  1:3711554‑M1/20‑36 (MQ=255)
                tttgcaaccggggccgcaaCCCATAAGCGCTACCtt                     >  1:2357130‑M1/20‑36 (MQ=255)
                tttgcaaccggggccgcaaCCCATAAGCGCTAACtt                     <  1:1563108‑M1/17‑1 (MQ=255)
                aaaaGTGGCGGGGTTGACTCCCTTAAGAGCGAAAtt                     >  1:6559300/1‑36 (MQ=12)
                aaaaGTGGCGGGGATCACTCCCATAAGCGCTAACtt                     <  1:4062296/36‑1 (MQ=255)
                 ttgcaaccggggccgcaaCCCATAAGCGCTAACTTa                    <  1:4365296‑M1/18‑1 (MQ=255)
                 ttgcaaccggggccgcaaCCCATAAGCGCTAACTTa                    >  1:481770‑M1/19‑36 (MQ=255)
                 aaaGTGGCGGGGGTCACTCCCCCAAGCGCTTACTTa                    >  1:2611131/1‑36 (MQ=255)
                  tgcaaccggggccgcaaCCCATAAGCGCTAACTTaa                   >  1:7783020‑M1/18‑36 (MQ=255)
                  aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTaa                   <  1:5316996/36‑1 (MQ=255)
                   gcaaccggggccgcaaCCCATAAGCGCTAACTTAAg                  <  1:605835‑M1/20‑1 (MQ=255)
                   aGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAg                  <  1:2411029/36‑1 (MQ=255)
                   aGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAg                  <  1:5327863/36‑1 (MQ=255)
                    caaccggggccgcaaCCCATAAGCGCTAACTTAAgg                 <  1:4119728‑M1/21‑1 (MQ=255)
                     tgtcGGGGATCACTCCCATAAGCGCTAACTTAAggg                <  1:477768/33‑1 (MQ=255)
                     tGGCGGGGATCACGCCCATAAGCGCTAACTTAAggg                >  1:2409040/1‑36 (MQ=255)
                     aaccggggccgcaaCCCATAAGCGCTAACTTAAggg                <  1:498486‑M1/22‑1 (MQ=255)
                     aaccggggccgcaaCCCATAAGCGCTAACTTAAggg                >  1:6794420‑M1/15‑36 (MQ=255)
                      ggCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt               >  1:689926/1‑36 (MQ=40)
                      ggCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt               >  1:8356848/1‑36 (MQ=40)
                       gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt              <  1:7242929/36‑1 (MQ=38)
                       gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt              <  1:5356032/36‑1 (MQ=38)
                       acGGGGATCACTCCCATCAGCTCTACCTTACGGGtt              <  1:6134512/35‑1 (MQ=12)
                        cggggccgcaaCCCATAAGCGCTAACTTAAGGGTtg             <  1:72552‑M1/25‑1 (MQ=255)
                        cGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg             >  1:6185948/1‑36 (MQ=38)
                        cGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg             >  1:1556715/1‑36 (MQ=38)
                         ggggATCACTCCCATAAGCGCTAACTTAAGGGTtgt            <  1:6041742/36‑1 (MQ=40)
                         ggggATCACTCCCATAAGCGCTAACTTAAGGGTtgt            >  1:2381019/1‑36 (MQ=40)
                          gggccgcaaCCCATAAGCGCTAACTTAAGGGTtgtg           >  1:3513277‑M1/10‑36 (MQ=255)
                          gggATCACTCCCATAAGCGCTAACTTAAGGGTtgtg           >  1:4570486/1‑36 (MQ=255)
                           ggATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg          <  1:5924076/36‑1 (MQ=255)
                           ggATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg          >  1:3789823/1‑36 (MQ=255)
                            gccgcaaCCCATAAGCGCTAACTTAAGGGTTGTGGt         >  1:5358454‑M1/8‑36 (MQ=255)
                            gccgcaaCCCATAAGCGCTAACTTAAGGGTTGTGGt         >  1:2011765‑M1/8‑36 (MQ=255)
                            gccgcaaCCCATAAGCGCTAACTTAAGGGTTGTGGt         >  1:1409268‑M1/8‑36 (MQ=255)
                            gATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGt         <  1:2358895/36‑1 (MQ=255)
                             ccgcaaCCCATAAGCGCTAACTTAAGGGTTGTGGTa        <  1:4647177‑M1/30‑1 (MQ=255)
                             ccgcaaCCCATAAGCGCTAACTTAAGGGTTGTGGTa        <  1:4212182‑M1/30‑1 (MQ=255)
                              tCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt       <  1:1660806/36‑1 (MQ=40)
                              cgcaaCCCATAAGCGCTAACTTAAGGGTTGTGGTAt       >  1:6013391‑M1/6‑36 (MQ=255)
                               gcaaCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt      >  1:4523484‑M1/5‑36 (MQ=255)
                               cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt      >  1:6484789/1‑36 (MQ=38)
                               cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt      <  1:4832122/36‑1 (MQ=38)
                                caaCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     <  1:2407216‑M1/33‑1 (MQ=255)
                                aCTCCCATAAGCGCTCACTTAAGGGTTGTGGTATTa     <  1:63531/36‑1 (MQ=255)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     >  1:8433688/1‑36 (MQ=35)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     >  1:2724458/1‑36 (MQ=35)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     <  1:6104892/36‑1 (MQ=35)
                                 aaCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc    <  1:7565401‑M1/34‑1 (MQ=255)
                                 aaCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc    <  1:5612767‑M1/34‑1 (MQ=255)
                                  aCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg   <  1:5356340‑M1/35‑1 (MQ=255)
                                   CCCCTAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:6583476‑M1/1‑36 (MQ=255)
                                   CCCATAAGCGCTAACTTACGGGTTGTGGTATTACGc  <  1:6359272‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:61716‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:5583079‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:633820‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:448494‑M1/1‑36 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:7711983‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:4199073‑M1/1‑36 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:8136677‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:1393397‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:1191941‑M1/1‑36 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:1119738‑M1/36‑1 (MQ=255)
                                   |                                   
GATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGC  >  REL606/15351‑15421

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 24 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG