Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,002,965 A→G 8.3% L493S (TTA→TCA)  ygeV ← predicted DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,002,9650AG 8.3% 201.2 / 7.9 72L493S (TTA→TCA) ygeVpredicted DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/31);  new base G (4/2);  total (39/33)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCAATGCCAGGAT  >  W3110S.gb/3002902‑3003034
                                                               |                                                                     
aGCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTcgcg                                                                >  1:737552/1‑71 (MQ=255)
aGCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTcgcg                                                                >  1:1236077/1‑71 (MQ=255)
aGCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTcgcg                                                                >  1:1360445/1‑71 (MQ=255)
aGCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTcgcg                                                                >  1:527323/1‑71 (MQ=255)
 gCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCATATATTCCATCAAATTGCTTGACTcgcg                                                                >  1:1031932/1‑70 (MQ=255)
 gCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTGACTcgcg                                                                >  1:598483/1‑70 (MQ=255)
 gCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTGACTcgcg                                                                >  1:2971892/1‑70 (MQ=255)
 gCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTGACTcgcg                                                                >  1:2063688/1‑70 (MQ=255)
         tCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTcc                                                       >  1:2202958/1‑71 (MQ=255)
           gATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTcgc                                                                 >  1:1743575/1‑59 (MQ=255)
           gATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCgg                                                     >  1:2451274/1‑71 (MQ=255)
           gATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCgg                                                     >  1:3531746/1‑71 (MQ=255)
               aCTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGcc                                                  <  1:502777/70‑1 (MQ=255)
               aCTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGcc                                                  <  1:3568640/70‑1 (MQ=255)
                 ttCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGttt                                                         <  1:499882/61‑1 (MQ=255)
                 ttCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGttt                                                         <  1:2620885/61‑1 (MQ=255)
                 ttCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGttt                                                         <  1:2060168/61‑1 (MQ=255)
                 ttCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGttt                                                         <  1:1691564/61‑1 (MQ=255)
                 ttCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGttt                                                         <  1:619773/61‑1 (MQ=255)
                  tCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAAc                                               >  1:1301223/1‑70 (MQ=255)
                  tCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAAc                                               >  1:529908/1‑70 (MQ=255)
                   cACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAAc                                               <  1:702915/69‑1 (MQ=255)
                   cACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGa                                             >  1:3533026/1‑71 (MQ=255)
                   cACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGa                                             >  1:2324337/1‑71 (MQ=255)
                       tGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAAc                                               <  1:68415/65‑1 (MQ=255)
                         aaGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTa                                       <  1:775162/71‑1 (MQ=255)
                         aaGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTa                                       <  1:1259215/71‑1 (MQ=255)
                          aGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTa                                       <  1:2734836/70‑1 (MQ=255)
                          aGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGAAg                                                             >  1:792980/1‑48 (MQ=255)
                             aaCCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCa                                   <  1:1608875/71‑1 (MQ=255)
                              aCCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCa                                   <  1:1313996/70‑1 (MQ=255)
                                cACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAata                                >  1:2789712/1‑71 (MQ=255)
                                cACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAata                                >  1:268020/1‑71 (MQ=255)
                                 aCGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGa                                             >  1:672738/1‑57 (MQ=255)
                                 aCGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATg                                           <  1:457337/59‑1 (MQ=255)
                                 aCGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATg                                           <  1:3247418/59‑1 (MQ=255)
                                 aCGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATg                                           <  1:2109002/59‑1 (MQ=255)
                                    ttAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTc                            >  1:1308977/1‑71 (MQ=255)
                                    ttAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTc                            >  1:1630281/1‑71 (MQ=255)
                                       aCCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTcgac                         <  1:3134844/71‑1 (MQ=255)
                                       aCCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTcgac                         <  1:573961/71‑1 (MQ=255)
                                          aGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCgacgac                      <  1:149482/71‑1 (MQ=255)
                                          aGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCgacgac                      <  1:1563488/71‑1 (MQ=255)
                                           gATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGTAATATTTCGACGACa                     <  1:565912/71‑1 (MQ=255)
                                           gATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACa                     <  1:1922896/71‑1 (MQ=255)
                                           gATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACa                     <  1:2799418/71‑1 (MQ=255)
                                             tatTCCATCAAATTGCTTGACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACAt                    <  1:106943/70‑1 (MQ=255)
                                             tatTCCATCAAATTGCTTGACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACAt                    <  1:266729/70‑1 (MQ=255)
                                             tatTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATc                   <  1:1219676/71‑1 (MQ=255)
                                             tatTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATc                   <  1:2589583/71‑1 (MQ=255)
                                             tatTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATc                   <  1:1331804/71‑1 (MQ=255)
                                              atTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCa                                   <  1:306937/54‑1 (MQ=255)
                                              atTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCa                                   <  1:2909136/54‑1 (MQ=255)
                                              atTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCg                  >  1:321150/1‑71 (MQ=255)
                                              atTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCg                  >  1:1381822/1‑71 (MQ=255)
                                                 ccATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGa                                             >  1:464936/1‑41 (MQ=255)
                                                 ccATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGa                                             >  1:153737/1‑41 (MQ=255)
                                                 ccATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCgg                 >  1:1165785/1‑69 (MQ=255)
                                                 ccATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCgg                 >  1:1308204/1‑69 (MQ=255)
                                                     cAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCaa           >  1:3077342/1‑71 (MQ=255)
                                                      aaaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAata                                >  1:3473489/1‑49 (MQ=255)
                                                      aaaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAata                                >  1:895717/1‑49 (MQ=255)
                                                      aaaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAata                                >  1:441822/1‑49 (MQ=255)
                                                      aaaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTcga                          >  1:1626820/1‑55 (MQ=255)
                                                      aaaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTcga                          >  1:3663947/1‑55 (MQ=255)
                                                      aaaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTcga                          >  1:112957/1‑55 (MQ=255)
                                                      aaaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCAAt          >  1:586235/1‑71 (MQ=255)
                                                      aaaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCAAt          >  1:594340/1‑71 (MQ=255)
                                                      aaaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCAAt          <  1:35951/71‑1 (MQ=255)
                                                       aaTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCaa                                  <  1:1645306/46‑1 (MQ=255)
                                                        aTTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACa                     >  1:553718/1‑58 (MQ=255)
                                                          tGCTTAACTCGCGCAGGTTTCAGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCAATGc        >  1:1675368/1‑69 (MQ=255)
                                                               aaCTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCAATGCCAGGAt  >  1:1401253/1‑70 (MQ=255)
                                                               aaCTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCAATGCCAGGAt  >  1:222711/1‑70 (MQ=255)
                                                               aaCTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCAATGCCAGGAt  >  1:2504312/1‑70 (MQ=255)
                                                               |                                                                     
AGCGTGCTGTCGATAACTTCACCTGAAGGAACCACGTTAACCAGATATTCCATCAAATTGCTTAACTCGCGCAGGTTTCCGGGCCAACGATGCTTACGCAATATTTCGACGACATCGGGAGCAATGCCAGGAT  >  W3110S.gb/3002902‑3003034

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