Predicted mutation
evidence seq id position mutation freq annotation gene description
JC JC REL606 4,524,522 IS186 (+) +6 bp 61.1% 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)16 (0.290) 14/64 NT 50.0% noncoding (1/1343 nt) IS186 repeat region
?REL606 = 4524527 16 (0.260)coding (499/549 nt) fimA major type 1 subunit fimbrin (pilin)
* ? REL606 = 16728NA (NA)28 (0.500) 20/64 NT 70.0% noncoding (1343/1343 nt) IS186 repeat region
?REL606 4524522 = 12 (0.200)coding (494/549 nt) fimA major type 1 subunit fimbrin (pilin)

GATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGC  >  REL606/15351‑15421
                                   |                                   
gttccaggcgcgttattttgcaaccggggccgcaac                                     <  1:1535038‑M1/1‑1 (MQ=255)
gttccaggcgcgttattttgcaaccggggccgcaac                                     <  1:3848795‑M1/1‑1 (MQ=255)
gttccaggcgcgttattttgcaaccggggccgcaac                                     <  1:3284932‑M1/1‑1 (MQ=255)
gATATGCGAGCACCTGAAAAGTGGCGGGGGTTACTc                                     >  1:1215120/1‑36 (MQ=37)
 atatGCGAGCACCTGAAAAGTGGCGGGGATCACTcc                                    <  1:1457552/36‑1 (MQ=255)
   atGCGAGCACCTGAAAAGCGGCGGGGGTCACTcccc                                  >  1:4243091/1‑35 (MQ=37)
    tGCGAGCACCTGAAAAGTGGCGGGGATCACTccccc                                 >  1:417748/1‑34 (MQ=255)
    tGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAt                                 <  1:6310381/36‑1 (MQ=255)
    tGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCAt                                 <  1:3752435/36‑1 (MQ=255)
    caggcgcgttattttgcaaccggggccgcaaCCCAt                                 <  1:2359398‑M1/5‑1 (MQ=255)
     gCGAGCACCTGAAAAGTGGCGGGGGTTACTCCCCTa                                >  1:4530041/1‑36 (MQ=25)
       gcgcgttattttgcaaccggggccgcaaCCCATAAg                              <  1:4529168‑M1/8‑1 (MQ=255)
       gAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAg                              <  1:2604215/36‑1 (MQ=255)
          cgttattttgcaaccggggccgcaaCCCATAAgcgc                           >  1:1004166‑M1/26‑36 (MQ=255)
          cACCTGAAAAGTGGCGGGGATCACTCCCCTAAgcgc                           >  1:874234/1‑36 (MQ=255)
           aCCTGAAAAGTGGCGGGGGTTACTCCCATTAGCGCt                          >  1:6331922/1‑36 (MQ=18)
            acTGAAAAGTGGCGGGGATGACTCCCCTAAgcgcga                         >  1:7610834/2‑34 (MQ=25)
             cTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTaa                        <  1:2592081/36‑1 (MQ=255)
             cTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTaa                        <  1:1325738/36‑1 (MQ=255)
              tGAAAAGTGGCGGGGGTCACTCCCCTAAGCGCTAAc                       >  1:6260042/1‑36 (MQ=25)
              tGAAAAGTGGCGGGGATCACTCCCCTAAGCGCTAAc                       >  1:6751898/1‑36 (MQ=255)
              attttgcaaccggggccgcaaCCCATAAGCGCTAAc                       >  1:5267936‑M1/22‑36 (MQ=255)
              attttgcaaccggggccgcaaCCCATAAGCGCTAAc                       <  1:1277342‑M1/15‑1 (MQ=255)
               ttttgcaaccggggccgcaaCCCATAAGCGCTAACt                      >  1:3026884‑M1/21‑36 (MQ=255)
               gAAAAGTGGCGGGGATCACTCCCCTAAGCGCTAACt                      >  1:6399410/1‑36 (MQ=255)
               gAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                      <  1:7806129/36‑1 (MQ=255)
               gAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACt                      <  1:3341309/36‑1 (MQ=255)
                tttgcaaccggggccgcaaCCCATAAGCGCTAACtt                     >  1:7105572‑M1/20‑36 (MQ=255)
                tttgcaaccggggccgcaaCCCATAAGCGCTAACtt                     >  1:233194‑M1/20‑36 (MQ=255)
                aaaaGTGGCGGGGATCACTCCCATAAGCGCTAACtt                     <  1:1755892/36‑1 (MQ=255)
                  aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTaa                   >  1:751385/1‑36 (MQ=255)
                  aaGTGGCGGGGATCACTCCCATAAGCGCTAACTTaa                   <  1:5170402/36‑1 (MQ=255)
                   aGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAg                  <  1:4647226/36‑1 (MQ=255)
                    gTGGCGGGGATCACTCCCATAAGCGCTAACTTAAgg                 <  1:3884710/36‑1 (MQ=255)
                    caaccggggccgcaaCCCATAAGCGCTAACTTAAgg                 >  1:650778‑M1/16‑36 (MQ=255)
                     tGGCGGGGATCACTCCCATAAGCGCTAACTTAAggg                >  1:5133702/1‑36 (MQ=255)
                     tGGCGGGGATCACTCCCATAAGCGCTAACTTAAggg                >  1:6914069/1‑36 (MQ=255)
                     aaccggggccgcaaCCCATAAGCGCTAACTTAAggg                <  1:5732424‑M1/22‑1 (MQ=255)
                      ggCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt               >  1:6005593/1‑36 (MQ=40)
                      ggCGGGGATCACTCCCATAAGCGCTAACTTAAGGGt               >  1:3292292/1‑36 (MQ=40)
                       gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt              <  1:5004862/36‑1 (MQ=38)
                       gCGGGGATCACTCCCATAAGCGCTAACTTAAGGGtt              <  1:7795936/36‑1 (MQ=38)
                        cGGGGATCACTCCCATAAGCGCTAACTTAAGGGTtg             <  1:419905/36‑1 (MQ=38)
                         ggggccgcaaCCCATAAGCGCTAACTTAAGGGTtgt            <  1:5018394‑M1/26‑1 (MQ=255)
                         ggggATCACTCCCATAAGCGCTAACTTAAGGGTtgt            <  1:7784455/36‑1 (MQ=40)
                         ggggATCACTCCCATAAGCGCTAACTTAAGGGTtgt            >  1:3039171/1‑36 (MQ=40)
                         ggggATCACTCCCATAAGCGCTAACTTAAGGGTtgt            >  1:6177146/1‑36 (MQ=40)
                           ggATCACTCCCATAAGCGCTAACTTAAGGGTTGTgg          <  1:788502/36‑1 (MQ=255)
                            gATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGt         <  1:5080852/36‑1 (MQ=255)
                             aTCACTCCCATAAGCGCTACCTTAAGGGTTGTGGTa        <  1:1787155/36‑1 (MQ=255)
                             aTCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTa        >  1:5391697/1‑36 (MQ=255)
                             aTCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTa        >  1:732417/1‑36 (MQ=255)
                              tCACTCCCCTAAGCGCTAACTTAAGGGTTGTGGTAt       <  1:3261679/36‑1 (MQ=255)
                               gcaaCCCATAAGCGCTAACTTAAGGGTTGTGGTAtt      <  1:5168277‑M1/32‑1 (MQ=255)
                                caaCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     >  1:2905070‑M1/4‑36 (MQ=255)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     >  1:4317257/1‑36 (MQ=35)
                                aCTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTa     <  1:417214/36‑1 (MQ=35)
                                 cTCCCATAAGCGCTAACTTAAGGGTTGTGGTAttgc    <  1:1527932/36‑3 (MQ=17)
                                 cTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc    <  1:2168972/36‑1 (MQ=32)
                                 aaCCCATAAGCGCTAACTTAAGGGTTGTGGTATTAc    <  1:7551936‑M1/34‑1 (MQ=255)
                                  aCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACg   <  1:1040859‑M1/35‑1 (MQ=255)
                                  aCCCATAAGCGCTAACTTAAGGGATGTGGTATTACg   <  1:3574979‑M1/35‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTAcga  <  1:5301247‑M1/36‑2 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:2263013‑M1/1‑36 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:3523401‑M1/1‑36 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:7604008‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:7509256‑M1/1‑36 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:1474564‑M1/1‑36 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:6919597‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:4459368‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  <  1:4463674‑M1/36‑1 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:1073628‑M1/1‑36 (MQ=255)
                                   CCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGc  >  1:5588069‑M1/1‑36 (MQ=255)
                                   |                                   
GATATGCGAGCACCTGAAAAGTGGCGGGGATCACTCCCATAAGCGCTAACTTAAGGGTTGTGGTATTACGC  >  REL606/15351‑15421

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