Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,095,396 G→T 16.1% S47S (TCC→TCA gntX ← gluconate periplasmic binding protein with phosphoribosyltransferase domain, GNT I system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,095,3960GT16.1% 83.4 / 3.7 31S47S (TCC→TCAgntXgluconate periplasmic binding protein with phosphoribosyltransferase domain, GNT I system
Reads supporting (aligned to +/- strand):  ref base G (6/20);  new base T (2/3);  total (8/23)
Fisher's exact test for biased strand distribution p-value = 5.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.33e-01

ATAGTCGGCAACCGTGACCAGTCTTTGCCAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCGGGTGGCGCGTGAACAGACCG  >  W3110S.gb/4095327‑4095466
                                                                     |                                                                      
aTAGTCGGCAACCGTGACCAGTCTTTGCCAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATggg                                                                       <  1:910968/71‑1 (MQ=255)
   gTCGGCAACCGTGACCAGTCTTTGCCAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAg                                                                     <  1:1546823/70‑1 (MQ=255)
   gTCGGCAACCGTGACCAGTCTTTGCCAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAg                                                                     <  1:1958304/70‑1 (MQ=255)
    tCGGCAACCGTGACCAGTCTTTGCCAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAgtg                                                                   <  1:783589/71‑1 (MQ=255)
       gCAACCTGACCAGTCTTTGCCAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGtgtgtg                                                               <  1:1252874/71‑1 (MQ=255)
              tGACCAGTCTTTGCCAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGt                                                         <  1:2511572/71‑1 (MQ=255)
                         tGCAAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGTAGTGTGTGGCGGGTAATCCACATtg                                              >  1:2916233/1‑71 (MQ=255)
                           ccAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAgtg                                                                   <  1:3668399/48‑1 (MQ=255)
                           ccAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAgtg                                                                   <  1:2421176/48‑1 (MQ=255)
                           ccAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAgtg                                                                   <  1:2319221/48‑1 (MQ=255)
                           ccAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACATtgt                                             <  1:694049/70‑1 (MQ=255)
                               ggcggcGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACAt                                                <  1:3105945/63‑1 (MQ=255)
                                gcggcgGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCAt                                       <  1:576566/71‑1 (MQ=255)
                                   gcggTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAAc                                    <  1:234120/71‑1 (MQ=255)
                                   gcggTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAAc                                    <  1:2040024/71‑1 (MQ=255)
                                     ggTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGt                                  <  1:189768/71‑1 (MQ=255)
                                      gTTTTTGCAGGCAGCGACCGCAGGGAAGATGTGAGTGTGTGGCGGGTAATCCACATtgg                                             >  1:3120379/1‑58 (MQ=255)
                                      gTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAgtg                                                                   <  1:3501705/37‑1 (MQ=255)
                                         tttGCAGGCAGCGACCGCAGGGAAGATGTGAGTGTGTGGCGGGTAATCCACATtggg                                            <  1:29703/57‑3 (MQ=255)
                                         tttGCAGGCAGCGACCGCAGGGAAGATGTGAGTGTGTGGCGGGTAATCCACATtggg                                            <  1:1461307/57‑3 (MQ=255)
                                         tttGCAGGCAGCGACAGCAGGGAAGATGTGAGTGTGTGGCGGGTAATCCACATtggg                                            <  1:372461/57‑3 (MQ=255)
                                                  agCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAAc                                    >  1:1642275/1‑56 (MQ=255)
                                                      aCCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTgg                                           >  1:3017985/1‑45 (MQ=255)
                                                             ggAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACg                                   <  1:3260201/46‑1 (MQ=255)
                                                              gAAGATGTGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCg                      >  1:3386054/1‑58 (MQ=255)
                                                              gAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACg                                   <  1:561478/45‑1 (MQ=255)
                                                              gAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCg                      >  1:1570907/1‑58 (MQ=255)
                                                              gAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCGGGTGGCGCGTGaa         <  1:2854228/71‑1 (MQ=255)
                                                               aaGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCGGGTGGCGCGTGaa         >  1:2939879/1‑70 (MQ=255)
                                                                aGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCGGGTGGCGCGTGAACa       <  1:933267/71‑1 (MQ=255)
                                                                 gATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCGTGTGGCGCGTGAACAg      <  1:503799/71‑1 (MQ=255)
                                                                   tGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCGGGTGGCGCGTGAACAGAc    >  1:3510436/1‑71 (MQ=255)
                                                                     ggAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCGGGTGGCGCGTGAACAGACCg  <  1:529972/71‑1 (MQ=255)
                                                                     |                                                                      
ATAGTCGGCAACCGTGACCAGTCTTTGCCAGGGCGGCGGTTTTTGCAGGCAGCGACCGCAGGGAAGATGGGAGTGTGTGGCGGGTAATCCACATTGTGGGCATAACGTTTTATCTGTGCGGGTGGCGCGTGAACAGACCG  >  W3110S.gb/4095327‑4095466

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: