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	93.3%	78.4 / ‑6.3	30	N10N (AAC→AAT)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (6/22);  minor base A (2/0);  total (8/22)
Fisher's exact test for biased strand distribution p-value = 6.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATG  >  NC_000913/3517‑3998                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             gcctggCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGTAAGGTTTTGTTCATATTTGCCGGCTGTATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGc                                                                                                                                                                                                                                          <  2:121885/247‑1 (MQ=255) gcctggCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAg                                                                                                                                                                                                                                           <  1:262801/246‑1 (MQ=255)               aggcgcggcaggcgGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTgcc                                                                                                                                                                                                                            <  2:1192/251‑2 (MQ=255)                     ggcaggcgGTTGCGGAAACCGGCGCGGTAGCTAGCGGTATCTCCGGTGCCTGCCCGACCTTGTTCGCTCTGTGTGACAAGCCAGATACCGCAAAGCGCGTTGCCGACTGGTTGGTTAAGAACTACCTGCAAAATCTGGAAGCTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCTCGAGTACTGGAAAACTAAATGAAACTCTACAATCTTAAAGATCACAATGAGCAGGTCAGCTTTGCGTAAGcc                                                                                                                                                                                                                      <  1:267657/251‑1 (MQ=255)                                                   cgcgcgGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTCCCTGCAAAAGCAGGAAGGTTGTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTCCAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAgggg                                                                                                                                                                                       <  1:80635/248‑1 (MQ=255)                                                    gagcgGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCGACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAgg                                                                                                                                                                                         <  1:122757/248‑1 (MQ=255)                                                               ccggctccggcCCGACCTTGTGCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTGTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGAGCACAATGAGCAGGTCAGCTGTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTcc                                                                                                                                                                            <  2:28358/250‑1 (MQ=255)                                                                          ccGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGcc                                                                                                                                                                <  2:171742/251‑1 (MQ=255)                                                                                 tGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTGGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAAGCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAAt                                                                                                                                                           <  1:60686/249‑1 (MQ=255)                                                                                     cGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGcc                                                                                                                                                     <  1:127454/251‑1 (MQ=255)                                                                                        tctGTGTGACAAGCCGGCTACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGGACTGGGAAACTAAATGAAAATCTCAAATCTGAAAAATCACAAAGAGCAGGGCAGCTTTGCGCCAACCGTAACCCAGGGG‑GTGGGGAAAAAACAAGGGGTGTTTTTTTCGCACGAACTGCCGGGATTCAGActg                                                                                                                                                   >  1:75136/1‑250 (MQ=255)                                                                                           gtgtgACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCctgactg                                                                                                                                               <  2:76589/251‑1 (MQ=255)                                                                                                                       gTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGTGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCGGAAGCTGGAttt                                                                                                                   <  1:165048/251‑1 (MQ=255)                                                                                                                              aCTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCAcc                                                                                                            <  1:82147/251‑1 (MQ=255)                                                                                                                                 ggTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTGTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGGCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCg                                                                                                          <  2:131961/250‑1 (MQ=255)                                                                                                                                           gAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTGTGTCACCCGCAGTGCGAAGa                                                                                               <  1:64144/251‑1 (MQ=255)                                                                                                                                                cccTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTTGATACGGCGGGCGCACGATTACTGGAAAACTAAATGAATCTCTACAATCTGAAAGATCACAATGATCAGGTCAGCTTTGCGCAAGCCTTAACCCATGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGTCCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCct                                                                                           <  1:237500/249‑1 (MQ=255)                                                                                                                                                        aaTCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGTGCACGAATTCTTGGTAACTAAATAAAAATCTCTAATCTGGAAGAACACACATAGCAGGGCAGCTTTTCGGAAGACGTAACCCAGGGGGTGGGGCAAAACAAAGGGG‑GGTTTTTTTCGCCCGACCTGCCGGGATTAAGCCTGTCTGTAATTTTTGAGGTGGTGAAGCTGGGGTTTGTCGCCCGCCGTGGGAGGATCCTCTCGGCGtt                                                                                  >  1:32529/1‑251 (MQ=255)                                                                                                                                                           cAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCTAAGATCCTCTCGGCGTTTAt                                                                               <  1:231210/251‑1 (MQ=255)                                                                                                                                                                                                      cACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTgcgc                                    <  1:7614/251‑1 (MQ=255)                                                                                                                                                                                                            tACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGACGATCCTCTCGGCGTTTCTTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCgcggcg                              >  2:19539/1‑251 (MQ=255)                                                                                                                                                                                                                 gAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGc                         >  1:254143/1‑251 (MQ=255)                                                                                                                                                                                                                   aaaCTAAATGAAACTCCACAATCTGAAAGATCACAATGAGCAGGTCAGCTGTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCTCGTGCGCGCGGCGTTTGCCt                       <  1:198353/251‑1 (MQ=255)                                                                                                                                                                                                                     gcTAAATGAAACTCTCCAATCTGCAAGATCACACTGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGACCTCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTc                     >  2:37719/2‑251 (MQ=255)                                                                                                                                                                                                                        aaaTGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTccc                   <  1:181509/250‑1 (MQ=255)                                                                                                                                                                                                                        aaaTGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTccc                   <  2:255985/250‑1 (MQ=255)                                                                                                                                                                                                                                  ctACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTCCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAAAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCg        >  2:193333/1‑251 (MQ=255)                                                                                                                                                                                                                                   tACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGc       >  1:1534/1‑251 (MQ=255)                                                                                                                                                                                                                                     cAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGcc      >  2:62384/1‑250 (MQ=255)                                                                                                                                                                                                                                           gAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTACCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAAGCCCGCAGGAAATCCTGGAATAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATg  <  2:110820/248‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                             GCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGTAACCCAGGGG‑TTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATG  >  NC_000913/3517‑3998

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