Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,111,654	C→T	14.7%	Q394* (CAG→TAG)	hisD →	histidinol dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,111,654	0	C	T	14.7%	64.5 / 2.8	27	Q394* (CAG→TAG)	hisD	histidinol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (9/14);  new base T (2/2);  total (11/16)
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.99e-01

GCTGGCGCGATCTTCATGGGCCGCCACACCAGCGAAGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGCGAGTCGCTGACCGCCCACGCCCGCAGCGCTGAATACCGTATTTTGA  >  NC_002947/1111528‑1111780                                                                                                                               |                                                                                                                               gCTGGCGCGATCTTCATGGGCCGCCACACCAGCGAAGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGtt                                                                                                                       <  1:479325/136‑1 (MQ=255)   tGGCGCGATCTTCATGGGCCGCCACACCAGCGAAGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTtcg                                                                                                                     >  2:361294/1‑136 (MQ=255)   tGGCGCGATCTTCATGGGCCGCCACACCAGCGAAGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTtcg                                                                                                                     >  2:282584/1‑136 (MQ=255)      cgcgATCTTCATGGGCCGCCACACCAGCGAAGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCTCTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTtcgtcg                                                                                                                  <  2:268750/136‑1 (MQ=255)      cgcgATCTTCATGGGCCGCCACACCAGCGAAGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTtcgtcg                                                                                                                  <  1:361294/136‑1 (MQ=255)       gcgATCTTCATGGGCCGCCACACCAGCGAAGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTATCGTCGCCGCTGGGGGTGTATGACTTCCAGCAGCGTTCGGCGa                                                                                                                 >  2:304886/1‑136 (MQ=255)                          acacCAGCGAAGCGCTGGGCTATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACa                                                                                              <  1:309460/136‑1 (MQ=255)                                   aaGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCTAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTc                                                                                     >  1:28831/1‑136 (MQ=255)                                   aaGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCTAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTc                                                                                     >  1:31449/1‑136 (MQ=255)                                    aGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGTCTTCCAGAAGCGTTCGTCGATCATCTTCTTCTCCGAACAGGGCGCGTcc                                                                                    <  1:427378/136‑1 (MQ=255)                                    aGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTcc                                                                                    >  1:312694/1‑136 (MQ=255)                                        cTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAAc                                                                                >  2:20559/1‑136 (MQ=255)                                                    tGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTg                                                                    <  1:164718/136‑1 (MQ=255)                                                    tGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTg                                                                    <  1:267317/136‑1 (MQ=255)                                                    tGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTg                                                                    <  1:244560/136‑1 (MQ=255)                                                                     cGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGt                                                   <  1:282584/136‑1 (MQ=255)                                                                     cGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGt                                                   <  1:213123/136‑1 (MQ=255)                                                                        gctgcCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGc                                                >  2:382670/1‑136 (MQ=255)                                                                        gctgcCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGc                                                >  2:456912/1‑136 (MQ=255)                                                                        gctgcCCACTTCAGGTACAGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGc                                                >  1:115021/1‑136 (MQ=255)                                                                                              cgTTTTTCGTCGCCGCTGGGGGTGTATGACTTCTAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGCGAGTCGCTGACCGCCCACGccc                          <  2:28831/136‑1 (MQ=255)                                                                                              cgTTTTTCGTCGCCGCTGGGGGTGTATGACTTCTAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGCGAGTCGCTGACCGCCCACGccc                          <  2:31449/136‑1 (MQ=255)                                                                                              cgTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGCGAGTCGCTGACCGCCCACGccc                          >  1:507658/1‑136 (MQ=255)                                                                                                 ttttCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGCGAGTCGCTGACCGCCCACGCCCGCa                       <  1:382670/136‑1 (MQ=255)                                                                                                                ggggTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGCGAGTCGCTGACCACCCACGCCCGCAGCGCTGAATACCGTa        <  1:20559/136‑1 (MQ=255)                                                                                                                     gtATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGCGAGTCGCTGACCGCCCACGCCCGCAGCGCTGAATACCGTATTTTg   <  1:529822/136‑1 (MQ=255)                                                                                                                      tATGACTTCCAGATGCGTTCGTCGATCTTCTTCTGCTCCGATCTTGGCGCGTCCGAACTGGGCCACACTGCCTGGGTCCTGGCCCGTGGCGAGTCGCTGACCGCCCACGCCCGCAGCGCTGAATACCGTATTTTGa  <  2:127368/136‑1 (MQ=255)                                                                                                                               |                                                                                                                               GCTGGCGCGATCTTCATGGGCCGCCACACCAGCGAAGCGCTGGGCGATTACTGCGCTGGCCCTAACCACGTGCTGCCCACTTCAGGTACCGCGCGTTTTTCGTCGCCGCTGGGGGTGTATGACTTCCAGAAGCGTTCGTCGATCATCTTCTGCTCCGAACAGGGCGCGTCCGAACTGGGCCACACTGCCTCGGTCCTGGCCCGTGGCGAGTCGCTGACCGCCCACGCCCGCAGCGCTGAATACCGTATTTTGA  >  NC_002947/1111528‑1111780

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