Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,046,876	T→C	G722G (GGA→GGG)	gcvP ←	glycine decarboxylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,046,876	0	T	C	96.7%	99.1 / NA	30	G722G (GGA→GGG)	gcvP	glycine decarboxylase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (13/16);  minor base A (1/0);  total (14/16)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGATTGGCAGGATAGAGGCGCTACCGAACGGTGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGTCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAAC  >  NC_000913/3046745‑3047010                                                                                                                                    |                                                                                                                                       cTGATTGGCAGGATAGAGGCGCTACCGAACGGTGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATAcc                                                                                                                                <  1:365762/140‑1 (MQ=255)                      cTACCGAACGGTGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCgg                                                                                                           >  1:150522/1‑140 (MQ=255)                        aCCGAACGGTGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGACCCATACCAGGACCACCACCGCCGTGCGGaa                                                                                                         >  2:263679/1‑140 (MQ=255)                             aCGGTGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAg                                                                                                    <  1:116638/140‑1 (MQ=255)                              cGGTGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATTCCCGGACCACCACCGCCGTGCGGAATGCAGa                                                                                                   <  1:613404/140‑1 (MQ=255)                              cGGTGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGa                                                                                                   <  2:537494/140‑1 (MQ=255)                                 taccGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAg                                                                                                >  2:76546/3‑140 (MQ=255)                                 tGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCGCCACACCGGGACCACCACCGCGGTGCGGAATGCAGAAAg                                                                                                >  2:208823/1‑140 (MQ=255)                                 tGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAg                                                                                                >  1:162511/1‑140 (MQ=255)                                 tGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAg                                                                                                >  1:9379/1‑140 (MQ=255)                                 tGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCCGACCACCACCGCCGTGCGGAATGCAGAAAg                                                                                                >  2:98776/1‑140 (MQ=255)                                       aGAAACCGCGCCCTGTCGGGTTGACTTGCCTTCGATTAGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTAt                                                                                          <  1:208823/140‑1 (MQ=255)                                       aGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTAt                                                                                          <  1:141140/140‑1 (MQ=255)                                       aGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTAt                                                                                          <  2:493604/140‑1 (MQ=255)                                       aGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTAt                                                                                          <  2:12523/140‑1 (MQ=255)                                              gcgcCCTGACGGGTTAACATGCGTTCGATTTGCCCCACGTTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGtt                                                                                   <  2:250762/140‑1 (MQ=255)                                              gcgcCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGtt                                                                                   <  2:391499/140‑1 (MQ=255)                                               cgcCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTa                                                                                  <  2:31418/140‑1 (MQ=255)                                                gcCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATTCAGAAACTTTTATGTAGGTTaa                                                                                 >  1:319600/1‑140 (MQ=255)                                                     gACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGAGCGTAATGCAGACAGTTTTATGTAGGTTAAGgtgt                                                                            >  2:220996/1‑140 (MQ=255)                                                             aaCATGCCTTCGATTGGCACCACGCTATGCCCCGGTACAACCGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTc                                                                    <  2:645694/140‑1 (MQ=255)                                                              aCATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCCCCGTGCGGTATGCAGAAAGTTTTATGTAGGGTAAGGTGTGAAACGTcc                                                                   >  2:275165/1‑140 (MQ=255)                                                                       cGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATa                                                          >  1:645694/1‑140 (MQ=255)                                                                       cGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATa                                                          >  1:554479/1‑140 (MQ=255)                                                                              caccacGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGAAGGTTAAAGTGTGAGACGTCCACACCAATAAAGCCCg                                                   >  1:226883/1‑140 (MQ=255)                                                                                ccacGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGc                                                 <  1:76546/140‑1 (MQ=255)                                                                                ccacGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCCGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGc                                                 <  1:98776/140‑1 (MQ=255)                                                                                          ccccGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGc                                       <  2:150522/139‑1 (MQ=255)                                                                                           cccGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGcc                                      <  2:394066/140‑1 (MQ=255)                                                                                                                               aTCGGCCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAAc  >  2:633425/1‑140 (MQ=255)                                                                                                                                    |                                                                                                                                       CTGATTGGCAGGATAGAGGCGCTACCGAACGGTGCCGCAGAAACCGCGCCCTGACGGGTTAACATGCCTTCGATTTGCACCACGCTATGACCCGGTACAAACGGTGCCAAATGCGCTTTCACGCCGATCGGTCCCATACCCGGACCACCACCGCCGTGCGGAATGCAGAAAGTTTTATGTAGGTTAAGGTGTGAAACGTCCGCACCAATAAAGCCCGGCGAGGTGATGCCAACCTGGGCGTTCATGTTCGCGCCATCAAGGTAAAC  >  NC_000913/3046745‑3047010

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