Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 2,655,307 T→C N92S (AAT→AGT)  clpB ← protein disaggregation chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6062,655,3070TC99.5% 923.6 / ‑4.4 216N92S (AAT→AGT) clpBprotein disaggregation chaperone
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (103/112);  total (104/112)
Fisher's exact test for biased strand distribution p-value = 4.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

Only 100 of 219 total aligned reads displayed.

AAGTGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTA  >  REL606/2655234‑2655379
                                                                         |                                                                        
 aGTGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTaa                                                                       <  2:6960254/76‑1 (MQ=255)
 aGTGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTaa                                                                       <  2:3537029/76‑1 (MQ=255)
  gTGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAg                                                                      <  2:5670537/76‑1 (MQ=255)
  gTGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAg                                                                      <  1:3869646/76‑1 (MQ=255)
  gTGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAg                                                                      <  2:6354529/76‑1 (MQ=255)
  gTGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAg                                                                      >  1:2372204/1‑76 (MQ=255)
   tGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGa                                                                     >  2:3815922/1‑76 (MQ=255)
   tGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGa                                                                     <  2:731857/76‑1 (MQ=255)
   tGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGa                                                                     <  2:5363374/76‑1 (MQ=255)
   tGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGa                                                                     <  1:2852722/76‑1 (MQ=255)
     ccgccAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAAc                                                                   <  1:6272539/76‑1 (MQ=255)
      cgccAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAAcg                                                                  <  2:3112523/76‑1 (MQ=255)
      cgccAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAAcg                                                                  <  2:6558649/76‑1 (MQ=255)
         cAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACgcgc                                                               >  2:6052700/1‑76 (MQ=255)
         cAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACgcgc                                                               >  2:5034465/1‑76 (MQ=255)
          aGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCa                                                              >  1:3750143/1‑76 (MQ=255)
          aGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCa                                                              >  1:5281050/1‑76 (MQ=255)
             acgaacAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCa                                                           <  2:1987421/76‑1 (MQ=255)
               gaacAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGa                                                         <  2:5773651/76‑1 (MQ=255)
               gaacAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGa                                                         <  1:5034504/76‑1 (MQ=255)
                aacAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGAt                                                        >  1:655830/1‑76 (MQ=255)
                aacAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGAt                                                        <  1:4476539/76‑1 (MQ=255)
                 acAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATc                                                       <  2:1440273/76‑1 (MQ=255)
                   aGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCt                                                     >  1:6133444/1‑76 (MQ=255)
                   aGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCt                                                     <  2:4939674/76‑1 (MQ=255)
                     ttCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCtgt                                                   >  1:240369/1‑76 (MQ=255)
                     ttCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCtgt                                                   >  2:8677/1‑76 (MQ=255)
                     ttCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCtgt                                                   >  1:3436886/1‑76 (MQ=255)
                      tCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCtgtg                                                  >  2:4060083/1‑76 (MQ=255)
                        tGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGAt                                                <  1:4223217/76‑1 (MQ=255)
                        tGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGAt                                                >  2:3348885/1‑76 (MQ=255)
                         gACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATg                                               <  1:35971/76‑1 (MQ=255)
                         gACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATg                                               <  2:4263403/76‑1 (MQ=255)
                         gACGAGATAAAGTTAGCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATg                                               <  1:1571791/76‑1 (MQ=255)
                          aCGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATgg                                              >  2:4213984/1‑76 (MQ=255)
                           agagaTAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGc                                             <  1:2024768/75‑1 (MQ=255)
                            gagaTAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCt                                            >  2:2291224/1‑76 (MQ=255)
                             agaTAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTg                                           >  2:539979/1‑76 (MQ=255)
                              gaTAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTgg                                          >  2:3709570/1‑76 (MQ=255)
                               aTAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGa                                         <  1:2911803/76‑1 (MQ=255)
                                tAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGAc                                        >  1:5244605/1‑76 (MQ=255)
                                 aaaGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACa                                       <  1:100496/76‑1 (MQ=255)
                                 aaaGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACa                                       >  2:6439797/1‑76 (MQ=255)
                                  aaGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACAt                                      >  1:6856914/1‑76 (MQ=255)
                                  aaGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACAt                                      <  2:1740672/76‑1 (MQ=255)
                                    gTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATca                                    >  2:515620/1‑76 (MQ=255)
                                    gTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATca                                    <  1:4289060/76‑1 (MQ=255)
                                    gTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATca                                    <  2:2859608/76‑1 (MQ=255)
                                     ttATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcac                                   >  2:4807634/1‑76 (MQ=255)
                                     ttATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcac                                   >  2:4455140/1‑76 (MQ=255)
                                     ttATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcac                                   >  1:6731784/1‑76 (MQ=255)
                                      tATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcacc                                  <  1:5902740/76‑1 (MQ=255)
                                      tATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcacc                                  >  2:6736714/1‑76 (MQ=255)
                                      tATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcacc                                  >  2:3020666/1‑76 (MQ=255)
                                       aTCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcacca                                 <  1:6009151/76‑1 (MQ=255)
                                       aTCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcacca                                 <  1:1700201/76‑1 (MQ=255)
                                        tCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcaccac                                <  2:3981651/76‑1 (MQ=255)
                                        tCACCACGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATcaccac                                >  2:5273325/1‑76 (MQ=255)
                                         caccacGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATTaccacc                               >  1:5752102/1‑76 (MQ=255)
                                         caccacGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCaccacc                               >  1:1110882/1‑76 (MQ=255)
                                         caccacGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCaccacc                               >  1:4621287/1‑76 (MQ=255)
                                         caccacGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCaccacc                               >  2:1261229/1‑76 (MQ=255)
                                            cacGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGt                            <  2:3442164/76‑1 (MQ=255)
                                            cacGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGt                            >  1:5583998/1‑76 (MQ=255)
                                              cGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTAc                          <  2:4538620/76‑1 (MQ=255)
                                              cGTTTTTGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTAc                          <  1:730261/76‑1 (MQ=255)
                                                 ttttGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCtt                       <  2:5737826/76‑1 (MQ=255)
                                                  tttGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTc                      >  2:2561836/1‑76 (MQ=255)
                                                  tttGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTc                      >  2:5865186/1‑76 (MQ=255)
                                                  tttGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTc                      >  1:308755/1‑76 (MQ=255)
                                                   ttGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCa                     >  1:2729537/1‑76 (MQ=255)
                                                   ttGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCGGGACATCACCACCAGTACCTTCa                     <  2:723031/76‑1 (MQ=255)
                                                    tGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCaa                    <  2:5426156/76‑1 (MQ=255)
                                                    tGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCaa                    <  1:276494/76‑1 (MQ=255)
                                                    tGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCaa                    >  1:982007/1‑76 (MQ=255)
                                                    tGCGCCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCaa                    <  2:1776280/76‑1 (MQ=255)
                                                      cgcCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAAcc                  >  1:1289268/1‑76 (MQ=255)
                                                       gcCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCt                 <  2:1363154/76‑1 (MQ=255)
                                                       gcCAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCt                 <  1:3709396/76‑1 (MQ=255)
                                                         cAGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGc               <  1:335631/76‑1 (MQ=255)
                                                          aGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCg              >  2:6838256/1‑76 (MQ=255)
                                                          aGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCg              <  1:1174486/76‑1 (MQ=255)
                                                          aGCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCg              >  2:2257275/1‑76 (MQ=255)
                                                           gCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCgg             >  2:169087/1‑76 (MQ=255)
                                                           gCTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCgg             >  2:2166458/1‑76 (MQ=255)
                                                            cTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGt            <  2:2783698/76‑1 (MQ=255)
                                                            cTTGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGt            >  1:4712093/1‑76 (MQ=255)
                                                             ttGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTa           >  1:3965503/1‑76 (MQ=255)
                                                              tGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTaa          <  2:3197605/76‑1 (MQ=255)
                                                              tGTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTaa          <  1:2140564/76‑1 (MQ=255)
                                                               gTCGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTaaa         <  2:5990938/76‑1 (MQ=255)
                                                                 cGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACg       >  1:4081125/1‑76 (MQ=255)
                                                                 cGCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACg       >  1:2241646/1‑76 (MQ=255)
                                                                  gCAAAGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGa      <  1:4633463/76‑1 (MQ=255)
                                                                    aaaGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGAtt    <  1:5110301/76‑1 (MQ=255)
                                                                    aaaGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGCAAACGAtt    >  2:5590673/1‑76 (MQ=255)
                                                                     aaGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGAttt   <  1:540650/76‑1 (MQ=255)
                                                                     aaGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGAttt   >  2:5537186/1‑76 (MQ=255)
                                                                      aGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCGCCACCAGTACCTTCAACCTGCGGTAAACGATTTa  <  2:6704097/76‑1 (MQ=255)
                                                                      aGACTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTa  <  1:78943/76‑1 (MQ=255)
                                                                         |                                                                        
AAGTGCCGCCAGAACGAACAGTTCTGACGAGATAAAGTTATCACCACGTTTTTGCGCCAGCTTGTCGCAAAGATTAAGAACGCGCACCAGATCCTGTGATGGCTGGACATCACCACCAGTACCTTCAACCTGCGGTAAACGATTTA  >  REL606/2655234‑2655379

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: