Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,842,809	(T)5→6	coding (439/1344 nt)	gdhA →	glutamate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,842,804	1	.	T	90.5%	57.0 / 0.1	23	coding (434/1344 nt)	gdhA	glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base T (10/9);  total (12/11)
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.90e-01

TCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCG‑TTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCC  >  NC_000913/1842687‑1842931                                                                                                                       |                                                                                                                                 tcAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGGAGGTGAACTGACGCGTTTTTT‑CCCAGgcgc                                                                                                                   >  1:176022/1‑133 (MQ=255)       tCCTCGGCGTTGAACGAATCTTTAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCtgatga                                                                                                             <  1:193243/133‑1 (MQ=255)         cTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACt                                                                                                           <  1:1157/133‑1 (MQ=255)          tCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTg                                                                                                          >  1:52010/1‑133 (MQ=255)           cGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCCAAGGTGAAATGAAGCGCTTTTT‑GCCAGGCGATGATGACTGa                                                                                                         >  2:193243/1‑133 (MQ=255)                        aCCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCAc                                                                                            >  2:128797/1‑133 (MQ=255)                        aCCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGAAGCG‑TTTTTAGCCAGGCGCTGATGACAGAACTGTATCGCCAc                                                                                            >  2:168909/1‑133 (MQ=255)                                              ctctGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTAACGCCACCTGGGCACGGAGACCGACGTTc                                                                      >  2:218347/1‑133 (MQ=255)                                               tctGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTATATGGCCCCCTGGGCGTGGATGCCGACGTTcc                                                                     >  2:197663/1‑133 (MQ=255)                                                    cGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAg                                                                >  2:340689/1‑133 (MQ=255)                                                     gATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAgg                                                               <  1:358488/133‑1 (MQ=255)                                                      aTGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGt                                                              >  1:429401/1‑133 (MQ=255)                                                      aTGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGt                                                              >  1:433907/1‑133 (MQ=255)                                                           cGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGatat                                                         <  1:340689/133‑1 (MQ=255)                                                                        ggcAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGAAGGGGAAATCGGGATTggtggt                                            >  1:237477/1‑133 (MQ=255)                                                                          cAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCg                                          <  2:285609/133‑1 (MQ=255)                                                                          cAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCg                                          >  2:112896/1‑133 (MQ=255)                                                                              gATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGaa                                      <  1:128797/133‑1 (MQ=255)                                                                                    gATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGc                                <  2:116793/133‑1 (MQ=255)                                                                                              gAAAAAGCGAAGGTGAAGTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCgg                      <  1:112896/133‑1 (MQ=255)                                                                                               aaaaaGCGAAGGTGATGTGATGCGATGTTT‑GCTAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCggg                     <  2:237477/133‑1 (MQ=255)                                                                                                               gTGATGCGTTTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGctc     <  2:429401/133‑1 (MQ=255)                                                                                                                  aTGCG‑TTTTTGGCCATGCGCTGATGACTTAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTcc  <  2:176022/133‑1 (MQ=255)                                                                                                                       |                                                                                                                                 TCAAATTCCTCGGCTTTGAACAAACCTTCAAAAATGCCCTGACTACTCTGCCGATGGGCGGTGGTAAAGGCGGCAGCGATTTCGATCCGAAAGGAAAAAGCGAAGGTGAAGTGATGCG‑TTTTT‑GCCAGGCGCTGATGACTGAACTGTATCGCCACCTGGGCGCGGATACCGACGTTCCGGCAGGTGATATCGGGGTTGGTGGTCGTGAAGTCGGCTTTATGGCGGGGATGATGAAAAAGCTCTCC  >  NC_000913/1842687‑1842931

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: ‑