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

GATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGC  >  REL606/15351‑15421
                                   |                                   
gttccaggcgcgttattttgcaaccggggccgcaac                                     <  1:4228150‑M1/1‑1 (MQ=255)
gttccaggcgcgttattttgcaaccggggccgcaac                                     >  1:6150680‑M1/36‑36 (MQ=255)
gATATGCGAGCACCTGAAAAGTGGCGGGGGTCACTc                                     >  1:3559699/1‑36 (MQ=255)
gATATGCGAGCACCTGAAAAGTGGCGGGGATTACTc                                     >  1:7006107/1‑36 (MQ=255)
 ttccaggcgcgttattttgcaaccggggccgcaacc                                    <  1:4662307‑M1/2‑1 (MQ=255)
  tccaggcgcgttattttgcaaccggggccgcaaccc                                   >  1:1498893‑M1/34‑36 (MQ=255)
  tatGCGAGCACCTGAAAAGTGGCGGGGATCNCTccc                                   >  1:6101917/1‑36 (MQ=255)
  tatGCGAGCACCTGAAAAGTGGAGGGGATCACTccc                                   <  1:3213389/36‑1 (MQ=255)
   atGCGAGCACCTGAAAAGTGGCGGGGGTCACGcccc                                  >  1:4933933/1‑35 (MQ=37)
    tGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAt                                 <  1:4950547/36‑1 (MQ=255)
      cGAGCACCTGAAAAGTGGCGGGGGTTACTCCCCTaa                               >  1:1898479/1‑36 (MQ=25)
      cGAGCACCTGAAAAGTGGCGGGGATCACTCCCATaa                               >  1:6978048/1‑36 (MQ=255)
      cGAGCACCTGAAAAGTGGCGGGGATCACTCCCATaa                               <  1:2932187/36‑1 (MQ=255)
       gAGCACCTGACAAGTGGCGGGGATCACTCCCATAAg                              <  1:395846/36‑1 (MQ=255)
       gAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAg                              <  1:5156058/36‑1 (MQ=255)
        aGCGCCTGAAAAGTGGCGGGGATCACTCCCATAAgc                             <  1:222513/36‑1 (MQ=255)
           aCCTGAAAAGTGGCGGGGATGACTCCCATAAGCGCt                          <  1:6532806/36‑1 (MQ=255)
           aCCTGAAAAGTGGCGGCGATAAATACCCCAAGCGCt                          >  1:211195/1‑36 (MQ=9)
             tattttgcaaccggcgccgcaaCCCATAAGCGCTaa                        >  1:4254583‑M1/23‑36 (MQ=255)
             cTGAAAAGTGGCGGGGATCACTCCCCTAAGCGCTaa                        >  1:6664718/1‑36 (MQ=255)
             cTGAAAAGTGGCGGGGATCACTCCCATAAGCGCtta                        >  1:823164/1‑34 (MQ=255)
             cTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTaa                        <  1:3095841/36‑1 (MQ=255)
              tGAAAAGTGGCTGGGATCACTCCCATAAGCGCTAAc                       <  1:2090775/36‑1 (MQ=255)
              tGAAAAGTGGCGGGGGTCACTCCCCTAAGCGCTAAc                       >  1:6337049/1‑36 (MQ=25)
              tGAAAAGTGGCGGGGGTCACTCCCATAAGCGCTAAc                       >  1:1736564/1‑36 (MQ=255)
              tGAAAAGTGGCGGGGGTCACTCCCATAAACGCTAAc                       >  1:6664433/1‑36 (MQ=25)
              tGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAAc                       <  1:5034538/36‑1 (MQ=255)
              attttgcaaccggggccccaaCCCATAAGCGCTAAc                       <  1:144910‑M1/15‑1 (MQ=255)
               gAAAAGTGGCGGGGGTCACTCCCCCAAACGCTCACt                      >  1:3815567/1‑36 (MQ=16)
               gAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                      <  1:2064755/36‑1 (MQ=255)
                tttgcaaccggggccgcaaCCCATAAGCGCTCACtt                     <  1:2686897‑M1/17‑1 (MQ=255)
                aaaaGTGGCGGGGATCACTCCCATAAGCGCTAACtt                     <  1:914535/36‑1 (MQ=255)
                 aaaGTGGCGGGGATCACTCCCATAAGCGCTAACTTa                    <  1:6424649/36‑1 (MQ=255)
                  aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTaa                   >  1:1251117/1‑36 (MQ=255)
                   aGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAg                  <  1:4996592/36‑1 (MQ=255)
                     aaccggggccgcaaCCCATAAGCGCTAACTTAAggg                <  1:696160‑M1/22‑1 (MQ=255)
                      ggCGGGGATCACTCCCATCAGCGCTAACTTAAGGGt               >  1:33167/1‑36 (MQ=255)
                      ggCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt               >  1:5190177/1‑36 (MQ=40)
                       gCGGGGATCACTCCCCTAAGCGCTAACTTAAGGGtt              >  1:1910780/1‑36 (MQ=255)
                       gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt              <  1:5493910/36‑1 (MQ=38)
                       gCGGGGATCACTCCCAGAAGCGCTAAGTTAAggggt              >  1:4559429/1‑34 (MQ=25)
                        cGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg             >  1:1019849/1‑36 (MQ=38)
                         ggggATCACTCCCATAAGCGCTAACTTAAGGGTtgt            <  1:2565065/36‑1 (MQ=40)
                         ggggATCACTCCCATAAGCGCTAACTTAAGGGTtgt            <  1:642617/36‑1 (MQ=40)
                          gggATCACTCCCATAAGCGCTAACTTAAGGGTtgtg           >  1:1777229/1‑36 (MQ=255)
                           ggATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg          >  1:495953/1‑36 (MQ=255)
                            gATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGt         <  1:6535136/36‑1 (MQ=255)
                             aTCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTa        >  1:5967242/1‑36 (MQ=255)
                              tCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt       <  1:6215573/36‑1 (MQ=40)
                              tCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAt       <  1:5487054/36‑1 (MQ=40)
                               gcaaCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt      <  1:1569491‑M1/32‑1 (MQ=255)
                               cACTCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt      >  1:1534059/1‑36 (MQ=38)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     >  1:3898607/1‑36 (MQ=35)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     <  1:6326821/36‑1 (MQ=35)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     >  1:488969/1‑36 (MQ=35)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     >  1:5642123/1‑36 (MQ=35)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGGGGTATTa     >  1:2561690/1‑36 (MQ=255)
                                 cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc    >  1:4978131/1‑36 (MQ=32)
                                 cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc    <  1:1636909/36‑1 (MQ=32)
                                 aaCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc    <  1:6830039‑M1/34‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:4015424‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:6591994‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:2394326‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:5704078‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:5028224‑M1/1‑36 (MQ=255)
                                   |                                   
GATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGC  >  REL606/15351‑15421

Base quality scores:  ATCG  < 3 ≤  ATCG  < 8 ≤  ATCG  < 19 ≤  ATCG  < 30 ≤  ATCG  < 37 ≤  ATCG