Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,821,424	G→C	G196G (GGC→GGG)	alaS ←	alanine‑‑tRNA ligase/DNA‑binding transcriptional repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,821,424	0	G	C	90.0%	59.7 / 0.2	20	G196G (GGC→GGG)	alaS	alanine‑‑tRNA ligase/DNA‑binding transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (10/8);  total (11/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.05e-01

TACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACG  >  NC_000913/2821297‑2821555                                                                                                                                |                                                                                                                                    tACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAAtg                                                                                                                         >  1:527938/1‑140 (MQ=255) tACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGCGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAAtg                                                                                                                         >  2:197649/1‑140 (MQ=255)   ccGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGGGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAAtgtg                                                                                                                       <  2:535210/140‑1 (MQ=255)                aCGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGcc                                                                                                                 >  2:145154/1‑133 (MQ=255)                aCGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGcc                                                                                                                 <  1:145154/133‑1 (MQ=255)                 cGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGt                                                                                                         >  1:322002/1‑140 (MQ=255)                                   cATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACg                                                                                       <  2:527938/140‑1 (MQ=255)                                          ccATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCa                                                                                <  2:113312/140‑1 (MQ=255)                                          ccATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCa                                                                                >  1:113312/1‑140 (MQ=255)                                          ccATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCa                                                                                <  1:130212/140‑1 (MQ=255)                                               ggCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGAc                                                                           <  2:44082/140‑1 (MQ=255)                                                 ccTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCa                                                                         <  2:322002/140‑1 (MQ=255)                                                                                   gATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGt                                       >  1:163608/1‑140 (MQ=255)                                                                                      ctcAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATg                                    >  2:497896/1‑140 (MQ=255)                                                                                         aaTGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGA‑GGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGgcgc                                >  2:370647/1‑140 (MQ=255)                                                                                                cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCg                                        >  1:260815/1‑126 (MQ=255)                                                                                                cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCg                                        <  2:260815/126‑1 (MQ=255)                                                                                                cGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCttt                          >  2:210010/1‑140 (MQ=255)                                                                                                            tcttcCGGGCTTCCCGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGAt              <  2:163608/140‑1 (MQ=255)                                                                                                                  gggCTTCCCGCA‑GGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGaata       <  1:370647/140‑1 (MQ=255)                                                                                                                        cccGGAGGCCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACg  >  1:537420/1‑140 (MQ=255)                                                                                                                                |                                                                                                                                    TACCGGTATCTACAGACGGCTTCGGCAGCGGTTCCATCGTGCCATCGGCCTGGCGGTTGAACTGCATGAAGACGATGTTCCAGATCTCAATGTAGCGGTCGCCGTCTTCTTCCGGGCTTCCCGGAGGGCCCCCCCAAATGTGGTCGCCGTGATCGTAGAAGATTTCGGTGCACGGGCCGCACGGACCAGTGTCACCCATCTGCCAGAAGTTGTCAGATGCGTATGGCGCACCTTTGTTATCGCCGATGCGAATAATACG  >  NC_000913/2821297‑2821555

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