Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,928,260 T→A I15N (ATC→AAC)  rho → transcription termination factor Rho

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,928,2600TA99.3% 1174.4 / ‑0.4 287I15N (ATC→AAC) rhotranscription termination factor Rho
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (124/160);  total (126/161)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.18e-01

Only 100 of 293 total aligned reads displayed.

TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928187‑3928332
                                                                         |                                                                        
 ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac                                                                       >  1:7681952/1‑76 (MQ=255)
 ttACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAacac                                                                       <  2:8958661/76‑1 (MQ=255)
  tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCAct                                                                      >  1:1078512/1‑76 (MQ=255)
  tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACAct                                                                      <  1:115677/76‑1 (MQ=255)
  tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACAct                                                                      <  2:5810975/76‑1 (MQ=255)
  tACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACAct                                                                      >  2:291544/1‑76 (MQ=255)
   aCCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActc                                                                     >  1:8620503/1‑76 (MQ=255)
    ccccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACActct                                                                    <  2:1228093/76‑1 (MQ=255)
     cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                   >  2:8805919/1‑76 (MQ=255)
     cccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACtctc                                                                   <  1:4127863/76‑1 (MQ=255)
      ccAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCg                                                                  <  2:4555178/76‑1 (MQ=255)
       cAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCgg                                                                 >  2:2051670/1‑76 (MQ=255)
        aaGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGc                                                                <  2:2377813/76‑1 (MQ=255)
         aGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCg                                                               <  2:7834261/76‑1 (MQ=255)
          gTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGa                                                              <  2:2133514/76‑1 (MQ=255)
           tttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGATTCTGAGCTGAACACTCTCGGCGaa                                                             >  2:2697761/1‑76 (MQ=255)
            ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaa                                                            >  1:132287/1‑76 (MQ=255)
            ttAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaa                                                            >  1:5874971/1‑76 (MQ=255)
             tAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                           >  1:5136780/1‑76 (MQ=255)
              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAGAat                                                          >  2:8122673/1‑76 (MQ=255)
              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          <  1:2743901/76‑1 (MQ=255)
              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          >  2:1658942/1‑76 (MQ=255)
              aaGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAat                                                          >  2:6576046/1‑76 (MQ=255)
               aGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAata                                                         >  2:4578236/1‑76 (MQ=255)
                  aCTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATgg                                                      >  1:3768912/1‑76 (MQ=255)
                   cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATggg                                                     >  2:6174922/1‑76 (MQ=255)
                   cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATggg                                                     >  1:7865781/1‑76 (MQ=255)
                   cTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATggg                                                     <  1:7530749/76‑1 (MQ=255)
                       caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                 <  2:8999709/76‑1 (MQ=255)
                       caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                 >  1:6903070/1‑76 (MQ=255)
                       caccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTg                                                 >  1:3210865/1‑76 (MQ=255)
                        accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                                >  2:224338/1‑76 (MQ=255)
                        accacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTgg                                                <  2:5943894/76‑1 (MQ=255)
                         ccacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGa                                               >  2:7310032/1‑76 (MQ=255)
                          cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa                                              >  2:2797234/1‑76 (MQ=255)
                          cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa                                              >  1:5367678/1‑76 (MQ=255)
                          cacTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaa                                              >  1:8997512/1‑76 (MQ=255)
                           acTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGaaa                                             <  2:6996876/76‑1 (MQ=255)
                             tATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGaaaa                                                           <  1:7354080/60‑1 (MQ=255)
                              aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                          >  1:5705646/1‑76 (MQ=255)
                              aTGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAAcc                                          >  1:5526541/1‑76 (MQ=255)
                               tGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCt                                         <  1:8499953/76‑1 (MQ=255)
                               tGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCt                                         >  1:4551921/1‑76 (MQ=255)
                               tGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCt                                         <  1:7077068/76‑1 (MQ=255)
                                 aaTCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTgg                                       >  1:1723951/1‑76 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    >  2:8084532/1‑76 (MQ=255)
                                    cTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTc                                    >  1:3859641/1‑76 (MQ=255)
                                     ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTcc                                   >  2:553781/1‑75 (MQ=255)
                                     ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg                                   <  1:2141979/76‑1 (MQ=255)
                                     ttACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCg                                   >  2:8672293/1‑76 (MQ=255)
                                      tACCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGt                                  <  1:5465474/76‑1 (MQ=255)
                                       aCCGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTa                                 <  1:8229666/76‑1 (MQ=255)
                                        ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                <  2:7633200/76‑1 (MQ=255)
                                        ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                <  2:232353/76‑1 (MQ=255)
                                        ccGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTAt                                <  2:6469884/76‑1 (MQ=255)
                                         cGAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCCTATg                               >  1:4861300/1‑76 (MQ=255)
                                          gAATTAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGc                              <  1:8805919/76‑1 (MQ=255)
                                             ttAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTa                           <  2:3700806/76‑1 (MQ=255)
                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          <  2:4171560/76‑1 (MQ=255)
                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          <  2:7114217/76‑1 (MQ=255)
                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTaa                          <  2:7199833/76‑1 (MQ=255)
                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGCATGCGTaa                          >  1:1534152/1‑76 (MQ=255)
                                              tAAAGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCCTATGCGGaa                          >  1:4373231/1‑76 (MQ=255)
                                               aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTCTGCGTAAg                         >  2:8599564/1‑76 (MQ=255)
                                               aaaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAg                         <  2:3845147/76‑1 (MQ=255)
                                                aaGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGc                        <  2:541709/76‑1 (MQ=255)
                                                 aGAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCa                       <  1:3590076/76‑1 (MQ=255)
                                                  gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                      <  2:6838675/76‑1 (MQ=255)
                                                  gAATACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAg                      <  2:305528/76‑1 (MQ=255)
                                                   aaTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGGAAGCAgg                     >  2:6765589/1‑76 (MQ=255)
                                                    aTACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGa                    <  1:4698855/76‑1 (MQ=255)
                                                     tACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGAc                   <  1:8142069/76‑1 (MQ=255)
                                                     tACGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGGAAGCAGGAc                   >  1:8935468/1‑76 (MQ=255)
                                                      aCGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACa                  <  2:5021675/76‑1 (MQ=255)
                                                       cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACat                 <  2:7672844/76‑1 (MQ=255)
                                                       cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGCATGCGTAAGCCGGACat                 >  1:6374694/1‑76 (MQ=255)
                                                       cGCCGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCCTATGCGTAAGCAGGGCat                 >  2:1323959/1‑76 (MQ=255)
                                                         ccGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACatta               <  1:7563803/76‑1 (MQ=255)
                                                          cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              <  1:398926/76‑1 (MQ=255)
                                                          cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              <  2:5057412/76‑1 (MQ=255)
                                                          cGGTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattat              <  1:1906529/76‑1 (MQ=255)
                                                           ggTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt             <  2:3440208/76‑1 (MQ=255)
                                                           ggTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACattatt             <  2:1149829/76‑1 (MQ=255)
                                                            gTTTCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCCCATGCGTAAACAGGGCAATAttt            >  1:3119596/1‑76 (MQ=255)
                                                             tttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAtttt           <  1:3478912/76‑1 (MQ=255)
                                                              ttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTCAGCAGGACATTAttttt          >  1:6774002/1‑76 (MQ=255)
                                                              ttCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTAttttt          <  1:5239809/76‑1 (MQ=255)
                                                               tCTGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGGAGGGAATTATTTTTg         >  2:1517785/1‑76 (MQ=255)
                                                                 tGAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGGGTAAGGAGGAAATTATTTTTTcc       >  1:1633277/1‑76 (MQ=255)
                                                                  gAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGGCCTTATTTTTGCCa      >  1:3370575/1‑76 (MQ=255)
                                                                  gAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa      <  1:8089361/76‑1 (MQ=255)
                                                                  gAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCa      <  1:2346996/76‑1 (MQ=255)
                                                                  gAGCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCCTGTGCGGAAGCCGGACATTATTTTTGTCa      >  2:3967259/1‑76 (MQ=255)
                                                                    gCTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTTCCATc    >  2:4781053/1‑76 (MQ=255)
                                                                     cTGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATcc   >  2:3595980/1‑76 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  <  1:2559140/76‑1 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  <  1:3083998/76‑1 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  <  2:8476054/76‑1 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  <  1:4300442/76‑1 (MQ=255)
                                                                      tGAACACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCt  <  1:5952232/76‑1 (MQ=255)
                                                                         |                                                                        
TTTACCCCAAGTTTAAGAACTCACACCACTATGAATCTTACCGAATTAAAGAATACGCCGGTTTCTGAGCTGATCACTCTCGGCGAAAATATGGGGCTGGAAAACCTGGCTCGTATGCGTAAGCAGGACATTATTTTTGCCATCCT  >  REL606/3928187‑3928332

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