Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,219,258 T→C 29.0% 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,2580TC29.0% 42.5 / 21.7 31F379S (TTC→TCC) ygjGputrescine:2‑oxoglutaric acid aminotransferase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base T (14/8);  new base C (3/6);  total (17/14)
Fisher's exact test for biased strand distribution p-value = 2.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGG  >  W3110S.gb/3219191‑3219312
                                                                   |                                                      
aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCg                                                     >  1:1104355/1‑71 (MQ=255)
aaTGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCg                                                     >  1:1714612/1‑71 (MQ=255)
              gCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCggg                                       <  1:1086526/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGaa                                     <  1:1161947/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGaa                                     <  1:1415637/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGaa                                     <  1:1506809/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGaa                                     <  1:735812/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGaa                                     <  1:1663716/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGaa                                     <  1:582596/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGaa                                     <  1:176169/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGaa                                     <  1:338094/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGaa                                     <  1:283831/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGaa                                     <  1:589474/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGaa                                     <  1:1801230/71‑1 (MQ=255)
                aGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGaa                                     <  1:694663/71‑1 (MQ=255)
                 gAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGaa                                     >  1:373785/1‑70 (MQ=255)
                   aCTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAatat                                  >  1:1610783/1‑71 (MQ=255)
                   aCTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAatat                                  >  1:1107028/1‑71 (MQ=255)
                   aCTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAatat                                  >  1:548440/1‑71 (MQ=255)
                   aCTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAatat                                  >  1:481924/1‑71 (MQ=255)
                   aCTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAata                                   >  1:310033/1‑70 (MQ=255)
                        ccGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCg                              <  1:791982/70‑1 (MQ=255)
                                           aaGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGt          >  1:748056/1‑71 (MQ=255)
                                           aaGGCGATATGTTGCTGGACGGTTCCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGt          >  1:125959/1‑71 (MQ=255)
                                                    tGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAgg  >  1:47832/1‑70 (MQ=255)
                                                    tGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAgg  >  1:1942193/1‑70 (MQ=255)
                                                    tGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAgg  >  1:1824034/1‑70 (MQ=255)
                                                    tGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAgg  >  1:595332/1‑70 (MQ=255)
                                                    tGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAgg  >  1:626924/1‑70 (MQ=255)
                                                    tGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAgg  >  1:1787257/1‑70 (MQ=255)
                                                    tGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAgg  >  1:1521233/1‑70 (MQ=255)
                                                                   |                                                      
AATGTGTTGCTGGAGCAGAACTTACCGGCTCAGGCTGAGCAAAAAGGCGATATGTTGCTGGACGGTTTCCGTCAACTGGCGCGGGAATATCCCGATCTGGTACAGGAAGCGCGTGGTAAAGG  >  W3110S.gb/3219191‑3219312

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: