Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,465,075 G→A S65L (TCG→TTG)  ompR ← osmolarity response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,465,0750GA97.9% 175.0 / NA 47S65L (TCG→TTG) ompRosmolarity response regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (27/19);  minor base T (1/0);  total (28/19)
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

ACTTCTTCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCGACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGA  >  REL606/3465003‑3465146
                                                                        |                                                                       
aCTTCTTCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACa                                                                      <  1:1150351/76‑1 (MQ=255)
 cttcttCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATcaacaa                                                                     >  2:845334/1‑76 (MQ=255)
 cttcttCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATcaacaa                                                                     <  1:804126/76‑1 (MQ=255)
 cttcttCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATcaacaa                                                                     <  1:622284/76‑1 (MQ=255)
  ttcttcCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAg                                                                    >  2:569720/1‑76 (MQ=255)
    cttcCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGcc                                                                  >  1:913959/1‑76 (MQ=255)
    cttcCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGcc                                                                  >  2:1105900/1‑76 (MQ=255)
        cccTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCt                                                              <  1:814456/76‑1 (MQ=255)
        cccTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCt                                                              <  1:964483/76‑1 (MQ=255)
          cTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTc                                                            >  2:735932/1‑76 (MQ=255)
           ttgcGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCa                                                           <  1:108807/73‑1 (MQ=255)
             tCGCCGTCNCCATAATGATCGGCATCGGGTTGCTCTGACCACGAAGACGTCGGCAAATCAACAAGCCATCTTCAcc                                                         <  1:1249692/76‑1 (MQ=255)
             tCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCAcc                                                         <  1:199975/76‑1 (MQ=255)
              cGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCTATCTTCACCa                                                        <  1:1318105/76‑1 (MQ=255)
               gCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAg                                                       <  1:367810/76‑1 (MQ=255)
                  gTCACCATAATGATCGGCATCGGGTTGCTTTGACTATGAAGACGTCGGCAAATCAACNAGCCATCTTCACCAGGTa                                                    <  1:1354247/76‑1 (MQ=255)
                   tCACCATAATGATCGGCATNGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTaa                                                   >  1:1335837/1‑76 (MQ=255)
                      ccATAATGATCGGCATCGGGTTGNTCTGANTACGAAGACGTCGGNAATTCATCAAGCCATCTTCACCAGGTAACAt                                                >  2:183831/1‑76 (MQ=255)
                        aTAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTa                                              >  2:154155/1‑76 (MQ=255)
                        aTAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTa                                              <  1:687112/76‑1 (MQ=255)
                         tAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTaa                                             >  2:502895/1‑76 (MQ=255)
                             gATCGGCATCGGGTTGCTCTGACTACGAAGCCGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATcc                                         >  2:271189/1‑76 (MQ=255)
                              gtCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCa                                        <  2:96224/75‑1 (MQ=255)
                               tCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAg                                       <  2:627837/76‑1 (MQ=255)
                               tCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAg                                       >  2:159706/1‑76 (MQ=255)
                                cGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGt                                      <  2:1042022/76‑1 (MQ=255)
                                  gCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTAc                                    <  1:987725/76‑1 (MQ=255)
                                  gCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTAc                                    <  1:374026/76‑1 (MQ=255)
                                    aTCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCa                                  <  1:1191297/76‑1 (MQ=255)
                                     tCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCAt                                 >  1:701109/1‑76 (MQ=255)
                                      cGGGTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATa                                >  2:773191/1‑76 (MQ=255)
                                        ggTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAg                              <  2:376404/76‑1 (MQ=255)
                                        ggTTGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATCAAATCCAGTACCATAAg                              >  1:161460/1‑76 (MQ=255)
                                          ttGCTCTGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAAACCAGTACCATAAGAt                            >  2:42888/1‑76 (MQ=255)
                                             ctctGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGa                         >  1:819265/1‑76 (MQ=255)
                                              tctGAATACGGAGACATCGGCAAAGCANCAAGCCATCTTCGCCATGTANCATCGNATCGAGTACCATAAGATGGaa                        <  1:1384182/76‑1 (MQ=255)
                                                tGACTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAg                      >  1:1120879/1‑76 (MQ=255)
                                                   cTACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGAtt                   >  2:1086058/1‑76 (MQ=255)
                                                    tACGNAGACGTCGGCAAATCAACCAGCCATCCTCACCAGNTAACATTAAATCCAGTACCATAAGATGGAAAGATTc                  <  1:171870/76‑1 (MQ=255)
                                                    tACGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTc                  >  1:703688/1‑76 (MQ=255)
                                                      cGAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCAc                >  2:1134388/1‑76 (MQ=255)
                                                       gAAGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACg               >  1:676719/1‑76 (MQ=255)
                                                        aaGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGa              >  2:175196/1‑76 (MQ=255)
                                                        aaGACGTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGa              >  2:870464/1‑76 (MQ=255)
                                                             gTCGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTcag         >  2:235185/1‑76 (MQ=255)
                                                               cGGCAAATCTACAAGCCATCTTCACCAGGAAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTcagca       >  2:1216156/1‑76 (MQ=255)
                                                               cGGCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTcagca       >  2:21788/1‑76 (MQ=255)
                                                                ggCAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTcagcag      >  1:1202500/1‑76 (MQ=255)
                                                                  cAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGc    >  2:137745/1‑76 (MQ=255)
                                                                  cAAATCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGc    >  2:1162142/1‑76 (MQ=255)
                                                                    aaTCAACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGa  <  1:1269430/76‑1 (MQ=255)
                                                                        |                                                                       
ACTTCTTCCCCTTTCGCCGTCACCATAATGATCGGCATCGGGTTGCTCTGACTACGAAGACGTCGGCAAATCGACAAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGA  >  REL606/3465003‑3465146

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