Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,794,102 G→A 100% L398L (TTG→TTA gabP → gamma‑aminobutyrate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,794,1020GA96.0% 139.3 / ‑1.4 50L398L (TTG→TTAgabPgamma‑aminobutyrate transporter
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base A (23/25);  total (23/27)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.74e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGT  >  W3110S.gb/2794035‑2794167
                                                                   |                                                                 
taGTCATCGCCGTTTGACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACg                                                                 <  1:551243/70‑1 (MQ=255)
taGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACg                                                                 <  1:2822366/70‑1 (MQ=255)
taGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACg                                                                 <  1:2814545/70‑1 (MQ=255)
taGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACg                                                                 <  1:569687/70‑1 (MQ=255)
        gCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCt                                                        >  1:2238415/1‑71 (MQ=255)
                  aGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTgg                                              <  1:2116118/71‑1 (MQ=255)
                    ttGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCg                                                 <  1:3200145/66‑1 (MQ=255)
                    ttGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCg                                                 <  1:2549916/66‑1 (MQ=255)
                    ttGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCg                                                 <  1:1441074/66‑1 (MQ=255)
                    ttGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCg                                                 <  1:3635399/66‑1 (MQ=255)
                    ttGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                            <  1:1785846/71‑1 (MQ=255)
                    ttGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                            <  1:865950/71‑1 (MQ=255)
                    ttGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                            <  1:960699/71‑1 (MQ=255)
                    ttGCGGATACGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                            <  1:3315508/71‑1 (MQ=255)
                      gCGGATGAGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTc                                           <  1:3560204/70‑1 (MQ=255)
                       cGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTc                                           <  1:611663/69‑1 (MQ=255)
                          aTGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTg                                      <  1:3672140/71‑1 (MQ=255)
                          aTGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTg                                      <  1:1153590/71‑1 (MQ=255)
                             cgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCg                                                 >  1:1074272/1‑57 (MQ=255)
                             cgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCg                                                 >  1:3358391/1‑57 (MQ=255)
                             cgcAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCg                                                 >  1:221241/1‑57 (MQ=255)
                                aaaaTTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGt                                >  1:898674/1‑71 (MQ=255)
                                aaaaTTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGt                                >  1:988846/1‑71 (MQ=255)
                                aaaaTTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGt                                >  1:3035266/1‑71 (MQ=255)
                                aaaaTTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGt                                >  1:2920038/1‑71 (MQ=255)
                                aaaaTTCTGCGAGCAGAAGGAAGCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACATGGCTAg                                 >  1:2528331/1‑70 (MQ=255)
                                               gaaggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                            <  1:1876482/44‑1 (MQ=255)
                                                 aggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCt               >  1:304395/1‑71 (MQ=255)
                                                 aggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCt               >  1:3333774/1‑71 (MQ=255)
                                                  ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCtt              >  1:1061141/1‑71 (MQ=255)
                                                  ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCtt              <  1:3533436/71‑1 (MQ=255)
                                                  ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCtt              >  1:1909764/1‑71 (MQ=255)
                                                  ggaagCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCtt              >  1:36869/1‑71 (MQ=255)
                                                     agCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt           >  1:3433912/1‑71 (MQ=255)
                                                     agCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt           <  1:119116/71‑1 (MQ=255)
                                                     agCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt           >  1:764060/1‑71 (MQ=255)
                                                     agCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTtgt           >  1:1806580/1‑71 (MQ=255)
                                                      gCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                            <  1:3310673/37‑1 (MQ=255)
                                                      gCGAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                            <  1:1118891/37‑1 (MQ=255)
                                                        gAAATTCGCTTACGAATGTGGCTTTATCCGTGGCt                                            >  1:714137/1‑35 (MQ=255)
                                                        gAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCa                                          >  1:2511883/1‑37 (MQ=255)
                                                        gAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCa                                          >  1:1533951/1‑37 (MQ=255)
                                                        gAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCa                                          >  1:2149420/1‑37 (MQ=255)
                                                        gAAATTCGCTTACGAATGTGGCTTTATCCGTGGCTCa                                          >  1:2501445/1‑37 (MQ=255)
                                                          aaTTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg      <  1:1274533/71‑1 (MQ=255)
                                                          aaTTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg      <  1:70952/71‑1 (MQ=255)
                                                          aaTTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg      <  1:1194814/71‑1 (MQ=255)
                                                          aaTTCGCTTACGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTgg      <  1:3174066/71‑1 (MQ=255)
                                                              cGCTTACGAATGTGGCTTTATCCGTGGCTCGCCTGGCTAGTCATAGGCTTTATTACCTTTGTGTTGGTAg   >  1:455281/1‑70 (MQ=255)
                                                               gCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGt  <  1:3048099/70‑1 (MQ=255)
                                                                   |                                                                 
TAGTCATCGCCGTTTCACAGTTGCGGATGCGCAAAATTCTGCGAGCAGAAGGAAGCGAAATTCGCTTGCGAATGTGGCTTTATCCGTGGCTCACCTGGCTAGTCATCGGCTTTATTACCTTTGTGTTGGTAGT  >  W3110S.gb/2794035‑2794167

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