Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,999,460	G→T	16.5%	T7T (ACC→ACA)	dcyD ←	D‑cysteine desulfhydrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,999,460	0	G	T	16.5%	87.6 / 6.2	41	T7T (ACC→ACA)	dcyD	D‑cysteine desulfhydrase
Reads supporting (aligned to +/- strand):  ref base G (19/14);  new base T (5/2);  total (24/17)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATC  >  NC_000913/1999319‑1999599                                                                                                                                              |                                                                                                                                            cTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACggg                                                                                                                                            <  1:686870/143‑1 (MQ=255)      gTATTTTATTGCCGCCCATTGCCCTGGTGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTaaa                                                                                                                                       <  2:579710/143‑1 (MQ=255)           ttAGTGCCGCCCATTGCCAGGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGTGTTAAATTATg                                                                                                                                  <  1:261265/143‑1 (MQ=255)                    cccATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGACAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACCCCCGCCGTGGAAAACGGGCTAACTTCTGCAGCGGCAt                                                                                                                         >  2:322306/1‑143 (MQ=255)                     ccATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATa                                                                                                                        >  1:835613/1‑143 (MQ=255)                            cATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGTGGTAGGTTATGCAGTGGCATAAGAGCCt                                                                                                                 >  1:23150/1‑143 (MQ=255)                               gggggTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCg                                                                                                              <  2:462239/143‑1 (MQ=255)                               gggcgTGACGTCATCCCGTTGGATGAAAATTTCCCGTCCTGGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCg                                                                                                              <  2:563654/143‑1 (MQ=255)                                ggggTGACGCCATCCCGTTTGATGAAACTTTCCCGTCCCAGATAATCAGAAAAGCGCGGCAGATATTCGAGCTGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGa                                                                                                             <  1:46578/143‑1 (MQ=255)                                    tGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGAtgtg                                                                                                         >  1:731514/1‑143 (MQ=255)                                     gACGTCATCCCGTTTGATGAAAATTCCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGTTTAAATTATGCAGTGGCATAACAGCCTCCGAtgtgt                                                                                                        <  1:612160/143‑1 (MQ=255)                                            tCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAAGTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGtgttgttg                                                                                                 >  1:452514/1‑143 (MQ=255)                                            tCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGtgctgttg                                                                                                 >  1:478528/1‑143 (MQ=255)                                               cGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGAGTTGTTGTGa                                                                                              >  1:767125/1‑143 (MQ=255)                                                    gaggaAAATTTCCCGGCCTAGATAATCAGAAAAGCGCGGCAGATATGCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTc                                                                                         <  2:465010/140‑1 (MQ=255)                                                     atgaAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGTGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGTGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCt                                                                                        >  2:894603/1‑143 (MQ=255)                                                          attgacccGGCCTGGATGAGGAGAAAAGCGCGGCAGATAGTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGGGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCttatta                                                                                   <  2:170948/138‑1 (MQ=255)                                                          aaTTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGTGTTAAATTATGCAGTGGCATCACCGCCTCCGATGTGCGTTGTTGTGACTTTCttatta                                                                                   >  2:195250/1‑143 (MQ=255)                                                            tttCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATCACAGCCTCCGCTGTGTCTTGTTGTGATTTTCTTATTTTg                                                                                 >  1:245139/1‑143 (MQ=255)                                                            tttCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTGGTTGTGATTTTCTTATTATg                                                                                 >  2:609623/1‑143 (MQ=255)                                                                   tCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGTGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTg                                                                          >  1:553972/1‑143 (MQ=255)                                                                      tAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTCGTTCTGATTTTCTTACTCTGCCCGCTGaaa                                                                       >  2:809577/1‑143 (MQ=255)                                                                          tAATCAGAAAGGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAAcgc                                                                   <  1:280779/143‑1 (MQ=255)                                                                             tCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGta                                                                >  2:519751/1‑143 (MQ=255)                                                                             tCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGACCGCTGAACACGCGta                                                                >  2:120216/1‑143 (MQ=255)                                                                             tCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTGCTTGTTATGCACGCTGAAAGCGCGta                                                                >  1:532860/1‑143 (MQ=255)                                                                             tCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTGAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGCTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGta                                                                >  2:463725/1‑143 (MQ=255)                                                                              cAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACCGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGtaa                                                               >  2:523525/1‑143 (MQ=255)                                                                               aGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGCTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGtaaa                                                              >  1:413249/1‑143 (MQ=255)                                                                                      ctcgGCAGATATTCGAGCGGCGTTGGCGATGCTATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGATGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATaaaaaa                                                       <  2:167068/141‑1 (MQ=255)                                                                                           cAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGgcgc                                                  <  2:835613/143‑1 (MQ=255)                                                                                                  tCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGCAAACGTGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTCGGTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAACAGGCGCCAGTGaa                                           >  1:102239/1‑143 (MQ=255)                                                                                                       tggcgTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAgcgc                                      <  2:606950/142‑1 (MQ=255)                                                                                                        ggcgTTGGCGCGCCGATAAACTCCAGCCGTGGATCGCGAGTGTATTGATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGcc                                     <  1:605547/143‑1 (MQ=255)                                                                                                              ggCGCGCCGATACACTCCAGCCGTGGAAAACGTGTGAAATCACCCCCCACCATCCCCGCGTCGGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCttttt                               <  1:202592/143‑1 (MQ=255)                                                                                                               gcgcgcCGAGAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCtttttt                              <  1:299056/143‑1 (MQ=255)                                                                                                                 gcgcCGATAAACTCCAGCCGTGGAAAACGGGTTGAATTATGCAGTGGCATAACAGCCTCCAATGTGTGTTGTTGTGCTTTACTTACTATTCACGCTGAAAACGCGTAAATAAACAAGGCGCTAGTGAAAGCGCCCTTTTTTGt                            >  1:79049/1‑143 (MQ=255)                                                                                                                 gcgcCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGt                            >  1:420074/1‑143 (MQ=255)                                                                                                                  cgcCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGAGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTc                           <  2:767125/143‑1 (MQ=255)                                                                                                                     cGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCAtt                        >  1:570899/1‑143 (MQ=255)                                                                                                                                           aCGGGTTAAATTCTGCAGTGGTATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAACCGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATc  <  2:235123/143‑1 (MQ=255)                                                                                                                                           aCGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGCAAGCGCCCTTTTTTGTCATTATGCTGACTATTTCGTCACATc  >  2:443723/1‑143 (MQ=255)                                                                                                                                              |                                                                                                                                            CTTACGTAATTTATTGCCGCCCATTGCCATGGGGGTGACGTCATCCCGTTTGATGAAAATTTCCCGTCCTAGATAATCAGAAAAGCGCGGCAGATATTCGAGCGGCGTTGGCGCGCCGATAAACTCCAGCCGTGGAAAACGGGTTAAATTATGCAGTGGCATAACAGCCTCCGATGTGTGTTGTTGTGATTTTCTTATTATGCACGCTGAAAACGCGTAAATAAAAAAGGCGCTAGTGAAAGCGCCCTTTTTTGTCATTATGCTGATTATTTGGTCACATC  >  NC_000913/1999319‑1999599

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑