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	90.5%	58.9 / NA	24	N10N (AAC→AAT)	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base T (7/12);  minor base G (0/2);  total (9/15)
Fisher's exact test for biased strand distribution p-value = 5.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

GGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGC‑TCTGTGTGACAAGCCGGAAACCGCCCAGC‑GCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCG  >  NC_000913/3521‑4007                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         ggCTGCCGGCAGGCGCGGCAGGCGGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGC‑TCTGTGTGACAAGCCGGATACCGCCCAGC‑GCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAgg                                                                                                                                                                                                                                                   <  1:194454/250‑1 (MQ=255)       cggcaggcgcggcaggcgGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGC‑TCTGTGTGACAAGCCGGATACCGCCCAGC‑GCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCtt                                                                                                                                                                                                                                            <  1:479116/251‑1 (MQ=255)                     ggcgGTTGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGC‑TCTGTGTGACAAGCCGGATACCGCCCAGC‑GCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAAGGAGCAGGTCAGCTTTGCGCAAGCCGT‑aa                                                                                                                                                                                                                             <  1:170277/251‑1 (MQ=255)                                 aaccggctcgggggcgtgcggttgctccggctcgggCCCGAC‑TTGGTCGCCTCTGGTTGACAAGCCCGATGCCGACC‑GCTGCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAGGATCACAAGGAGCAGGTGAGTTTTGCGCAAGCCGT‑AACCAAGGGGTTgg                                                                                                                                                                                                                 <  2:118650/227‑1 (MQ=255)                                         cgGTAGCGAGCGGTAGCTCCGGCTCCGGCCCGACCTTGTTCGC‑TCTGTGTGACAAGCCGGATACCGCCCAGC‑GCGTTGCCGACTGGTTGGTTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAAt                                                                                                                                                                                                         <  1:284576/251‑1 (MQ=255)                                                                 ccgtcccGGCCTTGTTCGT‑TGTGAGTAACAAGCCGGATACCGCCCAGG‑GCGTTGCCGACGGGTTGGGTAAGAACTACCTGCAAAAGCAGGAAGGTTTTGTCCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACAAAAGGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGCCTG‑TTTTTTCCCCACG‑Ac                                                                                                                                                                               <  2:307419/247‑1 (MQ=255)                                                                     cccGTCCTTGGTCGC‑GCTGGGTGTCAAGCCGGATACCGCGCAGC‑GCGTTGCCGACTGGTGGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCGGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGAGCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGTGCTG‑TTTTTTCCGCACC‑ACCTGc                                                                                                                                                                           <  2:113851/251‑1 (MQ=255)                                                                                 cGC‑TCTGTGTGACAAGCCGGATACCGCCCAGC‑GCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTTCTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGcc                                                                                                                                                              <  1:397058/251‑1 (MQ=255)                                                                                 cGC‑TCTGTGTGACAAGCCGGATACCGCCCAGC‑GCGTTGCCGACTGGTGGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCCCGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGc                                                                                                                                                               <  1:78141/250‑1 (MQ=255)                                                                                                      ataCCGCCCAGC‑GCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGAGGGCGGGGAATGGGAAAATAAAAGAAAATCTACAATCTGAAAGAACACAATGAGCAGGGCAGCTTTTCGCAAGCCGTGAACCCAGGGGTGGGGGAAAAA‑CAGGGGGTG‑GTTTTTTCGCCCGGCCCCCCCGG‑AATT‑AGCCTGACTGAAAatgatgtggtt                                                                                                                                          >  1:63813/3‑246 (MQ=255)                                                                                                        aCCGCCCAGC‑GCGTTGCCGACTGGTTGGGTAAGAACTACCGGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGc                                                                                                                                         <  2:287285/250‑1 (MQ=255)                                                                                                                 gcgggTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTgg                                                                                                                                <  2:149876/247‑1 (MQ=255)                                                                                                                            aCTGGTTGGGTAAGAACTACCTGCAACATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGGGCACGGGGACTGGGAAAAAAAAAAAAACTCTTCCAACTGAAAGAACACAAAGAGCAGGGCAGCTTTGCGCAAACCCG‑AACCCCGGGGGTGGGCAAAAATTAGGGGGTGGTTTTTTCC‑CCCG‑CCCCGCCCG‑GATTCAGCCTGGCTGAAATTGATGAGAGGGTGGAGCTGGGtttttttacc                                                                                                                     >  1:153623/1‑245 (MQ=255)                                                                                                                                   ggttAAGATCTACCTTCAAAATCAGGAAGGTTTTTTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCTGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTg                                                                                                              <  2:266239/248‑1 (MQ=255)                                                                                                                                                               ggTTTTGTTCATATTTGCCGGCTGGATACGGCGGGAGCAAGCGGACTGGGAAAATAAATGGAAATCTACAAACTGAAAGAACACAAAGAGCAGGTCCGCTTTGCGCAAACCGT‑AACCCCGGGGGTTGGCGAAAAAAAGGGGCTG‑TTTTTTTCCCACG‑ACCCGCCGGGAATT‑AGCCTGACTGAAATTTGTGGGATGCTGAAACTTGGTTTTGTCACCCCCCGGGCGGAGGTCCTCTCGGCGTTTATTgggtg                                                                                  >  1:176566/1‑248 (MQ=255)                                                                                                                                                                                                          tACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGCAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCgcggcg                                       >  1:89309/1‑251 (MQ=255)                                                                                                                                                                                                                  aaCTAAATTAAACTTTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCGG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCtt                               <  1:251027/251‑1 (MQ=255)                                                                                                                                                                                                                          gAAACTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTc                       >  1:112126/1‑251 (MQ=255)                                                                                                                                                                                                                            aaCTCTACAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCg                     <  1:337189/251‑1 (MQ=255)                                                                                                                                                                                                                                  aCAATCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCCGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGcc               >  1:474427/1‑251 (MQ=255)                                                                                                                                                                                                                                     aTCTGAAAGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAAt            >  1:100661/1‑251 (MQ=255)                                                                                                                                                                                                                                      ccggaaagacctaaaCGAGCAGGCCCGCTTTGCGCAAGCCGT‑AACCCAGGGGTCGGGCAAAAATCCGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGCGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCGCGGCTTTGAGTGGTGATGAAAGCCCGCAGGAAATCCTGGAAGATCGCGTGCGCGCGGCGTTTGCCTTCCCCGCTCCGGTCGCCattg           <  2:68236/238‑4 (MQ=255)                                                                                                                                                                                                                                           aaGATCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGa       >  2:377908/1‑250 (MQ=255)                                                                                                                                                                                                                                               tCACAATGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCGCAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAGGCg  >  1:29602/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         GGCTTCCGGCAGGCGCGGCAGGCGGTCGCGGAAATCGGCGCGGTAGCGAGCGGTATCTCCGGCTCCGGCCCGACCTTGTTCGC‑TCTGTGTGACAAGCCGGAAACCGCCCAGC‑GCGTTGCCGACTGGTTGGGTAAGAACTACCTGCAAAATCAGGAAGGTTTTGTTCATATTTGCCGGCTGGATACGGCGGGCGCACGAGTACTGGAAAACTAAATGAAACTCTACAATCTGAAAGATCACAACGAGCAGGTCAGCTTTGCGCAAGCCGT‑AACCCAGGGGTTGGGCAAAAATCAGGGGCTG‑TTTTTTCCGCACG‑ACCTGCCGG‑AATTCAGCCTGACTGAAATTGATGAGATGCTGAAGCTGGATTTTGTCACCCGCAGTGCGAAGATCCTCTCGGCGTTTATTGGTGATGAAATCCCACAGGAAATCCTGGAAGAGCGCGTGCGCGCGGCGTTTGCCTTCCCGGCTCCGGTCGCCAATGTTGAAAGCG  >  NC_000913/3521‑4007

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