Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,317,129	G→A	Q32* (CAG→TAG)	trpA ←	tryptophan synthase subunit alpha

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,317,129	0	G	A	86.7%	35.8 / NA	15	Q32* (CAG→TAG)	trpA	tryptophan synthase subunit alpha
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (7/6);  minor base C (2/0);  total (9/6)
Fisher's exact test for biased strand distribution p-value = 4.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCA  >  NC_000913/1316967‑1317340                                                                                                                                                                   |                                                                                                                                                                                                                    ccAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGCATCGTCGCGCCATCAGCCAGTGGGTCGGAGAAGGGGATACATAAATCCCGCGCGTCAGCACCGGCTTCAATTAGAGTATCGATAATTT‑CCATGACTCCTCAATGCCCGGATCCCCGAGCGTGACGACAGGCACGAACGCCACCTCTTTGCGCTCCTTCAAATGGGCAAACCCAGCTTTGTAGCGtt                                                                                                                            >  1:325181/1‑251 (MQ=255)                             cacaCCTGCCGCAAAGGCGCGCAGCGTGGCGTTTTGACTCGTCGGGCCCTCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCCGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTcc                                                                                                               >  1:413467/1‑237 (MQ=255)                                                         gCGTTTTGAAGCGTCGGGCCAGCCGCCAGTGGGTCGGAGAAGGGGATTCCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCg                                                                                          <  1:77294/230‑1 (MQ=255)                                                             tttGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATcg                                                                  <  1:459807/250‑1 (MQ=255)                                                                          gCCATCCGNCNNNGNNTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGCCTTCNATTAGCGTCTCGATANTTNTCACNGACNACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTgc                                                                                                                                                                >  2:348883/1‑143 (MQ=255)                                                                          gCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACGACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTgc                                                                                                                                                                <  1:348883/143‑1 (MQ=255)                                                                               ccGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCg                                                                                                                                                                                              <  1:10426/108‑1 (MQ=255)                                                                                 gCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCg                                                                                          >  1:213474/1‑206 (MQ=255)                                                                                 gCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAacc                                             <  1:436609/251‑1 (MQ=255)                                                                                                         ccTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTCTCGATAATTTTCAATGCCTCCTCAATGCCCGGATCAc                                                                                                                                                                                                    >  1:79328/1‑76 (MQ=255)                                                                                                                   gcgcGTCAGCACCGGCTTCAAGTAGCGTATCGATAATGTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCACTCCTTCAGCCGGGCAGACAGAGATTCGTAACTTTCCATCAGAGCTCCCCGCGTGCTTTCAAAAGATCGTCAACGGTGAAGACGTCGTTAACGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGGTTTCGCGca           <  1:486143/251‑1 (MQ=255)                                                                                                                             aCCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCa  >  1:219738/1‑250 (MQ=255)                                                                                                                              ccGGCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAacca                                            >  1:204960/1‑207 (MQ=255)                                                                                                                                ggCTTCAATTAGCGTATCGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGcatca        >  1:368157/1‑241 (MQ=255)                                                                                                                                                 cGATAATTTTCAATGACTACTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAAt                                                                                                                                                                             >  1:204396/1‑59 (MQ=255)                                                                                                                                                                   |                                                                                                                                                                                                                    CCAGCATTTCAAAACATTGTGCCGGAGTCACACCTGCCGCAAAGGCGCGCAGAGTGGCGTTTTGAATCGTCGGGCCATCCGCCAGTGGGTCGGAGAAGGGGATACCTAACTCCAGCGCGTCAGCACCGGCTTCAATTAGCGTATCGATAATTTTCAATGACTGCTCAATGCCCGGATCACCGAGCGTGACGAAAGGAACGAATGCGCCTTCTTTGCGCTCCTTCAACTGGGCAAACAGAGATTCGTAGCGTTCCATCAGATTTCCCCTCGTGCTTTCAAAATATCGTGAACGGTGAAGATGTCTTTATCGCCGCGACCGGAAAGGTTAACCACCAGTAGCTGCTCTTTATCCGGGTTTTCGCGCATCATTTTCA  >  NC_000913/1316967‑1317340

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