Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,715,661	C→T	D77N (GAT→AAT)	pdxY ←	pyridoxal kinase 2

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,715,661	0	C	T	93.8%	48.3 / NA	16	D77N (GAT→AAT)	pdxY	pyridoxal kinase 2
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (6/9);  minor base G (0/1);  total (6/10)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCT  >  NC_000913/1715418‑1715896                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            ccGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGgtg                                                                                                                                                                                                                                      >  1:57668/1‑251 (MQ=255)   gCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAAGGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGAAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGgtgtg                                                                                                                                                                                                                                    <  1:167219/251‑1 (MQ=255)       gCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTaa                                                                                                                                                                                                                                 <  2:186812/250‑1 (MQ=255)                             acccAGATTTGGCGCAATGATATCGCTGGCAGTCAAACCGTGTCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTGTTTCTGGATGACCCATTCCCGGATCGCAAAAAGATTTCGCCTGGTGATTCGCGGCTTTCACCGTACGGACTATACCGAGGATATTTTCTCCTTGCTCCGCCGATCCCAGATTGCCACTTAATACGGCATTACAGGTGTGTAATTGATCAATGG‑GGCAATGCCTTg                                                                                                                                                                                                         <  1:138806/249‑1 (MQ=255)                                 aGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGAGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATAACGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTCATTTATCAATGGCGGCAAAGCCTTGcac                                                                                                                                                                                                      >  2:80310/1‑251 (MQ=255)                                                       tGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggc                                                                                                                                                                                >  2:257074/1‑251 (MQ=255)                                                        ggcaggcaAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGCCGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGggcg                                                                                                                                                                               >  1:24037/1‑251 (MQ=255)                                                                                       gCGACCCCCGGTGAAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGTATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAAGGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCATCCAGTCCATTTGCCGTATTg                                                                                                                                                <  2:189889/251‑1 (MQ=255)                                                                                       gCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTAtt                                                                                                                                                 <  1:257074/250‑1 (MQ=255)                                                                                                      aCGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGaaa                                                                                                                                  >  2:45714/1‑250 (MQ=255)                                                                                                        gATACAGCCTTTTCCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCGGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCGGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCGTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAAtt                                                                                                                               <  2:154300/251‑1 (MQ=255)                                                                                                             aGCCTTTTTCCGGATGACCCATTACCGGACCGCAAAAATATTGCGCCTGCGGATTCGCGGCTTTCCCCTGACGGACGATACCGAGGATATGTTCCCCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACg                                                                                                                          <  1:58618/251‑1 (MQ=255)                                                                                                                 tttttCCGGATGACCCATTACCGGATCGCAAAAATATTACGCCTGAGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATGGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGgtgt                                                                                                                      >  1:175257/1‑251 (MQ=255)                                                                                                                                               aaaaTATTGCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATGACATGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTTGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACTTTCGCGCCCAggcggc                                                                                        <  2:24037/251‑1 (MQ=255)                                                                                                                                                                                                     attttCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGAc                                   <  2:66405/247‑1 (MQ=255)                                                                                                                                                                                                                                     cACTTAATACGGCATTACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCt  <  1:259263/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                            CCGCATGCTCACAGAGTATTTCCAGCTCAACCAGATTTGGCGCAATGATATCGCTGGCAGGCAAACCGTGCCGCACATGAAACTCTGCGACACCCGGTGCAACGATACAGCCTTTTTCCGGATGACCCATTACCGGATCGCAAAAATATTTCGCCTGCGGATTCGCGGCTTTCACCTGACGGACGATACCGAGGATATGTTCACCCTGCTCCGCCGATCCCAGATAGCCACTTAATACGGCATCACAGGTGTGTAATTTATCAATGGCGGCAATGCCTTGCACAATTTCGGTTAAATGGCTGGGCGGCATCACGCAGCCAGTCCATTTGCCGTATTGGGTGTGATTAGAAAATTGAACGGTGTTCAGCGGCCAGACGTTCGCGCCCAGGCGGCGCATCGGAAACTCTGCCGCACTGTTACCCGCATGACCATAAACAACGTGAGACTGGATAGCGAGAATATTTTTCATCATGTTCCCT  >  NC_000913/1715418‑1715896

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