Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 633,015 G→T 51.5% V69V (GTG→GTT ybdL → methionine aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb633,0150GT51.5% 13.5 / 58.6 33V69V (GTG→GTTybdLmethionine aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base G (6/10);  new base T (12/5);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 8.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.40e-01

TACAGGAGCGGCTGGCGCACCACGTTGCACAGGGGGCAAACCAATACGCGCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCAGATGCCGATAGC  >  W3110S.gb/632951‑633084
                                                                |                                                                     
tACAGGAGCGGCTGGCGCACCACGTTGCACAGGGGGCAAACCAATACGCGCCCATGACCGGCGTTCAGGcc                                                                 >  1:2037178/1‑71 (MQ=255)
tACAGGAGCGGCTGGCGCACCACGTTGCACAGGGGGCAAACCAATACGCGCCCATGACCGGCGTTCAGGcc                                                                 >  1:2042187/1‑71 (MQ=255)
tACAGGAGCGGCTGGCGCACCACGTTGCACAGGGGGCAAACCAATACGCGCCCATGACCGGCGTGCAGGcc                                                                 >  1:2160033/1‑71 (MQ=255)
                   ccacGTTGCACAGGGGGCAAACCAATACGCGCCCATGACCGGCGTGCAGGc                                                                  <  1:3214901/51‑1 (MQ=255)
                   ccacGTTGCACAGGGGGCAAACCAATACGCGCCCATGACCGGCGTGCAGGc                                                                  <  1:957233/51‑1 (MQ=255)
                   ccacGTTGCACAGGGGGCAAACCAATACGCGCCCATGACCGGCGTGCAGGc                                                                  <  1:3542252/51‑1 (MQ=255)
                              aGGGGGCAAACCAATACGCGCCCCTGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGaa                                   >  1:2300613/1‑71 (MQ=255)
                              aGGGGGCAAACCAATACGCGCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGaa                                   >  1:4060674/1‑71 (MQ=255)
                              aGGGGGCAAACCAATACGCGCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGaa                                   >  1:225510/1‑71 (MQ=255)
                              aGGGGGCAAACCAATACGCGCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGaa                                   >  1:3651619/1‑71 (MQ=255)
                               gggggCAAACCAATACGCGCCCATGACCGGCGTTCAGGCCTTGCGCGAGGCGa                                                    >  1:3476830/1‑53 (MQ=255)
                               gggggCAAACCAATACGCGCACATGACCGGCGTTCAGGCCTTGCGCGAGGCGa                                                    >  1:2815670/1‑53 (MQ=255)
                                    cAAACCAATACGCGCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTg                             <  1:1464019/71‑1 (MQ=255)
                                     aaaCCAATACGCGCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTg                             <  1:981428/70‑1 (MQ=255)
                                     aaaCCAATACGCGCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTg                             <  1:2430510/70‑1 (MQ=255)
                                        ccAATACGCGCCCATGACCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATg                         >  1:2694269/1‑71 (MQ=255)
                                        ccAATACGCGCCCATGACCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATg                         >  1:3837509/1‑71 (MQ=255)
                                        ccAATACGCGCCCATGACCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATg                         >  1:898035/1‑71 (MQ=255)
                                        ccAATACGCGCCCATGACCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATg                         >  1:2962754/1‑71 (MQ=255)
                                            tACGCGCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTg                             <  1:2726242/63‑1 (MQ=255)
                                              cgcgCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTAt                    >  1:1292163/1‑70 (MQ=255)
                                                   ccATGACCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCa              <  1:1200115/71‑1 (MQ=255)
                                                   ccATGACCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCa              <  1:2511970/71‑1 (MQ=255)
                                                   ccATGACCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCa              <  1:2359922/71‑1 (MQ=255)
                                                      tGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCa              <  1:4037060/68‑1 (MQ=255)
                                                        aCCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTg                             >  1:3777105/1‑51 (MQ=255)
                                                        aCCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTa                           >  1:1995319/1‑53 (MQ=255)
                                                        aCCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCAGATGc         >  1:1806560/1‑71 (MQ=255)
                                                        aCCGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTCTATGGCTATCAACCAGATGc         >  1:1393963/1‑71 (MQ=255)
                                                         ccGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCGGAAAACGGAACGTTTGTATGGCTATCAACCAGATGcc        <  1:4068070/71‑1 (MQ=255)
                                                         ccGGCGTTCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCAGATGcc        <  1:2894150/71‑1 (MQ=255)
                                                               tGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCAGATGCCGATAGc  <  1:470326/71‑1 (MQ=255)
                                                               tGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCAGATGCCGATAGc  <  1:1719774/71‑1 (MQ=255)
                                                                |                                                                     
TACAGGAGCGGCTGGCGCACCACGTTGCACAGGGGGCAAACCAATACGCGCCCATGACCGGCGTGCAGGCCTTGCGCGAGGCGATTGCTCAGAAAACGGAACGTTTGTATGGCTATCAACCAGATGCCGATAGC  >  W3110S.gb/632951‑633084

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: