Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 433,668 C→T 100% R55C (CGT→TGT)  hupB → HU, DNA‑binding transcriptional regulator, beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL606433,6680CT100.0% 283.5 / NA 83R55C (CGT→TGT) hupBHU, DNA‑binding transcriptional regulator, beta subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (47/36);  total (47/36)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAC  >  REL606/433620‑433714
                                                |                                              
ggggATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAgtgt                                              <  1:2240439/51‑1 (MQ=255)
ggggATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAgtgt                                              <  1:7473450/51‑1 (MQ=255)
 gggATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGtgtg                                             >  1:7085997/1‑51 (MQ=255)
  ggATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgc                                            <  1:9743536/51‑1 (MQ=255)
  ggATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgc                                            >  1:10113937/1‑51 (MQ=255)
  ggATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgc                                            <  1:1357150/51‑1 (MQ=255)
    atgatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgct                                           <  1:1174076/50‑1 (MQ=255)
     tgatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgctgc                                         >  1:6363718/1‑51 (MQ=255)
     tgatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgctgc                                         >  1:1947379/1‑51 (MQ=255)
     tgatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGtgctgc                                         >  1:6695412/1‑51 (MQ=255)
      gatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGcc                                        >  1:3354298/1‑51 (MQ=255)
      gatgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGcc                                        >  1:1809454/1‑51 (MQ=255)
       atgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGccc                                       <  1:2652207/51‑1 (MQ=255)
       atgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGccc                                       >  1:3076715/1‑51 (MQ=255)
       atgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGccc                                       <  1:1333190/51‑1 (MQ=255)
       atgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGccc                                       >  1:10590722/1‑51 (MQ=255)
        tgTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCg                                      >  1:1459699/1‑51 (MQ=255)
         gTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGt                                     <  1:10576064/51‑1 (MQ=255)
         gTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGt                                     <  1:3048945/51‑1 (MQ=255)
          tAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTa                                    <  1:6005412/51‑1 (MQ=255)
          tAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTa                                    >  1:10380805/1‑51 (MQ=255)
           aGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTAc                                   >  1:7546426/1‑51 (MQ=255)
           aGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTAc                                   >  1:3718498/1‑51 (MQ=255)
            gCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACt                                  >  1:8004074/1‑51 (MQ=255)
            gCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACt                                  >  1:8490469/1‑51 (MQ=255)
            gCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACt                                  >  1:8034748/1‑51 (MQ=255)
            gCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACt                                  <  1:4245605/51‑1 (MQ=255)
             cACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTg                                 >  1:9505645/1‑51 (MQ=255)
             cACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTg                                 <  1:8057487/51‑1 (MQ=255)
              aCTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTgg                                >  1:10898182/1‑51 (MQ=255)
                tGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGcc                              <  1:9993134/51‑1 (MQ=255)
                tGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGcc                              <  1:2793128/51‑1 (MQ=255)
                tGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGcc                              <  1:10560377/51‑1 (MQ=255)
                tGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGcc                              >  1:2317134/1‑51 (MQ=255)
                tGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGcc                              <  1:5029057/51‑1 (MQ=255)
                 ggTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCg                             >  1:5574504/1‑51 (MQ=255)
                  gTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCg                             >  1:9578387/1‑50 (MQ=255)
                  gTAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGc                            <  1:5635032/51‑1 (MQ=255)
                   tAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCa                           <  1:8003411/51‑1 (MQ=255)
                   tAGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCa                           >  1:9665507/1‑51 (MQ=255)
                    aGGTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCaa                          <  1:5631264/51‑1 (MQ=255)
                     ggTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAc                         <  1:8096912/51‑1 (MQ=255)
                      gTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAcc                        >  1:6474118/1‑51 (MQ=255)
                      gTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAcc                        >  1:7239597/1‑51 (MQ=255)
                      gTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAcc                        <  1:7196513/51‑1 (MQ=255)
                      gTTTTGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAcc                        <  1:1211948/51‑1 (MQ=255)
                       ttttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCaa                          <  1:8819123/48‑1 (MQ=255)
                       ttttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAccc                       >  1:2377045/1‑51 (MQ=255)
                       ttttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAAccc                       >  1:1505425/1‑51 (MQ=255)
                        tttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCg                      >  1:3386107/1‑51 (MQ=255)
                         ttGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGc                     >  1:2098535/1‑51 (MQ=255)
                          tGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCa                    <  1:6829674/51‑1 (MQ=255)
                          tGGTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCa                    <  1:9326454/51‑1 (MQ=255)
                            gTACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGa                  >  1:6354949/1‑51 (MQ=255)
                             tACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAc                 <  1:1885912/51‑1 (MQ=255)
                             tACTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAc                 >  1:10064354/1‑51 (MQ=255)
                              aCTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAcc                <  1:140024/51‑1 (MQ=255)
                              aCTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAcc                >  1:8783907/1‑51 (MQ=255)
                              aCTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAcc                <  1:10851329/51‑1 (MQ=255)
                               cTTTTGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAcc                <  1:7193907/50‑1 (MQ=255)
                                ttttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCgg              >  1:10474999/1‑51 (MQ=255)
                                ttttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCgg              >  1:7925532/1‑51 (MQ=255)
                                ttttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCgg              >  1:2255600/1‑51 (MQ=255)
                                ttttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCgg              >  1:4723916/1‑51 (MQ=255)
                                ttttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCgg              >  1:7453678/1‑51 (MQ=255)
                                 tttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGAcc                >  1:10819996/1‑48 (MQ=255)
                                 tttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGt             >  1:3710637/1‑51 (MQ=255)
                                  ttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTa            <  1:4059886/51‑1 (MQ=255)
                                  ttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTa            <  1:6920751/51‑1 (MQ=255)
                                  ttGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTa            >  1:6463545/1‑51 (MQ=255)
                                   tGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTaa           <  1:1138980/51‑1 (MQ=255)
                                   tGCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTaa           >  1:3767191/1‑51 (MQ=255)
                                    gCCGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTaaa          >  1:2957922/1‑51 (MQ=255)
                                     ccGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAag         >  1:6119421/1‑51 (MQ=255)
                                     ccGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAag         <  1:3069226/51‑1 (MQ=255)
                                     ccGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAag         <  1:8989177/51‑1 (MQ=255)
                                      cGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAaga        <  1:10040331/51‑1 (MQ=255)
                                      cGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAaga        >  1:9453693/1‑51 (MQ=255)
                                      cGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAaga        >  1:1319526/1‑51 (MQ=255)
                                      cGTTAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAaga        <  1:4180672/51‑1 (MQ=255)
                                        ttAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAgaga      >  1:1170244/1‑51 (MQ=255)
                                         tAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAt     <  1:4555544/51‑1 (MQ=255)
                                         tAAAGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGAt     <  1:9095470/51‑1 (MQ=255)
                                           aaGAGTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCa   >  1:1317413/1‑51 (MQ=255)
                                            agagTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAc  >  1:6700984/1‑51 (MQ=255)
                                            agagTGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAc  >  1:6771611/1‑51 (MQ=255)
                                                |                                              
GGGGATGATGTAGCACTGGTAGGTTTTGGTACTTTTGCCGTTAAAGAGCGTGCTGCCCGTACTGGCCGCAACCCGCAGACCGGTAAAGAGATCAC  >  REL606/433620‑433714

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 29 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: