Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	635,896	A→G	L225P (CTC→CCC)	ybdN ←	PAPS reductase‑like domain protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	635,896	0	A	G	95.5%	68.2 / ‑3.9	22	L225P (CTC→CCC)	ybdN	PAPS reductase‑like domain protein
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (10/11);  total (10/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

CGATCCGGTTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAA  >  NC_000913/635770‑636029                                                                                                                               |                                                                                                                                      cGATCCGGTTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTgatgat                                                                                                                           >  2:57867/1‑139 (MQ=255)  gATCCGGTTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGTCTTtgaggttt                                                                                                                          <  1:337321/139‑6 (MQ=255)           cGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCa                                                                                                                 >  2:320829/1‑139 (MQ=255)                                 tCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTc                                                                                           <  2:386699/139‑1 (MQ=255)                                           tGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTaaa                                                                                 <  1:466435/139‑1 (MQ=255)                                            gCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAAt                                                                                >  1:131235/1‑139 (MQ=255)                                              tCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATgg                                                                              <  1:57867/139‑1 (MQ=255)                                                        aaGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCTCTTTCCAGTCGTAAATGGGGTAAATGTa                                                                    <  1:320829/139‑1 (MQ=255)                                                        aaGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTCGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa                                                                    <  1:413053/139‑1 (MQ=255)                                                                                  agagGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGGTACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCg                                          <  2:452992/139‑1 (MQ=255)                                                                                                                  ctgTGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATTGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACgttgt          <  1:386497/137‑1 (MQ=255)                                                                                                                     tGGGTTGCAGGGACTTTGCTGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCATGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGttta       >  1:15852/1‑139 (MQ=255)                                                                                                                     tGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTc                                                                                           <  1:372201/55‑1 (MQ=255)                                                                                                                     tGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTc                                                                                           >  2:372201/1‑55 (MQ=255)                                                                                                                     tGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTa                                                                          <  1:47479/72‑1 (MQ=255)                                                                                                                     tGGGTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTa                                                                          >  2:47479/1‑72 (MQ=255)                                                                                                                      gggTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa                                                                    >  2:351424/1‑77 (MQ=255)                                                                                                                      gggTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa                                                                    <  1:491059/77‑1 (MQ=255)                                                                                                                      gggTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa                                                                    <  1:351424/77‑1 (MQ=255)                                                                                                                      gggTTGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTa                                                                    >  2:491059/1‑77 (MQ=255)                                                                                                                          tGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTAt                                                            >  2:123267/1‑81 (MQ=255)                                                                                                                          tGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTAt                                                            <  1:123267/81‑1 (MQ=255)                                                                                                                          tGCAGGGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTaa  >  1:126747/1‑139 (MQ=255)                                                                                                                               |                                                                                                                                      CGATCCGGTTCGATAACGTGATAGAGCCATAATCCTTGTCGTTGCTCCGGGCCAAAAGGTTCGCAAATTCGCATATGACGCAGAGGAACGCCTGCCTGATACATTAGGTTATACAGTGGGTTGCAGAGACTTTGATGATTTGCATACCAGGTCCAGATATCAGCCACTTTCCAGTCGTAAATGGGGTAAATGTACCAGCTATGACCGCCTGGTGCGGCCGTAGTCCAGGGTTTATCGTCGGCAAAACGTTGTTTATTTAA  >  NC_000913/635770‑636029

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