Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	789,133	T→C	18.0%	L123P (CTG→CCG)	gdhA →	NADP‑dependent glutamate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	789,133	0	T	C	18.0%	56.0 / 4.4	27	L123P (CTG→CCG)	gdhA	NADP‑dependent glutamate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (13/9);  new base C (2/3);  total (15/12)
Fisher's exact test for biased strand distribution p-value = 6.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGCCACATCGGTGCTGACTGCGACGTACCGGCCGGTGACATCGG  >  NC_002947/788985‑789277                                                                                                                                                     |                                                                                                                                                 aaGGTGCAGGTCAACCGCGGCTACAGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTg                                                                                                                                                 >  1:158340/1‑150 (MQ=255)                cgcgGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCCGCCCATGggcagcgag                                                                                                                                  <  2:664969/150‑3 (MQ=255)                    gCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCaacg                                                                                                                              <  2:553712/150‑3 (MQ=255)                        ccGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCCGCCCATGggcagcgag                                                                                                                                  <  1:810285/142‑3 (MQ=255)                        ccGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCCGCCCATGggcagcgag                                                                                                                                  >  2:810285/1‑140 (MQ=255)                              tCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGa                                                                                                                    <  1:512339/150‑1 (MQ=255)                                agagGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGTGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACt                                                                                                                  >  2:213886/1‑150 (MQ=255)                                  aTGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACttt                                                                                                                <  1:372357/150‑1 (MQ=255)                                                  gCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAaga                                                                                                >  2:456370/1‑150 (MQ=255)                                                    ccGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGc                                                                                              <  2:491541/150‑1 (MQ=255)                                                    ccGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCCGCCCATgg                                                                                                                                          <  2:527295/106‑1 (MQ=255)                                                    ccGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCCGCCCATgg                                                                                                                                          >  1:527295/1‑106 (MQ=255)                                                           aGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCAGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCg                                                                                       <  1:348554/150‑1 (MQ=255)                                                                             aCCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGcc                                                                     >  1:639306/1‑150 (MQ=255)                                                                              cccGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCa                                                                    <  2:376457/150‑1 (MQ=255)                                                                              cccGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCa                                                                    >  1:376457/1‑150 (MQ=255)                                                                                     gTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTc                                                             <  1:213886/150‑1 (MQ=255)                                                                                                       aggTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGc                                           <  2:163083/148‑1 (MQ=255)                                                                                                                               ttCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGCCACATCGGTGCTGACTGCGACGTa                   >  2:702119/1‑150 (MQ=255)                                                                                                                                 cAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGCCACATCGGTGCTGACTGCGACGTAcc                 >  2:498147/1‑150 (MQ=255)                                                                                                                                  aaGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGCCACATCGGTGCTGACTGCGACGTAccg                >  1:325867/1‑150 (MQ=255)                                                                                                                                  aaGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGCCACATCGGTGCTGACTGCGACGTAccg                >  1:766161/1‑150 (MQ=255)                                                                                                                                  aaGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGCCACATCGGTGCTGACTGCGACGTAccg                >  1:172938/1‑150 (MQ=255)                                                                                                                                     aaCTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGa                                                                                                                    >  2:402399/1‑47 (MQ=255)                                                                                                                                     aaCTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGa                                                                                                                    <  1:402399/47‑1 (MQ=255)                                                                                                                                          gCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGCCACATCGGTGCTGACTGCGACGTACCGGCCGGTGa        >  2:713588/1‑150 (MQ=255)                                                                                                                                                cTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGCCACATCGGTGCTGACTGCGACGTACCGGCCGGTGACATCgg  >  1:296370/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AAGGTGCAGGTCAACCGCGGCTACCGCATTCAGATGAGCAGTGCCATTGGCCCGTACAAGGGCGGGCTGCGCTTCCACCCGTCGGTGAACCTCAGCGTGCTGAAGTTCCTGGCGTTCGAGCAGGTGTTCAAGAACTCGCTGACCTCGCTGCCCATGGGCGGCGGCAAGGGCGGCTCGGACTTTGACCCGAAAGGCAAGAGCGACGCCGAAGTGATGCGCTTCTGCCAGGCGTTCATGAGCGAGCTGTATCGCCACATCGGTGCTGACTGCGACGTACCGGCCGGTGACATCGG  >  NC_002947/788985‑789277

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