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	88.9%	56.6 / ‑1.1	29	N10N (AAC→AAT)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (5/19);  total (8/21)
Fisher's exact test for biased strand distribution p-value = 5.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.72e-01

GCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGC  >  NC_000913/3517‑4006                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       gcctggCTTCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAg                                                                                                                                                                                                                                                     <  1:18553/246‑1 (MQ=255)     ggttcccGGCAGGCGCGGCAGGCGGTTGCTGAAATCGGCGCGGTGGCGAGCGGTATCTCCGGCGCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGTCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAAGCAGGAAGGTTTTGGTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGTAAAACTAAATGAAACTCTACAATCTGAATGATCACAATGAGCAGGt                                                                                                                                                                                                                                                <  2:196253/246‑1 (MQ=255)               cggcgcggcagtagGTTGCTGAAATCGGCGCGGGAGCGATCGGTATCGCCGTCTCCGGCCCGACCTTGGTCGCTCTGTGTGCCAAGCCGGATACCGCCCATTGCGTTTCCGACTGGTTGGGTAAGAACTACCTTCAACATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTAATGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTgcg                                                                                                                                                                                                                                      <  2:302839/250‑1 (MQ=255)                             gttgcgGCAATCGGCCCGGTAGCGAGGGTTACCCCCGGCTCCGGCCCGACATGGTTTGTTCTGTGTGACAATCCGGATACCGCTCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACGAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCa                                                                                                                                                                                                                       <  2:4234/248‑1 (MQ=255)                                        cGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCa                                                                                                                                                                                                            <  1:208635/251‑1 (MQ=255)                                        cGGCGCGGTAGCGAGCGGTAGCTCCGGCTCCTGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGTCCAGCGCGTTGCCGACTGGTTGGGTAATAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTCCAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCa                                                                                                                                                                                                            <  1:48066/251‑1 (MQ=255)                                                               ccggctccggcCCGACCTTGTTCGCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTc                                                                                                                                                                                       <  1:68650/249‑1 (MQ=255)                                                               ccggctccggcACGGCCTTGTTCTCTCTGTGTGACAAGCCGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGttttt                                                                                                                                                                                         <  1:217421/247‑1 (MQ=255)                                                                                                      cGGATACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGACGAGGGCTGGGAAACTAAAAGAAACTCTACCAACTTAAAGAGTACAAAAAACAGGTGCGCTTTTCGCAAGACGTCAACCCAGGGGTGGGGGCACAA‑CACGGGGTGTTTTTTCCGCCCGCCCTGCCGGAAAT‑CAGcccgtcttgaatttttggg                                                                                                                                             >  1:65061/1‑234 (MQ=255)                                                                                                       ggtacccGCCCAGCGAGTTGCCGACTTGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTGAATATTGGCCGGCTGGATACGGCGGGCGTACGATTACTGGAAAATGAAATGAAACGATACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGTTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCTTGACTGAAATTGATgaga                                                                                                                                            <  2:40322/246‑1 (MQ=255)                                                                                                              gcacagggcgTTGACGACTTGTTGAGGAAGAACTACCGGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATAGGGCGGGCGCACGAGTACTGGAAAACTAAATGAACCTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGaa                                                                                                                                     <  2:122527/244‑1 (MQ=255)                                                                                                                  aGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTg                                                                                                                                 <  1:153611/251‑1 (MQ=255)                                                                                                                         tGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGAtttt                                                                                                                           <  1:222469/250‑1 (MQ=255)                                                                                                                          gCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGt                                                                                                                                                                                                                              <  1:178771/152‑1 (MQ=255)                                                                                                                           ccGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTc                                                                                                                        <  1:353408/251‑1 (MQ=255)                                                                                                                             gACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGGCCGAGGCCTTGCAAACTAACTGAAAATCTACAATTTGAAAGAATACAAATAGCAGGGCAGCTTTTTGCCAACCGT‑AAACCCGCGGGTGGGCAAAAATCAGGGGCTGTTTTTTTCGCCCCAACTGCCGCGGTT‑CAGCCTTGCAGATATTGATGAGATTCTTAAGCTGGAtttttttcc                                                                                                                      >  1:176597/1‑246 (MQ=255)                                                                                                                                gtgTTGGGTAAGAACGACCTGCAAACACAGGAAGGGTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACGGGAAATCTAAATGAAACTCTCCAACCTGAAAGATCTCAATTAGCAGGGCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGTACGACCTGCCGGAATT‑CATCCTGACTGAAATTGATGAGATGTGGAAGCTGGATTTTGTCACCCg                                                                                                                   <  2:370125/249‑1 (MQ=255)                                                                                                                                          aGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAg                                                                                                         <  1:77237/251‑1 (MQ=255)                                                                                                                                           aaaCTACATGCAAAATCAGGAAGGTTGTGTTAATATTTGCCGGCTGGAAACGGCGGGCGCACGAGTACTGGAAAAATAAATGAAACTCTACAATCTGAAAGAGCACAATGAGAAGGGCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGAACTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGa                                                                                                        <  2:340357/250‑1 (MQ=255)                                                                                                                                                             ggAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACCCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGAGTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTg                                                                                      <  2:176597/251‑1 (MQ=255)                                                                                                                                                                                   cGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAAt                                                                <  1:194377/251‑1 (MQ=255)                                                                                                                                                                                                     caaccgtgaatGGAAATCAAAAAGAAACTCTAACTTCTGAACGAACACAACGATAAGGTCATCTTTGCGCTAGCCGT‑AACCCAGTGGTTGGGTAAAACTCAGGGTCTGTTGTTTCCTCACGTTCTGCCGTAATT‑CAGCTTGAATGAAATTGATGATAGGCTGAAGCTGGATTTTGTCACCCGCAGTGCTAAGATCCGCTCGGCGTGTATTGGTTATGAAATCTCTCAGTAATTCCTGAAATAGCGCGTgcg                                              <  2:253234/241‑1 (MQ=255)                                                                                                                                                                                                                ggAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACCGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTg                                  >  1:368723/1‑251 (MQ=255)                                                                                                                                                                                                                             aaaCTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTcc                     >  1:209003/1‑251 (MQ=255)                                                                                                                                                                                                                                ctctACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGATGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGCATT‑CAGCCTGACTGAAATTGATGAGATGCTGCAGCTGGATTTTGTCACCCGCAGCGCGAAGATACTCTCGGCGTTTATTGGTGATGCAATCCCGCAGGAAATCATGGCAGAGCGCGTGCGCGC‑GGCATTTGCCTTCACGGCACCGGt                  >  2:337411/1‑251 (MQ=255)                                                                                                                                                                                                                                  ccaCAATCTGAAAGATCACAATGAGCAGGGCAGCGTTGGGCAAGCCGT‑AGCCCAGGGGTTGGGCAAAAAGCAGGGGCGGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCATTCCCGGCTCCGGt                  <  2:144861/247‑1 (MQ=255)                                                                                                                                                                                                                                               gATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGaa     >  1:149506/1‑249 (MQ=255)                                                                                                                                                                                                                                                aTCACAATGTGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGTGGATGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGACGGAATTACA‑CCTGACTGAAATTGCAGAAATGCTGAAGCTGGATTTTGTCAACCGCCGTGCGCAGATCCTCTCGGCGTTTATTAGTGATGAAATCCCGTAGGAAATCCTGCATGAGCCCCTGAGCGA‑GGCGTTTGCCTTGCGGGCTCCGTTCGCCaaggtggcaggc  >  2:287415/1‑241 (MQ=255)                                                                                                                                                                                                                                                aTCACAATGAGAAGGTCCGATTTGCGCAAGCCGT‑ACCACAGGGGTTAGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTACCGGAATT‑CAGACTGAATGAAATGAATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGGGCGAAGATCCTCTCTGCGTTTATTGGTGGTGAAATCACGCGGGAAATCCTGGAAAAGCGCGTGC‑CGCTGGCGTTTCCATTCCCGACTCCGGTCGCCAATggtgcaagg  >  2:287844/1‑243 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       GCCAGGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGCTCTGTGTGACAAGCCGGAAACCGCCCAGCGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTGTTTTTTCCGCACGACCTGCCGGAATT‑CAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGC‑GGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGC  >  NC_000913/3517‑4006

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