Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,763	C→T	N10N (AAC→AAT)	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,763	0	C	T	87.1%	67.7 / 3.8	33	N10N (AAC→AAT)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (6/21);  total (9/24)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.49e-01

TCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGAT  >  NC_000913/3525‑4009                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        tCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAg                                                                                                                                                                                                                                              <  1:89776/250‑1 (MQ=255)                  gcgGTTGCGGAAATCGGCGCGGTAGCGGGCGGTACCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAAc                                                                                                                                                                                                                            <  2:344302/251‑1 (MQ=255)                                                ggTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGCCGAGGACTGGGAAACTAAATGAAAATCTACAATCTGAAAGATAACAAATAGCAGGGACGCTTTTCGCACACCGTAACCCAAGGGGTGGGGGAAAAA‑CAgggggtggt                                                                                                                                                                                             >  1:68115/1‑246 (MQ=255)                                                     ccccggcgccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTc                                                                                                                                                                                        <  2:264586/249‑1 (MQ=255)                                                       ccggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTc                                                                                                                                                                                        <  1:130608/249‑1 (MQ=255)                                                        cggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGtt                                                                                                                                                                                                                  <  1:380208/222‑1 (MQ=255)                                                                gcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTg                                                                                                                                                                             <  1:124718/251‑1 (MQ=255)                                                                         tGTTCGATTTGTGTGACAAGCCGGATTCCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAAtt                                                                                                                                                                     <  1:58581/250‑1 (MQ=255)                                                                                 gtgtgtgACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgac                                                                                                                                                            <  2:237840/250‑1 (MQ=255)                                                                                           aGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGc                                                                                                                                                                            <  1:14997/225‑1 (MQ=255)                                                                                             ccGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTg                                                                                                                                                                                               <  2:264470/204‑1 (MQ=255)                                                                                                      gcccatggggTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAACCAGGAAGGTTTTGTTCATATTTGCCGGCTGGAGACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGACCACAATGGGCAGGTCAGGTGTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTGCCGCGCGACCTGCCGGAATTCAGCCTGACTGAAATTGAGGAGAGGCTGaa                                                                                                                                       <  2:198553/242‑1 (MQ=255)                                                                                                                 tGCCGACTGGTTGGGTAAGAACTACCTGCAAAAACAGGAGGGTTTTGTTCATTTTTGCCGTCTGGATACGGCGGGCGCACGAGTACGGGAAAACTAAATGAAACTCTACAATCGGAAAGACCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGTTTGGTCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCATCCTGACTGAAATTGATGAGATGCTGAAGCTGTAtttt                                                                                                                             <  2:305779/250‑1 (MQ=255)                                                                                                                         ggTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCAccc                                                                                                                      <  2:98969/249‑1 (MQ=255)                                                                                                                                    gtgtgctcttcgggtctggAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGAt                                                                                                         <  1:440499/234‑1 (MQ=255)                                                                                                                                                 aTCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGAAAACTAAATGAAAATCTACAAACTGAAAGATCACAATGAGCAGGGCCGCTTTTCGCAAGCCGTAAACCA‑GGGGGTGGGGAAAAATAAGGGGCTGTTTTTTTCGCCCCACCTGCCGGGATTCAGCCTGACTGAAATTGTTGAGATGCTGAAGCTGGAGTTTTTCTCCCGCCGAGGGGAGAACCTCTCGGCGttt                                                                                            >  1:364465/1‑251 (MQ=255)                                                                                                                                                      gAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACGCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTgg                                                                                       <  2:204517/251‑1 (MQ=255)                                                                                                                                                       aaGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGt                                                                                      <  1:108617/251‑1 (MQ=255)                                                                                                                                                               gTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGaaa                                                                               <  2:364465/250‑1 (MQ=255)                                                                                                                                                                             gcgggATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAAGGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATcc                                                                <  2:68115/248‑1 (MQ=255)                                                                                                                                                                                ggATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTgg                                                             <  1:269836/251‑1 (MQ=255)                                                                                                                                                                                     cggcggGCGCCCGCGTACTTGAAAACTAAATGAAACTCTACAAGCTGGAAGATCACCATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCC‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGACTTAAGCCTGGCTTAAATTTATGAGGTTCTTAAGCTGGAGTTTGCTACCCGCAATGCGAAGAACCTTTCGGCGTTTATTTGTGATGAAATCCCGCAGGacacccagtaaa                                                          >  2:74165/1‑238 (MQ=255)                                                                                                                                                                                      ggcgggcgCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgac                                                                                                                                                            <  1:104515/150‑1 (MQ=255)                                                                                                                                                                                      ggcgggcgCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgac                                                                                                                                                            >  2:104515/1‑150 (MQ=255)                                                                                                                                                                                                agAGTACTCGAAAACTACTTGAATCTCTACAATCTGAAAGATAACAATGAGCAGGGCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATAAGGGGCTGTTTTCTCCGCCAGACATGCCGGACTTCAGCCTGACTGCAATTGATGCGATGCTGACGCTGGCTTTTGTCACCCGATGTGCGAAGTTCCTCTCGGCGTTTATTGGTGATGCAATCCCGCAGGCAATCCTGGAAAAACGCTTgcgcgc                                             >  2:406348/2‑251 (MQ=255)                                                                                                                                                                                                     aCTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGt                                        >  1:433245/1‑251 (MQ=255)                                                                                                                                                                                                        ggAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGttt                                      <  1:256240/250‑1 (MQ=255)                                                                                                                                                                                                            aaCTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCtt                                 <  1:91782/251‑1 (MQ=255)                                                                                                                                                                                                                      aaccctACAATCGGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTTGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCg                       <  1:383597/247‑1 (MQ=255)                                                                                                                                                                                                                          ctACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCg                   >  1:187611/1‑251 (MQ=255)                                                                                                                                                                                                                                        accacaATGAGCAGGCCAGCTTTGCGCTAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGc     <  2:328586/249‑1 (MQ=255)                                                                                                                                                                                                                                           acaATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGAt  >  1:16663/1‑251 (MQ=255)                                                                                                                                                                                                                                            caACGAGCAGGTCAGCGTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGGATTCCGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAGAATCCTCTCGGCGTTTATTGGTGATGAACTCCCACAGGAAATCCTGGAAACGCGCGTGCGCGCGACGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTgcaggcga   >  2:190582/1‑242 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                        TCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCA‑GGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCGAT  >  NC_000913/3525‑4009

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