Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,219,258 T→C 18.9% F379S (TTC→TCC)  ygjG → putrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,219,2580TC18.9% 115.8 / 22.5 53F379S (TTC→TCC) ygjGputrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (24/19);  new base C (2/8);  total (26/27)
Fisher's exact test for biased strand distribution p-value = 7.64e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.24e-01

AATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTT  >  W3110S.gb/3219191‑3219318
                                                                   |                                                            
aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGttt                                                              >  1:2849511/1‑68 (MQ=255)
aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTc                                                             >  1:92978/1‑69 (MQ=255)
aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTc                                                             >  1:889797/1‑69 (MQ=255)
  tgtgTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTc                                                             >  1:1971379/1‑67 (MQ=255)
  tgtgTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTc                                                             >  1:3514414/1‑67 (MQ=255)
    tgtTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGt                                                          >  1:1871086/1‑68 (MQ=255)
    tgtTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTc                                                         <  1:894198/69‑1 (MQ=255)
    tgtTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTc                                                         >  1:1423665/1‑69 (MQ=255)
    tgtTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTc                                                         <  1:2857468/69‑1 (MQ=255)
      ttGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCa                                                        <  1:211573/68‑1 (MQ=255)
                  aaCTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGaa                                           <  1:306570/69‑1 (MQ=255)
                  aaCTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCggg                                             <  1:2574450/67‑1 (MQ=255)
                   aCTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAat                                          >  1:1273251/1‑69 (MQ=255)
                    cTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGa                                            >  1:2447013/1‑66 (MQ=255)
                    cTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGAata                                         <  1:2922288/69‑1 (MQ=255)
                     ttACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAatat                                        <  1:438150/69‑1 (MQ=255)
                     ttACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAatat                                        <  1:2369992/69‑1 (MQ=255)
                     ttACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGAata                                         <  1:1885707/68‑1 (MQ=255)
                      tACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATc                                       >  1:2364897/1‑69 (MQ=255)
                      tACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATc                                       >  1:3097024/1‑69 (MQ=255)
                      tACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATc                                       >  1:894536/1‑69 (MQ=255)
                      tACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATc                                       >  1:1165009/1‑69 (MQ=255)
                       aCCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATc                                       >  1:188649/1‑68 (MQ=255)
                          ggCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGa                                   >  1:986613/1‑69 (MQ=255)
                          ggCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGa                                   >  1:2459702/1‑69 (MQ=255)
                               aGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATccc                                     >  1:3243285/1‑62 (MQ=255)
                                 gCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGAATATCCCGATCTgg                              <  1:2863422/67‑1 (MQ=255)
                                   tGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTAc                           >  1:1394615/1‑68 (MQ=255)
                                   tGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATATGGTAc                           >  1:1401874/1‑68 (MQ=255)
                                        aaaaaGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAg                     >  1:1400171/1‑69 (MQ=255)
                                            aGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGcgcg                 <  1:689120/69‑1 (MQ=255)
                                            aGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGcgcg                 >  1:2349823/1‑69 (MQ=255)
                                            aGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGcgcg                 >  1:167306/1‑69 (MQ=255)
                                                 atatGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTa            <  1:1786750/69‑1 (MQ=255)
                                                 atatGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTa            <  1:1254240/69‑1 (MQ=255)
                                                      ttGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAggg       >  1:779735/1‑69 (MQ=255)
                                                      ttGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAggg       >  1:437106/1‑69 (MQ=255)
                                                      ttGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAggg       >  1:1459855/1‑69 (MQ=255)
                                                      ttGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAggg       >  1:3103282/1‑69 (MQ=255)
                                                      ttGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAggg       >  1:3443072/1‑69 (MQ=255)
                                                        gCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGa      <  1:3423210/68‑1 (MQ=255)
                                                           ggACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:56458/69‑1 (MQ=255)
                                                           ggACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:1458015/69‑1 (MQ=255)
                                                           ggACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:1871242/69‑1 (MQ=255)
                                                           ggACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:152877/69‑1 (MQ=255)
                                                           ggACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:1688146/69‑1 (MQ=255)
                                                           ggACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:2828735/69‑1 (MQ=255)
                                                           ggACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:2782833/69‑1 (MQ=255)
                                                           ggACGGTTCCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:1014457/69‑1 (MQ=255)
                                                           ggACGGTTCCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:1018632/69‑1 (MQ=255)
                                                              cGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:3082800/66‑1 (MQ=255)
                                                              cGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:2492452/66‑1 (MQ=255)
                                                              cGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:1391366/66‑1 (MQ=255)
                                                               ggTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGtt  <  1:1489365/65‑1 (MQ=255)
                                                                  ttCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTa            <  1:2924698/52‑1 (MQ=255)
                                                                   |                                                            
AATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGGGATGTT  >  W3110S.gb/3219191‑3219318

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: