Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,000,742 C→T 58.1% M292I (ATG→ATA yedO ← D‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,000,7420CT58.1% 11.0 / 36.6 31M292I (ATG→ATAyedOD‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (7/6);  new base T (14/4);  total (21/10)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.74e-01

TGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTC  >  W3110S.gb/2000677‑2000810
                                                                 |                                                                    
tGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCg                                                                  <  1:2260848/70‑1 (MQ=255)
  aaTAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCtt                                                               >  1:1091053/1‑71 (MQ=255)
  aaTAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtt                                                               >  1:1756071/1‑71 (MQ=255)
   aTAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCttt                                                              >  1:1357169/1‑71 (MQ=255)
      aaCAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCttt                                                              >  1:415709/1‑68 (MQ=255)
      aaCAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTcc                                                           >  1:250965/1‑71 (MQ=255)
        cAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCgg                                                         <  1:1197007/71‑1 (MQ=255)
             tCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCGGtata                                                     >  1:969521/1‑70 (MQ=255)
               ggCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCGGTATaca                                                   <  1:2251498/70‑1 (MQ=255)
                   cTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCgg                                                         <  1:740114/60‑1 (MQ=255)
                      cATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATc                                            >  1:2439539/1‑70 (MQ=255)
                          tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAg                                         >  1:1014893/1‑69 (MQ=255)
                          tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAg                                         >  1:1375825/1‑69 (MQ=255)
                          tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATc                                                                   <  1:1946737/43‑1 (MQ=255)
                          tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATc                                                                   <  1:372575/43‑1 (MQ=255)
                          tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATc                                                                   <  1:193101/43‑1 (MQ=255)
                          tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATc                                                                   <  1:119692/43‑1 (MQ=255)
                          tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGta                                                       >  1:2138104/1‑55 (MQ=255)
                               aGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAGCAGAATg                                  >  1:1555249/1‑71 (MQ=255)
                                     tCTGACTGATACCGTCAATCAGCCCCGCTATCGCtt                                                               >  1:1100602/1‑36 (MQ=255)
                                          ctgaTACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAg                                         >  1:568478/1‑53 (MQ=255)
                                            gaTACCGTCAATCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCcgcg                     >  1:2666159/1‑71 (MQ=255)
                                                  gTCAATCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGc                >  1:568653/1‑70 (MQ=255)
                                                     aaTCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAAtt            >  1:882928/1‑71 (MQ=255)
                                                     aaTCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAAtt            >  1:1176693/1‑71 (MQ=255)
                                                      aTCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGc                                 <  1:1078919/49‑1 (MQ=255)
                                                         aGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCAc        >  1:1069674/1‑71 (MQ=255)
                                                            cccGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAg                          >  1:701880/1‑50 (MQ=255)
                                                            cccGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAg                          >  1:2681392/1‑50 (MQ=255)
                                                            cccGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGc     >  1:2367500/1‑71 (MQ=255)
                                                               gCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTc  <  1:2142410/71‑1 (MQ=255)
                                                                 |                                                                    
TGAATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTC  >  W3110S.gb/2000677‑2000810

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: