Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 1,588,010 T→G 48.0% D64A (GAC→GCC)  ydeP ← predicted oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb1,588,0100TG48.0% 1.4 / 40.1 25D64A (GAC→GCC) ydePpredicted oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (9/4);  new base G (5/7);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 2.38e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAT  >  W3110S.gb/1587941‑1588080
                                                                     |                                                                      
aagaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTc                                                                       >  1:1127211/1‑71 (MQ=255)
aagaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTc                                                                       >  1:570659/1‑71 (MQ=255)
aagaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTc                                                                       >  1:196718/1‑71 (MQ=255)
aagaaagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTc                                                                       >  1:1950529/1‑71 (MQ=255)
    aagaGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTCGCCGTTTTCACAAATGGCGAAt                                                                   >  1:1360857/1‑70 (MQ=255)
              tACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACtgt                                                         >  1:391960/1‑71 (MQ=255)
              tACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACtgt                                                         >  1:1761527/1‑71 (MQ=255)
                       aTCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGAt                                                >  1:279247/1‑71 (MQ=255)
                       aTCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGAt                                                >  1:411950/1‑71 (MQ=255)
                       aTCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGAt                                                >  1:1832991/1‑71 (MQ=255)
                            tGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGc                                           <  1:1081013/71‑1 (MQ=255)
                            tGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGc                                           <  1:336646/71‑1 (MQ=255)
                            tGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTACTTAGGATCTGGc                                           <  1:500803/71‑1 (MQ=255)
                            tGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGc                                           <  1:1658998/71‑1 (MQ=255)
                            tGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGc                                           <  1:1461063/71‑1 (MQ=255)
                            tGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGc                                           <  1:581130/71‑1 (MQ=255)
                            tGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGc                                           <  1:1180943/71‑1 (MQ=255)
                            tGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTAACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGc                                           <  1:1715888/71‑1 (MQ=255)
                              aCTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGc                                           <  1:1721988/69‑1 (MQ=255)
                                  cccAGGCGATTGCTTTTGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCAATg                                       <  1:726054/69‑1 (MQ=255)
                                                  tGCGCCGTTTTCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCa                      <  1:60250/70‑1 (MQ=255)
                                                           ttCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACaa                                  >  1:760085/1‑49 (MQ=255)
                                                           ttCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACaa                                  >  1:1720634/1‑49 (MQ=255)
                                                           ttCACAAATGGCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACaa                                  >  1:1455899/1‑49 (MQ=255)
                                                           ttCACAAATGGCGAATGAAGCACTGTGCTTAGGATCTGGCCATGCACaa                                  >  1:420981/1‑49 (MQ=39)
                                                                     tCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAt  >  1:753853/1‑71 (MQ=255)
                                                                     tCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAt  >  1:828080/1‑71 (MQ=255)
                                                                     tCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAt  >  1:69253/1‑71 (MQ=255)
                                                                     tCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAt  >  1:601254/1‑71 (MQ=255)
                                                                     tCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAt  >  1:500301/1‑71 (MQ=255)
                                                                     tCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAt  >  1:1370601/1‑71 (MQ=255)
                                                                     tCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCAGACAGTCAAAGCCCTCTGGCTTATTCAt  >  1:1216205/1‑71 (MQ=255)
                                                                     |                                                                      
AAGAAAGAGGCGTTTACCTGCTTATCCGTGACTTCCCAGGCGATTGCTTTTGCGCCGTTTTCACAAATGTCGAATGACGCACTGTGCTTAGGATCTGGCCATGCACAACCCGGACAGTCAAAGCCCTCTGGCTTATTCAT  >  W3110S.gb/1587941‑1588080

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: