Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,643,554	C→T	E251K (GAA→AAA)	ydfU ←	protein YdfU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,643,554	0	C	T	96.5%	83.2 / NA	29	E251K (GAA→AAA)	ydfU	protein YdfU
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (17/11);  minor base G (0/1);  total (17/12)
Fisher's exact test for biased strand distribution p-value = 4.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTATGCAGCTCGTTGTGATGCTTTCTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGC  >  NC_000913/1643347‑1643790                                                                                                                                                                                                                |                                                                                                                                                                                                                                             ccGTATGCAGCTCGTTGTGATGCTTTCTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTc                                                                                                                                                                                                    <  1:37139/250‑1 (MQ=255)    tATGCAGCTCGTTGTGATGCTTTCTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCATGCCCGTGACCAATCAGATGAAGCGGATCGTCGGCTGGCATACCGCCGCAAGCACACGGCTGTGTCTTAACCCATCGTGTGTATTTCTCCTTAACCCAGAGGAGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTcttgg                                                                                                                                                                                                >  1:503021/1‑247 (MQ=255)                   gATGCTTTCTGCACAAAGGCAATACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTCCCCGTGACCCATCAGATGATGCGGATCGTCGGCAGGCATCACGCAGCCAGCACACGCCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCCACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGCTCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGc                                                                                                                                                                                 >  1:464073/1‑251 (MQ=255)                             gCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGt                                                                                                                                                                       <  1:160177/251‑1 (MQ=255)                                                                                nccGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGACTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACACTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCgg                                                                                                                    >  1:94426/2‑251 (MQ=255)                                                                                     gACCAATCAGATGATGCGGATCGTCGGCTTGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATGTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTac                                                                                                                                                                                                                <  1:255494/154‑1 (MQ=255)                                                                                                     cGGATCGACGGCTGGCATCCCGCAGAAAGCACACGGCTGTGTCTTAACACAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGGGCGTTGTTTCAATATGCAGGTGGCGGTCTGCTCTCCCGGTACGATGt                                                                                               >  1:139087/1‑251 (MQ=255)                                                                                                           gtcgGCTGGCATACCGCAGCAAGCACACGGCTGGGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGGCTTGTCCTTTGGTGATTTTTGTTGCTGGTGGGCGTCAGGCAACGGTGCAAGATTTTGTGTACGTTGTTGCAATATTCTGGTGGCGGTCTGCTCTCCCGGTACGATg                                                                                                <  1:192728/244‑1 (MQ=255)                                                                                                                             gcaagcaCACNGNTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTNTAGGCAGCTTNATGNNAGATNTNGNAGACTCTGGATCAACGGTGATGCTTACCANCGTCTTTTCCTGTGGTGATTTTTNNTGNNNGTNNGCGTa                                                                                                                                                                      <  2:369727/156‑1 (MQ=255)                                                                                                                             gcaagcaCACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTa                                                                                                                                                                      >  1:369727/1‑156 (MQ=255)                                                                                                                             gcaagcaCACGGCTGTGGCTTAACCCACCGTGTGTATTTATCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGCCTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCtctc                                                                                                           >  1:10132/1‑215 (MQ=255)                                                                                                                                 gnacacagatgtgTCTTAACCCAACGTGCGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGAtttttt                                                                                                                                                >  1:393270/10‑174 (MQ=255)                                                                                                                                 gcaCACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGCGACTCTGGATCCACGGTGATGCTTACCACCGTCTTTTACTGTGGTGATTTttgt                                                                                                                                                                                    >  1:334288/1‑138 (MQ=255)                                                                                                                                    cacGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTACCCCAGCGGAGACGTTTAGGCCGCTTCATGAAAGCTTTCGGCGACTCTGGATCAACGGTGCTGCTTACCACCGTCTTTTCCTGTGGTGATTTttgttg                                                                                                                                                                                  >  1:492624/1‑137 (MQ=255)                                                                                                                                        gCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTCACCCAGCGGCGACGATTAGGCAGCTTCCTGAACGATTTCGGAGACTATGGAAAAACGGTGATGCTTAACACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTACGGCAACGGTGCAAGATTTTTTGTGCGTTGttt                                                                                                                                     >  1:451167/1‑178 (MQ=255)                                                                                                                                           gtgtCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGCGGGCGTAAGGCAACGGTGCAAGATTGTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTg                                                           <  1:461226/249‑1 (MQ=255)                                                                                                                                             gtCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCCTGAAAGATTTCGGCCACTCCGGATCAGCGGCAATGCTTACCACCGTCATTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGATGGCGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGCGCAGATTTTTTCCGCTGGTa                                                        >  1:288290/1‑250 (MQ=255)                                                                                                                                                                                 nggcgACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTATGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTACGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCACAGCGAac                                            >  1:320588/2‑226 (MQ=255)                                                                                                                                                                                  ggcgANNTTNNNGANNCTTCATGAACGATTTCGGAGACTNTGGATCAACGGTGATGCTTAACAACGTCTNTNCNTGTGNNGCTNTTTGTTGCTGNTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTCGTTTCAAAATGNTNGTGGCGGTAAGCACTCCCGGTCCGATGTCGCNATCGCGGTACACCGGGCGGCTTTTNTCCTNNGGNannnnn                                                   >  2:379788/1‑213 (MQ=255)                                                                                                                                                                                  ggcgACGTTTAGGCAGCTTNATGAAAGATTNCGGNNACTCNGNANCAACGGTGATGCTTACCACCGTCTTTTCCNGTGGTGATTTTTGTTGCTGGTGGNNGTNNNGCNNCGGTg                                                                                                                                                           <  2:378234/114‑1 (MQ=255)                                                                                                                                                                                  ggcgACGTTTAGGCAGCTTCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTg                                                                                                                                                           >  1:378234/1‑114 (MQ=255)                                                                                                                                                                                  ggcgACGGTTAGGCAGCTGCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGTCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTGTCAATATTCTGGTGGCGGTCTTCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCg                                               <  1:379788/222‑1 (MQ=255)                                                                                                                                                                                     gACGTTTAGGCAGCTTCATGAAAGATTGCGTAGACTCTGGATCAACGTTGATGCTTACCACCGTCTTTTCCGGCGGTGATTGTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGc                                                                                                                      <  1:258164/148‑1 (MQ=255)                                                                                                                                                                                                    tCATGAAAGATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGc  <  1:246184/249‑1 (MQ=255)                                                                                                                                                                                                           ngATTTCGGAGACTCTGGATCAACGGCGATGCTTACCACAGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGt                                                                                   >  1:99507/2‑161 (MQ=255)                                                                                                                                                                                                            gATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGAtttt                                                                                                                                                                                       <  1:192856/60‑1 (MQ=255)                                                                                                                                                                                                            gATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTgn                                                                                                                                                          >  2:13204/1‑88 (MQ=255)                                                                                                                                                                                                            gATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGc                                                                                                                                                          <  1:13204/89‑1 (MQ=255)                                                                                                                                                                                                            gATTTCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGGGGCGGTCTGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAg                                                                           >  1:109844/1‑168 (MQ=255)                                                                                                                                                                                                            gATTTCGGAGACTCNNGATCAACGGNGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGNCAACGGGGCNNGANTTTTGNTNNNTNNNTTNAATATGCTGGTGGCGGTCTNCTNNCCCGGTNCGATGTCGCTTTCGCGGTACACCGAg                                                                           <  2:109844/168‑1 (MQ=255)                                                                                                                                                                                                                |                                                                                                                                                                                                                                             CCGTATGCAGCTCGTTGTGATGCTTTCTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTGCAAGATTTTTTGTGCGTTGTTTCAATATGCTGGTGGCGGTCAGCTCTCCCGGTACGATGTCGCTTTCGCGGTACACCGAGCAGATTTTTTCCGCTGGTAATCCCAGCGAACGACGCGATACAGCCTCAGGTAGTGCATCCACCACCTGATTGC  >  NC_000913/1643347‑1643790

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