Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 501,525 A→G 18.2% L313P (CTG→CCG)  ybaL ← predicted transporter with NAD(P)‑binding Rossmann‑fold domain

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb501,5250AG18.2% 123.2 / 21.1 55L313P (CTG→CCG) ybaLpredicted transporter with NAD(P)‑binding Rossmann‑fold domain
Reads supporting (aligned to +/- strand):  ref base A (25/20);  new base G (6/4);  total (31/24)
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

TAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCC  >  W3110S.gb/501458‑501595
                                                                   |                                                                      
tAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGa                                                                      <  1:1847786/70‑1 (MQ=255)
  aTGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGaata                                                                   <  1:3155437/71‑1 (MQ=255)
     ccGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGaataat                                                                 >  1:4409878/1‑70 (MQ=255)
      cGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCg                                                               <  1:1841647/71‑1 (MQ=255)
       gTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGc                                                              <  1:3465102/71‑1 (MQ=255)
          cGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAg                                                           >  1:2193863/1‑71 (MQ=255)
               ggAGTTACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTc                                                       <  1:1911050/70‑1 (MQ=255)
                  gTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGc                                                     >  1:968638/1‑69 (MQ=255)
                   tGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAg                                                  <  1:902021/71‑1 (MQ=255)
                    gACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAg                                                  <  1:1143768/70‑1 (MQ=255)
                     actaaaCAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGc                                                 >  1:33958/4‑70 (MQ=255)
                     aCCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGaataa                                                                  >  1:2789788/1‑53 (MQ=255)
                     aCCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGaa                                                                     >  1:4165248/1‑50 (MQ=255)
                     aCCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCa                                                            >  1:4632620/1‑59 (MQ=255)
                     aCCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAg                                                  <  1:2232603/69‑1 (MQ=255)
                     aCCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCa                                                >  1:1910612/1‑71 (MQ=255)
                     aCCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCa                                                >  1:2239923/1‑71 (MQ=255)
                      ccAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGa                                                                      <  1:1853125/48‑1 (MQ=255)
                       cAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCATCa                                                >  1:4181784/1‑69 (MQ=255)
                        aaaCAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACGGAATAATCGCCAGCGTCGc                                                     >  1:758607/1‑63 (MQ=255)
                        aaaCAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCg                                                         <  1:3406202/59‑1 (MQ=255)
                        aaaCAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTc                                                       <  1:2653161/61‑1 (MQ=255)
                         aaCAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACGGaata                                                                   <  1:1890369/48‑1 (MQ=255)
                         aaCAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGc                                            <  1:557600/71‑1 (MQ=255)
                          aCAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACGGaat                                                                    <  1:3347993/46‑1 (MQ=255)
                          aCAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGc                                            >  1:2779617/1‑70 (MQ=255)
                          aCAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAAACGc                                                              >  1:400166/1‑52 (MQ=255)
                           cAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGaataat                                                                 >  1:4275383/1‑48 (MQ=255)
                            aGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACGGaataa                                                                  >  1:752989/1‑46 (MQ=255)
                            aGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACGGaataa                                                                  >  1:4014274/1‑46 (MQ=255)
                            aGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACGGaataa                                                                  >  1:162221/1‑46 (MQ=255)
                              tCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGa                                                                      <  1:1120360/40‑1 (MQ=255)
                               cGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCa                                                   >  1:2633545/1‑58 (MQ=255)
                                 cACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGc                                                 <  1:2957323/58‑1 (MQ=255)
                                 cACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAg                                         <  1:449434/66‑1 (MQ=255)
                                     aGGAAAAATGCGGCTAACGACTTACCAAACGGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCt                                 >  1:845477/1‑70 (MQ=255)
                                     aGGAAAAATGCGGCTAACGACTTACCAAACAGaataa                                                                  >  1:2493915/1‑37 (MQ=255)
                                     aGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCt                                 >  1:4277683/1‑70 (MQ=255)
                                      ggAAAAATGCGGCTAACGACTTACCAAACAGaataat                                                                 >  1:2514467/1‑37 (MQ=255)
                                      ggAAAAATGCGGCTAACGACTTACCAAACAGaataat                                                                 >  1:1583667/1‑37 (MQ=255)
                                       gAAAAATGCGGCTAACGACTTACCAAACGGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGaa                              <  1:3152401/71‑1 (MQ=255)
                                       gAAAAATGCGGCTAACGACTTACCAAACAGaataa                                                                  >  1:3243790/1‑35 (MQ=255)
                                         aaaaTGCGGCTAACGACTTACCAAACAGAATAATCGCCAga                                                          >  1:2204403/1‑40 (MQ=255)
                                          aaaTGCGGCTAACGACTTACCAAACAGAATAATAGCCAGCGTCGCCAGCACt                                              >  1:3001937/1‑52 (MQ=255)
                                               cGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAAt                       <  1:223229/70‑1 (MQ=255)
                                                       gACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCtta                               <  1:121847/54‑3 (MQ=255)
                                                          ttACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAAtt                      <  1:2645547/60‑1 (MQ=255)
                                                            aCCAAACGGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAaac          >  1:3237886/1‑70 (MQ=255)
                                                            aCCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCgg                                      >  1:1551589/1‑42 (MQ=255)
                                                            aCCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGt                                     >  1:2930234/1‑43 (MQ=255)
                                                            aCCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAtac          <  1:3746361/70‑1 (MQ=255)
                                                            aCCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAaac          <  1:2085906/70‑1 (MQ=255)
                                                             ccAAACGGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAaacaa        <  1:112891/71‑1 (MQ=255)
                                                             ccAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAaaca         >  1:2095597/1‑70 (MQ=255)
                                                              cAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAaacaa        >  1:4629791/1‑70 (MQ=255)
                                                                   aGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCa             <  1:3029850/60‑1 (MQ=255)
                                                                   aGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATccc  <  1:1850572/71‑1 (MQ=255)
                                                                   |                                                                      
TAATGCCGTACGTTGGGAGTGACCAAACAGTCGCACCAGGAAAAATGCGGCTAACGACTTACCAAACAGAATAATCGCCAGCGTCGCCAGCACTGCCAGCGGTTGCTGAATCAGAATTAACGGATCAAACAACATCCC  >  W3110S.gb/501458‑501595

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