Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,464,461 T→G Y31S (TAT→TCT)  envZ ← osmolarity sensor protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,464,4610TG98.5% 258.1 / NA 66Y31S (TAT→TCT) envZosmolarity sensor protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  major base G (27/37);  minor base C (0/1);  total (27/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCACTTCGTACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGATAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCAAATGAACTTCGTGGCGA  >  REL606/3464389‑3464534
                                                                        |                                                                         
cGCACTTCGTACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGaa                                                                         <  1:742478/75‑1 (MQ=255)
 gCACTTCGTACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAg                                                                        <  2:1345413/75‑1 (MQ=255)
  cACTTCGTACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAgt                                                                       >  1:428329/1‑75 (MQ=255)
   aCTTCGTACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAgtc                                                                      >  1:1380891/1‑75 (MQ=255)
        gTACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCAc                                                                 <  1:1236668/75‑1 (MQ=255)
         tACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTGGGCAcc                                                                >  1:860707/1‑75 (MQ=255)
         tACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCAcc                                                                >  2:1294138/1‑75 (MQ=255)
         tACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCAcc                                                                >  1:2114262/1‑75 (MQ=255)
          aCGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCa                                                               >  1:1716855/1‑75 (MQ=255)
             cgAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGc                                                            <  1:1440943/75‑1 (MQ=255)
              gAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCt                                                           >  1:1426845/1‑75 (MQ=255)
              gAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCt                                                           <  1:1944956/75‑1 (MQ=255)
                    tttATTAAACTGCTGGCGGCTCAGCAAACTCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCggc                                                     >  1:1166013/1‑73 (MQ=255)
                    tttATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGaa                                                     >  1:1166009/1‑75 (MQ=255)
                    tttATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGaa                                                     >  2:1850491/1‑75 (MQ=255)
                     ttattaAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAAc                                                    <  2:1651309/75‑1 (MQ=255)
                      tattaAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACa                                                   <  1:1826102/75‑1 (MQ=255)
                       attaAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGACGTCGCCACCAGGCTGGCGAACAg                                                  >  2:1563241/1‑75 (MQ=255)
                       attaAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAg                                                  <  2:1102034/75‑1 (MQ=255)
                        ttaAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCAACATGATGGCGAAAAGc                                                 >  1:443531/1‑75 (MQ=255)
                          aaaCTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCaa                                               <  1:206630/75‑1 (MQ=255)
                            aCTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAgg                                             <  2:1590132/75‑1 (MQ=255)
                              tgctgGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCCGGCTGGGGAACAGCCAgggg                                           >  1:584792/1‑73 (MQ=255)
                              tgctgGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAgggg                                           >  2:209906/1‑73 (MQ=255)
                              tgctgGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTg                                           <  2:1953865/75‑1 (MQ=255)
                               gctgGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGa                                          <  1:1869823/75‑1 (MQ=255)
                                 tgGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACg                                        >  2:1834328/1‑75 (MQ=255)
                                 tgGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAAAAGCAAGGGGACg                                        >  2:862988/1‑75 (MQ=255)
                                 tgGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCACGCTGGGGAAAAACAAGGTGACg                                        >  2:1701718/1‑75 (MQ=255)
                                  ggAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGcgacga                                       >  2:1523134/1‑75 (MQ=255)
                                  ggAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGa                                       <  2:1348556/75‑1 (MQ=255)
                                   gAGGCTCAGCAAAATCGCGAGGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGAt                                      <  2:1837770/75‑1 (MQ=255)
                                   gAGGCTAAGCAACATCGCGGCCCTCCGCAACCCGAGAGCAGTCGTCACCCGGCGGTCGAACGGCAAGGTGACGAt                                      <  2:565741/75‑1 (MQ=255)
                                    aGGCTCAGCAAAATCGCGAAGTTCAGCACCACCGGAGAAGTCGTCTCCAGGCTGGCGAACAGCAAGGTGACGATg                                     <  2:1967288/75‑1 (MQ=255)
                                    aGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATg                                     <  2:1206198/75‑1 (MQ=255)
                                    aGGCTCAGCAAAATCGCAAAGTTCAGCACCACCAGACAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATg                                     <  2:839112/75‑1 (MQ=255)
                                     ggCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGa                                    <  2:2052304/75‑1 (MQ=255)
                                      gCTCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAg                                   <  2:1114987/75‑1 (MQ=255)
                                        tCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCCAGGTGACGCTGAGCa                                 >  2:1819069/1‑75 (MQ=255)
                                        tCAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCa                                 >  1:2093861/1‑75 (MQ=255)
                                        tCAGCAAAAACGCGAAGTTCAGCCCCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCa                                 <  2:1584585/75‑1 (MQ=255)
                                         cAGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCCGGCGGAAAGCAAGGGGGCGATGAGCaa                                >  1:53414/1‑75 (MQ=255)
                                          aGCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAAt                               <  2:814943/75‑1 (MQ=255)
                                           gCAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCaaga                              <  2:95482/75‑3 (MQ=255)
                                            cAAAATCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATaa                             <  2:1802547/75‑1 (MQ=255)
                                              aaaTCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATaaag                           >  1:130425/1‑73 (MQ=255)
                                              aaaTCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACg                           <  2:1866574/75‑1 (MQ=255)
                                              aaaTCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACg                           <  2:709136/75‑1 (MQ=255)
                                              aaaTCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACg                           <  2:1597903/75‑1 (MQ=255)
                                              aaaTCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACg                           <  2:1032128/75‑1 (MQ=255)
                                               aaTCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGt                          <  2:218198/75‑1 (MQ=255)
                                                aTCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTa                         <  1:1910751/75‑1 (MQ=255)
                                                aTCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTa                         <  1:1882570/75‑1 (MQ=255)
                                                 tCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTAc                        >  2:1827494/1‑75 (MQ=255)
                                                 tCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTAc                        <  2:1771157/75‑1 (MQ=255)
                                                 tCGCGAAGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTAc                        <  2:517311/75‑1 (MQ=255)
                                                   gcgAACTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAAAAGCAAGGTGACGATGAGCAATAACGTACgg                      <  2:1229504/75‑1 (MQ=255)
                                                      aaGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCa                   >  2:552497/1‑75 (MQ=255)
                                                       aGTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAAAAAGAAGGGGACGGTGAGCAATAACGTACGGGCaa                  >  2:1730570/1‑75 (MQ=255)
                                                        gTTCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCaaa                 <  1:1064665/75‑1 (MQ=255)
                                                         tccAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCAAAt                <  2:329574/73‑1 (MQ=255)
                                                          tCAGCACCACCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGAGACGATGAGCAATAACGTACGGGCAAATg               >  2:55726/1‑75 (MQ=255)
                                                              caccacCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCAAATGAACt           <  1:2127935/75‑1 (MQ=255)
                                                              caccacCAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCAAATGAACt           <  2:1379177/75‑1 (MQ=255)
                                                               accaccAGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCAAATGAACtt          >  2:242343/1‑75 (MQ=255)
                                                                  accaGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCAAATGAACTTCGt       <  2:432631/75‑1 (MQ=255)
                                                                  accaGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCAAATGAACTTCGt       <  2:186476/75‑1 (MQ=255)
                                                                   ccaGAGAAGTCGTCACCAGGCTGGCGAACAGCAAGGGGACGATGAGCAATAACGTACGGGCAAATGAACTTCGTg      >  1:1392011/1‑75 (MQ=255)
                                                                       aTAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCAAATGAACTTCGTGGCga  <  2:538948/75‑1 (MQ=255)
                                                                        |                                                                         
CGCACTTCGTACGCGAGGACTTTATTAAACTGCTGGAGGCTCAGCAAAATCGCGAAGTTCAGCACCACCAGATAAGTCGTCACCAGGCTGGCGAACAGCAAGGTGACGATGAGCAATAACGTACGGGCAAATGAACTTCGTGGCGA  >  REL606/3464389‑3464534

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