Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,617,499	T→A	T59T (ACT→ACA)	deoC →	deoxyribose‑phosphate aldolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,617,499	0	T	A	96.7%	77.9 / NA	30	T59T (ACT→ACA)	deoC	deoxyribose‑phosphate aldolase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (18/11);  minor base C (1/0);  total (19/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAA  >  NC_000913/4617253‑4617745                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        gCGGATTAGATGTTAAAATAATAACAAACTCG‑AAAGTGAAATTTTATGGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCC‑CCCGG‑AATGACGACGACCCCGACGAGAAAGTAATTGCTCTGTGACATCAGGCAAAACCCCCGG‑CCGGCAATTCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGa                                                                                                                                                                                                                                                     <  1:207878/251‑1 (MQ=255)      gTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATAACTGATCTGAAAGCACGCAGCCTCCGTGCACTGAAATTGCTAGACCT‑ACCAACCCTG‑AATGCCGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑CCGGCAATACCGCCGCCATCTGAATCTATCCTCGCTTTATCCCGATTGCTCCCAAACCCCTGaaaagg                                                                                                                                                                                                                                                >  1:352849/1‑248 (MQ=255)       tAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGc                                                                                                                                                                                                                                               >  1:220891/1‑251 (MQ=255)           tGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGATTGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGCCGACACCGACGAGAAAGTAATTGCCCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGAAggg                                                                                                                                                                                                                                           >  1:326711/1‑251 (MQ=255)                gAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCAccc                                                                                                                                                                                                                                      >  1:322283/1‑251 (MQ=255)                  aTACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCC‑CCCTG‑AATGACGACGACACCGACGAGATAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCACCGCTATCTGTATCTATCCTCGCTTTATCACGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCg                                                                                                                                                                                                                                    >  1:372026/1‑251 (MQ=255)                  aTACTAACAAACTCGCAAGGTGAA‑TTTTATTGGAGACAAGCCAGGAGAATGAAATGCCTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATAGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCTACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAAACCCGG‑TCGGCAGTACCGCCGCTATCTGTATCTATCCTAGCTTGATCCCGATTGCTCGAAAAACACTGAAAGGGCAGGGCAccc                                                                                                                                                                                                                                      >  2:211004/1‑249 (MQ=255)                  aTACTAACAAAATCGCAAGGTGAA‑TTTTATTGGAGACAAGCCAGGAGAAGGAAATGACTGATCTGAAAGAAAGCAGCCTGCGAGAAATGAAATTGAAGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGACATCAGGCCAAAACCCAGG‑TCGGCAATAACGCCGCTATCTGTATCAATCCTAGCTATATCACGATTGCTCGCAACATACTGAAAGAGCAGGGAAccc                                                                                                                                                                                                                                      >  1:297125/1‑249 (MQ=255)                   tACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCACAACCCCGG‑ACGGCAATACCGCCGCTATATGTCTCTATCCTCGCTTTATCCCGGTTGCTCGCAAAACACTGAAAGAGCCGGGCAcccc                                                                                                                                                                                                                                     >  1:262782/1‑249 (MQ=255)                          aaaCTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCg                                                                                                                                                                                                                            >  1:14706/1‑251 (MQ=255)                                           ttttATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAAc                                                                                                                                                                                                             >  1:65004/1‑250 (MQ=255)                                                 gggcaaCAAGCCAGGAAAATGAAATGACTGATCTGAAAGAAAGCAGCCTGCGTGCACTAAAATTGATGGACCTGCCCA‑CCCGG‑AATGACGACGACACCGACGAAAAAGTAATTGCTCTGTGTCATCAGGCAAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTc                                                                                                                                                                                                      <  1:278295/246‑1 (MQ=255)                                                       aaaGCAAGGAGAATGAATTATCTGATCTAAAAGCAAGAAGCTTGCGGGCACTGAAATTGATGGACCTGACC‑‑CCCTGTAATGACGACGACCCCGACGAGAAAGTAATTGCTCGGTGTCACCAGCCCAAACCCCCGG‑GCGACAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCacac                                                                                                                                                                                                <  1:272010/250‑1 (MQ=255)                                                                agaATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAacgat                                                                                                                                                                                       >  1:120141/1‑250 (MQ=255)                                                                    tgaaatgaCTGATCTGAAAGCAAGCAGACTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTGCCCACACGGTAacgatgcc                                                                                                                                                                                    >  1:82832/1‑246 (MQ=255)                                                                                 cTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCgcgc                                                                                                                                                                       >  1:366178/1‑250 (MQ=255)                                                                                                   cGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGc                                                                                                                                                    >  2:112977/1‑251 (MQ=255)                                                                                                       cACTGAAATTGAGGGCCCTGCCCC‑CCCTG‑AAGGACGACGACACCGACGAGAATGTAATTGCTTTGTTTCTCCGGGCCAAAACCCCGG‑TCGGAAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATc                                                                                                                                                <  1:221741/251‑1 (MQ=255)                                                                                                                       ccatgcccccccTG‑AATGACGACGACCCCGACGAAAAAGTAATTGCTCTGTGTCATCGGGCCAAACCCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCcgat                                                                                                                                 <  1:112977/242‑1 (MQ=255)                                                                                                                                                  cccGACGATAAAGTAATTGCGCTGTGACATCAGCCCAAACCCCCGG‑TCGCAAATCCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACGTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCTTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACcg                                                                                                       <  1:152271/250‑1 (MQ=255)                                                                                                                                                    cgacgaGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACcgcg                                                                                                     >  1:271164/1‑251 (MQ=255)                                                                                                                                                       cgaGAAAGTAATTGCTCTGTGTCATCAGGCCAAAACCCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCgcgcgc                                                                                                  >  1:344025/1‑251 (MQ=255)                                                                                                                                                                              caccagtcAAAACCCCGGGTCGGCAATTCCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAggt                                                                            <  1:262823/244‑1 (MQ=255)                                                                                                                                                                                                   ggCAATACCGCCGCTATCTGTATCTATCCTAGATTTCTCCCGTTTGCTCGCAAAACACTGCACGAGCAGGCCCCCCAGGAAATCCGTATTCCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGAGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGCCGTTGACGTGGTGTTCCCGTACCGCGCGCTGCTGGCGGTTTACGCGCAggttcgcttttacctctgcgactc                                                       >  2:182439/1‑231 (MQ=255)                                                                                                                                                                                                                  atcggTATCTAACCTCGCTTTATCCCGCTTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAAACCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAAGCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGAGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTg                                       <  2:262054/247‑1 (MQ=255)                                                                                                                                                                                                                        atctatcCTCGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGATGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTGAAGCGCAGGTTGGGTTT‑ACACTGGTGAAAGCCttgtaagaggcgaggcgagg                                 >  2:64970/1‑232 (MQ=255)                                                                                                                                                                                                                                 cGCTTTATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCggaagt                        >  1:382374/1‑246 (MQ=255)                                                                                                                                                                                                                                     ttATCCCGATTGCTCGCAAAACACTGAAAGAGCAGGGCACCCCGGAAATCCGTATTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAATGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTAct                     <  1:338816/250‑1 (MQ=255)                                                                                                                                                                                                                                                gCTCGCAAAACACTGAAAGAGCCGGGCACCCCGGAAATCCGTATTGCGACGGTAACCAACTGCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGa          <  2:245686/250‑1 (MQ=255)                                                                                                                                                                                                                                                       aaaCACTGAAAGAGCAGGGCACCCCGGAAATCCGTAGTGCTACGGTAACCAACTTCCCACACGGTAACGATGACATCGAAATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCCGATGAAGTTGACGTGGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGaa  <  1:182439/251‑1 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                        GCGGAGTAGATGTTAGAATACTAACAAACTCGCAAGGTGAA‑TTTTATTGGCGACAAGCCAGGAGAATGAAATGACTGATCTGAAAGCAAGCAGCCTGCGTGCACTGAAATTGATGGACCTGACCA‑CCCTG‑AATGACGACGACACCGACGAGAAAGTGATCGCCCTGTGTCATCAGGCCAAAACTCCGG‑TCGGCAATACCGCCGCTATCTGTATCTATCCTCGCTTTATCCCGATTGCTCGCAAAACTCTGAAAGAGCAGGGCACCCCGGAAATCCGTATCGCTACGGTAACCAACTTCCCACACGGTAACGACGACATCGACATCGCGCTGGCAGAAACCCGTGCGGCAATCGCCTACGGTGCTGATGAAGTTGACGTTGTGTTCCCGTACCGCGCGCTGATGGCGGGTAACGAGCAGGTTGGTTTTGAC‑CTGGTGAAAGCCTGTAAAGAGGCTTGCGCGGCAGCGAATGTACTGCTGAAAGTGATCATCGAA  >  NC_000913/4617253‑4617745

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