Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 1,588,754 A→C F50L (TTT→TTG eamA ← cysteine and O‑acetyl‑L‑serine efflux system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL6061,588,7540AC97.0% 261.8 / ‑1.9 67F50L (TTT→TTGeamAcysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (31/34);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.63e-01

AAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTATGAAGC  >  REL606/1588681‑1588826
                                                                         |                                                                        
aagCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACcaa                                                                        >  2:213602/1‑76 (MQ=255)
 agCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACcaaa                                                                       <  1:263070/76‑1 (MQ=255)
  gCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACcaaaa                                                                      <  2:892673/76‑1 (MQ=255)
   cAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACcaaaaa                                                                     >  1:472902/1‑76 (MQ=255)
    aaaCTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAg                                                                    >  2:170602/1‑76 (MQ=255)
     aaCTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGa                                                                   >  1:359544/1‑76 (MQ=255)
      aCTGCGCAAAACTGATGGTTAATCCATACCCCGGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGAt                                                                  <  2:239372/76‑1 (MQ=255)
      aCTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGAt                                                                  <  1:1548183/76‑1 (MQ=255)
      aCTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGAt                                                                  <  2:836925/76‑1 (MQ=255)
       cTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATa                                                                 >  2:20865/1‑76 (MQ=255)
        tGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAg                                                                <  2:471717/76‑1 (MQ=255)
        tGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAg                                                                <  2:1786090/76‑1 (MQ=255)
         gcgcAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGc                                                               <  1:119574/76‑1 (MQ=255)
             aaaaCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGa                                                           >  1:937818/1‑76 (MQ=255)
             aaaaCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGa                                                           >  1:1884089/1‑76 (MQ=255)
              aaaCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGaa                                                          <  2:525639/76‑1 (MQ=255)
              aaaCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGaa                                                          >  2:1803795/1‑76 (MQ=255)
               aaCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGaaa                                                         >  2:1239211/1‑76 (MQ=255)
               aaCTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGaaa                                                         <  2:775846/76‑1 (MQ=255)
                 cTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAg                                                       >  1:450637/1‑76 (MQ=255)
                   gATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCg                                                     >  1:674986/1‑76 (MQ=255)
                    aTGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGa                                                    >  2:1338722/1‑76 (MQ=255)
                     tGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGAc                                                   <  2:783732/76‑1 (MQ=255)
                      ggTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGAcc                                                  >  2:1043908/1‑76 (MQ=255)
                      ggTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGAcc                                                  >  2:498119/1‑76 (MQ=255)
                       gTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCa                                                 <  1:466830/76‑1 (MQ=255)
                        ttAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCaa                                                <  1:10820/76‑1 (MQ=255)
                        ttAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCaa                                                <  2:964462/76‑1 (MQ=255)
                        ttAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCaa                                                >  1:1743617/1‑76 (MQ=255)
                         tAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAAc                                               >  1:1806111/1‑76 (MQ=255)
                         tAACCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAAc                                               >  1:1875320/1‑76 (MQ=255)
                          aaTCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACa                                              <  2:1882983/76‑1 (MQ=255)
                           aTCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACAt                                             >  1:1575207/1‑76 (MQ=255)
                           aTCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAACAt                                             <  1:1073708/76‑1 (MQ=255)
                            tCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATa                                            >  2:822003/1‑76 (MQ=255)
                             ccATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATaa                                           <  2:178883/76‑1 (MQ=255)
                              cATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATaaa                                          >  1:1760572/1‑76 (MQ=255)
                              cATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATaaa                                          >  2:1547888/1‑76 (MQ=255)
                               aTACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAg                                         <  2:611709/76‑1 (MQ=255)
                                 aCCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAgcg                                       <  2:1905735/76‑1 (MQ=255)
                                 aCCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAgcg                                       >  1:1917546/1‑76 (MQ=255)
                                     cagcagCAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCaaa                                   >  1:1468134/1‑76 (MQ=255)
                                      agcagcAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAAc                                  <  2:1292711/76‑1 (MQ=255)
                                      agcagcAAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAAc                                  >  2:879087/1‑76 (MQ=255)
                                       gcagcaAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAAcc                                 >  2:1819741/1‑76 (MQ=255)
                                        cagcaAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCg                                <  2:1397031/76‑1 (MQ=255)
                                          gcaAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCgac                              >  2:1581324/1‑74 (MQ=255)
                                           caAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGcc                             <  1:1041164/76‑1 (MQ=255)
                                           caAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGcc                             <  1:1343598/76‑1 (MQ=255)
                                           caAATTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGcc                             <  1:711718/76‑1 (MQ=255)
                                             aaTTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAg                           <  2:411182/76‑1 (MQ=255)
                                             aaTTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAg                           >  1:1042641/1‑76 (MQ=255)
                                             aaTTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAg                           <  1:1337580/76‑1 (MQ=255)
                                              aTTCAGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGc                          >  1:693426/1‑76 (MQ=255)
                                                  aGTGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTa                      <  2:488638/76‑1 (MQ=255)
                                                    tGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGc                    >  1:1291932/1‑76 (MQ=255)
                                                    tGGTACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGc                    >  1:1329446/1‑76 (MQ=255)
                                                       tACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCggg                 >  2:386309/1‑75 (MQ=255)
                                                       tACTTTCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGt                 <  1:622647/76‑1 (MQ=255)
                                                           ttCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCa             >  2:548001/1‑76 (MQ=255)
                                                           ttCGGTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCa             >  2:1942009/1‑76 (MQ=255)
                                                              ggTCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGt          <  1:950541/76‑1 (MQ=255)
                                                                tCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTa        <  2:1036576/76‑1 (MQ=255)
                                                                tCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTa        <  2:1254282/76‑1 (MQ=255)
                                                                tCGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTa        <  1:553643/76‑1 (MQ=255)
                                                                 cGTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTAt       <  2:1016287/76‑1 (MQ=255)
                                                                  gTGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTATg      <  2:278937/76‑1 (MQ=255)
                                                                   tGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTATGa     >  2:1098788/1‑76 (MQ=255)
                                                                   tGCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTATGa     <  2:1340147/76‑1 (MQ=255)
                                                                    gCGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTATGaa    <  1:1687428/76‑1 (MQ=255)
                                                                     cGACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTATGAAg   <  2:1968392/76‑1 (MQ=255)
                                                                      gACCAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTATGAAGc  <  2:832729/76‑1 (MQ=255)
                                                                         |                                                                        
AAGCAAACTGCGCAAAACTGATGGTTAATCCATACCCCAGCAGCAAATTCAGTGGTACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAACATAAAGCGCAAACCGGCCAGCATTAGCGGTGGCATGTTATGAAGC  >  REL606/1588681‑1588826

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: