Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 3,465,147 T→A D41V (GAT→GTT)  ompR ← osmolarity response regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6063,465,1470TA96.8% 245.8 / ‑3.8 62D41V (GAT→GTT) ompRosmolarity response regulator
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base A (31/29);  total (32/30)
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

AAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGATCCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTCCAGCAGCG  >  REL606/3465078‑3465216
                                                                     |                                                                     
aaGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATc                                                                  <  2:1521750/75‑1 (MQ=255)
   ccATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGt                                                               <  1:798867/75‑1 (MQ=255)
     aTCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTc                                                             <  2:623493/75‑1 (MQ=255)
      tCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCt                                                            >  1:1148648/1‑75 (MQ=255)
      tCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCt                                                            <  1:768622/75‑1 (MQ=255)
      tCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCt                                                            >  1:1397297/1‑75 (MQ=255)
        ttCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGc                                                          >  2:1078926/1‑75 (MQ=255)
           aCCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCAtt                                                       <  1:8771/75‑1 (MQ=255)
            ccAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTa                                                      >  2:1369246/1‑75 (MQ=255)
            ccAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTa                                                      >  2:264637/1‑75 (MQ=255)
             cAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAg                                                     <  1:891834/75‑1 (MQ=255)
              aGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGATCCATCTGTTCTGCATTAGc                                                    >  2:1538166/1‑75 (MQ=255)
              aGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGc                                                    >  1:520684/1‑75 (MQ=255)
              aGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGc                                                    >  2:703547/1‑75 (MQ=255)
                 tAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGAc                                                 <  1:819755/75‑1 (MQ=255)
                   aCATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGc                                               >  1:512467/1‑75 (MQ=255)
                     aTTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCtt                                             >  2:333652/1‑75 (MQ=255)
                     aTTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCtt                                             >  2:333663/1‑75 (MQ=255)
                      ttAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTc                                            <  2:1148648/75‑1 (MQ=255)
                      ttAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTc                                            >  1:498381/1‑75 (MQ=255)
                        aaaTCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGa                                          <  1:134213/75‑1 (MQ=255)
                        aaaTCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGa                                          <  2:876765/75‑1 (MQ=255)
                         aaTCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGaa                                         >  1:971047/1‑75 (MQ=255)
                         aaTCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGaa                                         <  1:670148/75‑1 (MQ=255)
                         aaTCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGaa                                         <  1:1086256/75‑1 (MQ=255)
                          aTCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAAc                                        >  1:1096123/1‑75 (MQ=255)
                           tCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAAcc                                       <  1:1185947/75‑1 (MQ=255)
                            ccACTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCt                                      <  2:1461423/75‑1 (MQ=255)
                              aGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTgg                                    <  2:1379714/75‑1 (MQ=255)
                               gTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGa                                   >  1:1500382/1‑75 (MQ=255)
                                tcccAGAAGATGGAAAGATTCACGAGTCAGCAGGCGATCCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGaa                                  <  2:614412/73‑1 (MQ=255)
                                tACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGaa                                  <  1:1371517/75‑1 (MQ=255)
                                  ccATAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGc                                >  2:485601/1‑75 (MQ=255)
                                    aTAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCt                              >  1:794727/1‑75 (MQ=255)
                                    aTAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCt                              >  2:597891/1‑75 (MQ=255)
                                     tAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCtt                             <  2:819016/75‑1 (MQ=255)
                                     tAAGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCtt                             >  1:645430/1‑75 (MQ=255)
                                       aGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGt                           <  1:840668/75‑1 (MQ=255)
                                       aGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGt                           <  2:222817/75‑1 (MQ=255)
                                       aGATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGt                           >  2:1510542/1‑75 (MQ=255)
                                        gATGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGtt                          <  1:600698/75‑1 (MQ=255)
                                          tGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCg                        <  1:1403882/75‑1 (MQ=255)
                                          tGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCg                        >  1:571069/1‑75 (MQ=255)
                                          tGGAAAGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCg                        >  2:1599918/1‑75 (MQ=255)
                                               aGATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTgag                   >  1:1297155/1‑75 (MQ=255)
                                                gATTCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTgaga                  >  1:992437/1‑75 (MQ=255)
                                                  ttCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATa                <  2:964481/75‑1 (MQ=255)
                                                  ttCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATa                >  1:788385/1‑75 (MQ=255)
                                                  ttCACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATa                <  2:1186051/75‑1 (MQ=255)
                                                  tgcACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATa                <  1:166679/73‑1 (MQ=255)
                                                    cACGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAAc              >  1:39804/1‑75 (MQ=255)
                                                      cGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGt            >  1:1587630/1‑75 (MQ=255)
                                                      cGAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGt            >  2:1311927/1‑75 (MQ=255)
                                                       gAGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGtt           >  1:536027/1‑75 (MQ=255)
                                                        aGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTc          >  2:1501705/1‑75 (MQ=255)
                                                        aGTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTc          <  1:964751/75‑1 (MQ=255)
                                                         gTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTcc         >  1:171397/1‑75 (MQ=255)
                                                         gTCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTcc         <  2:1110584/75‑1 (MQ=255)
                                                          tCAGCAGGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTCca        <  2:498487/75‑1 (MQ=255)
                                                             gcagGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTCcagca     >  1:43462/1‑75 (MQ=255)
                                                              cagGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTCcagcag    <  1:1458672/75‑1 (MQ=255)
                                                              cagGCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTCcagcag    <  1:787070/75‑1 (MQ=255)
                                                                ggCGAACCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTCCAGCAgcg  <  2:1285978/75‑1 (MQ=255)
                                                                     |                                                                     
AAGCCATCTTCACCAGGTAACATTAAATCCAGTACCATAAGATGGAAAGATTCACGAGTCAGCAGGCGATCCATCTGTTCTGCATTAGCGACGCTTCGAACCTGGAAGCCTTGTTCGGTGAGATAACGTTCCAGCAGCG  >  REL606/3465078‑3465216

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