Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,000,742 C→T 13.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,7420CT13.1% 103.8 / 7.1 38M292I (ATG→ATAyedOD‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (17/16);  new base T (2/3);  total (19/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

AATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTC  >  W3110S.gb/2000679‑2000810
                                                               |                                                                    
aaTAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTTATACCGTCAATCAGCCCCGCTATCGc                                                                 <  1:3211655/69‑1 (MQ=255)
aaTAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGc                                                                 >  1:2470100/1‑69 (MQ=255)
aaTAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGc                                                                 >  1:1861645/1‑69 (MQ=255)
       aGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCtttt                                                             >  1:2852552/1‑66 (MQ=255)
         aaTCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt                                                        <  1:600554/69‑1 (MQ=255)
         aaTCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt                                                        <  1:437617/69‑1 (MQ=255)
         aaTCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt                                                        <  1:300887/69‑1 (MQ=255)
          aTCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGta                                                       <  1:2279421/69‑1 (MQ=255)
          aTCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGt                                                        <  1:296592/68‑1 (MQ=255)
            cGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGtata                                                     <  1:3110024/69‑1 (MQ=255)
            cGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGtata                                                     >  1:2767967/1‑69 (MQ=255)
                    cATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGa                                              >  1:2315981/1‑68 (MQ=255)
                     aTCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGcc                                                                      <  1:3032896/43‑1 (MQ=255)
                     aTCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGcc                                                                      <  1:1596686/43‑1 (MQ=255)
                     aTCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGcc                                                                      <  1:1128681/43‑1 (MQ=255)
                      tCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGcc                                                                      <  1:1252645/42‑1 (MQ=255)
                        tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCgg                                                         >  1:260310/1‑53 (MQ=255)
                        tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCGGTATaca                                                   <  1:1580651/59‑1 (MQ=255)
                        tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCa                                           >  1:3434878/1‑67 (MQ=255)
                        tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCa                                           >  1:3440333/1‑67 (MQ=255)
                        tttGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCa                                           >  1:402427/1‑67 (MQ=255)
                            aaGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAg                                      <  1:1139182/68‑1 (MQ=255)
                            aaGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAg                                      <  1:1510414/68‑1 (MQ=255)
                                   tCTGACTGATACCGTCAATCAGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCt                               >  1:2382947/1‑68 (MQ=255)
                                   tCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATacac                                                  <  1:3523175/49‑1 (MQ=255)
                                   tCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCa                                       <  1:2908034/60‑1 (MQ=255)
                                   tCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAg                                      <  1:1745024/61‑1 (MQ=255)
                                   tCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCt                               >  1:1524361/1‑68 (MQ=255)
                                   tCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATAAACAGGATCAAGCAg                                      <  1:3091697/61‑1 (MQ=255)
                                              ccGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCa                                       <  1:1793047/49‑1 (MQ=255)
                                              ccGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGcc                   >  1:2225875/1‑69 (MQ=255)
                                              ccGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGcc                   >  1:223445/1‑69 (MQ=255)
                                                   aaTCAGCCCCGCCATCGCTTTTCCGGTATTCACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCaa              >  1:122907/1‑69 (MQ=255)
                                                   aaTCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAg                                      <  1:1976641/45‑1 (MQ=255)
                                                   aaTCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAg                                      >  1:2457359/1‑45 (MQ=255)
                                                   aaTCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAg                                      <  1:2726101/45‑1 (MQ=255)
                                                   aaTCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCaa              >  1:3021567/1‑69 (MQ=255)
                                                   aaTCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCaa              >  1:2889490/1‑69 (MQ=255)
                                                      cAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAAttt           <  1:1325279/69‑1 (MQ=255)
                                                       aGCCCCGCTATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTc          <  1:3529932/69‑1 (MQ=255)
                                                       aGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCa               >  1:1381249/1‑64 (MQ=255)
                                                       aGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAAttt           >  1:1717552/1‑68 (MQ=255)
                                                               cATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTc  <  1:1877358/69‑1 (MQ=255)
                                                               cATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTc  <  1:2215065/69‑1 (MQ=255)
                                                               |                                                                    
AATAAACAGAATCGGCCCTTCATCTTTGAAGCGTTTCTGACTGATACCGTCAATCAGCCCCGCCATCGCTTTTCCGGTATACACAGGATCAAGCAGAATGCCTTCCAGCCGCGCCAGCAATTTCACTGCTTC  >  W3110S.gb/2000679‑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: