Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,999,315 C→A 8.9% E56* (GAA→TAA)  dcyD ← D‑cysteine desulfhydrase, PLP‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,999,3150CA 8.9% 154.0 / 3.4 45E56* (GAA→TAA) dcyDD‑cysteine desulfhydrase, PLP‑dependent
Reads supporting (aligned to +/- strand):  ref base C (27/14);  new base A (2/2);  total (29/16)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

ATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCG  >  NC_000913/1999178‑1999436
                                                                                                                                         |                                                                                                                          
aTTTTCCAGCAGCGCCACGCAGTGCTGACCTAGGTTGGCAGCGACTGCGGCAGTCTGCCGCACATCGTTAGAATGGATCGCCCCGGCCCTAATCCGCGTAACGGCACCTTCACTCAGATCATCAGCCCCTAGAACATCCa                                                                                                                          >  2:332348/1‑140 (MQ=255)
    tCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCtt                                                                                                                      >  2:335450/1‑140 (MQ=255)
    tCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCtt                                                                                                                      >  1:955579/1‑140 (MQ=255)
       agcagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACg                                                                                                                   <  2:202261/140‑1 (MQ=255)
         cagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTa                                                                                                                 <  1:733536/140‑1 (MQ=255)
         cagcGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTa                                                                                                                 >  2:490096/1‑140 (MQ=255)
                        cAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCa                                                                                                 >  2:290093/1‑140 (MQ=255)
                        cAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCa                                                                                                 <  2:1033867/140‑1 (MQ=255)
                         aGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCAt                                                                                                <  1:899697/140‑1 (MQ=255)
                         aGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCAt                                                                                                <  1:518173/140‑1 (MQ=255)
                           acagAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTACAGCTTACGTAATTTATT‑GCCGCCCATTg                                                                                              >  1:982124/4‑140 (MQ=255)
                           aCCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTg                                                                                              >  1:658342/1‑140 (MQ=255)
                           aCCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTACAGCTTACGTAATTTATT‑GCCGCCCATTg                                                                                              >  1:982789/1‑140 (MQ=255)
                             cGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGcc                                                                                            >  1:618108/1‑140 (MQ=255)
                                gTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATg                                                                                         >  1:1052744/1‑140 (MQ=255)
                                gTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATg                                                                                         >  1:808922/1‑140 (MQ=255)
                                gTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATg                                                                                         >  1:40995/1‑140 (MQ=255)
                                gTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATg                                                                                         >  1:1152526/1‑140 (MQ=255)
                                     gCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGt                                                                                    <  1:898670/140‑1 (MQ=255)
                                                     tCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTtg                                                                    >  2:400167/1‑140 (MQ=255)
                                                               aTGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAAttt                                                          >  2:423589/1‑140 (MQ=255)
                                                                tGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTACAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTc                                                         <  2:982124/140‑1 (MQ=255)
                                                                tGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTACAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTc                                                         <  2:982789/140‑1 (MQ=255)
                                                                tGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCGTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCAGTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTc                                                         <  2:955744/140‑1 (MQ=255)
                                                                           gATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGGCCTTTGCCGTGGGGGTGTCGTCTTCGCGTTTGATGAATGTTTCAGGTCCTTGAt                                              >  2:1220386/1‑140 (MQ=255)
                                                                           gATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGAt                                              >  2:284246/1‑140 (MQ=255)
                                                                           gATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGAt                                              >  1:347168/1‑140 (MQ=255)
                                                                           gATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGAt                                              >  1:346906/1‑140 (MQ=255)
                                                                           gATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGAt                                              >  1:269380/1‑140 (MQ=255)
                                                                           gATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGAt                                              >  2:339237/1‑140 (MQ=255)
                                                                           gATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGAt                                              >  1:295872/1‑140 (MQ=255)
                                                                           gATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGAt                                              >  1:295771/1‑140 (MQ=255)
                                                                                 cccGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAg                                        >  1:356689/1‑140 (MQ=255)
                                                                                 cccGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAg                                        >  1:356605/1‑140 (MQ=255)
                                                                                        gTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAgcg                                 <  1:212371/140‑1 (MQ=255)
                                                                                           aTCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCgg                              >  1:767837/1‑140 (MQ=255)
                                                                                               gCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGa                          <  1:499898/140‑1 (MQ=255)
                                                                                                       gCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGGTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAg                  <  2:880691/140‑1 (MQ=255)
                                                                                                       gCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAg                  <  2:457327/140‑1 (MQ=255)
                                                                                                       gCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAg                  <  2:457385/140‑1 (MQ=255)
                                                                                                       gCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAg                  <  1:345380/140‑1 (MQ=255)
                                                                                                       gCACCGTCCCGCAGAGCAACTGCCGCGAGAACATCCAGCTTACG‑AACTTATTGGCCGCCCATAGCCATGGGGGTGACGTCAACCCGTTAGAAGACCAGTCCCCGTCCAAGATAATCCGAAAAGCGCGGCAGATATTCgcg                  <  2:937681/140‑3 (MQ=255)
                                                                                                                     agCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGgcgcgc    >  1:1125197/1‑140 (MQ=255)
                                                                                                                       cATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCg  >  1:482093/1‑140 (MQ=255)
                                                                                                                       cATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCg  >  1:482369/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                          
ATTTTCCAGCAGCGCCACGCAGTGCAGACCGAGTTTGGCAGCGACTGCGGCAGTCTGCCGCACATGGTTAGACTGGATCGCCCCGGCAGTAATCAGCGTATCGGCACCTTCACGCAGAGCATCTGCCGCGAGAAATTCCAGCTTACGTAATTTATT‑GCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCG  >  NC_000913/1999178‑1999436

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑