Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,931,805	C→T	R251W (CGG→TGG)	ygdG →	flap endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,931,805	0	C	T	94.7%	45.3 / NA	19	R251W (CGG→TGG)	ygdG	flap endonuclease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (3/15);  minor base G (0/1);  total (3/16)
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

GGGGACTTGCGGGGATCAGCAGTTCAAAGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTTGCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAATGTGGCGGAATACCGACATCACGGTTGAGAGCAAACA  >  NC_000913/2931560‑2932037                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         ggggACTTGCGGGGATCAGCAGTTCAAAGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTTGACGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTtggtatggta                                                                                                                                                                                                                                      >  1:255839/1‑250 (MQ=255)                                aCCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTTGCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGTCGGGAACGCCGCacacaccgtat                                                                                                                                                                                                      >  2:194200/1‑241 (MQ=255)                                                 tCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTTGCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGAAACGCCCCCAACGACGTATCCGGcctttttacatcc                                                                                                                                                                                     >  1:153090/1‑239 (MQ=255)                                                          aaaGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTTGCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTa                                                                                                                                                                            <  2:11785/251‑1 (MQ=255)                                                                 cACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTTGCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAACCGGGTAAGCCGCCAACGCCGCATCCGGCATTATctcacccacgactttaccacgcg                                                                                                                                                                     >  2:83032/1‑229 (MQ=255)                                                                         ggctgggtgAGTTTCAGAGTCTGGAAGGGAAATATGAGAATCTGGATGCGGTTGCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTGTCTGTGCCGCGATATTGCCCGCTTACAAACCGAGTTGCATAGCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCTTAGCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTa                                                                                                                                                             <  2:268209/243‑1 (MQ=255)                                                                         ggctgGTCGAGTTTCAGAGTATGGAAGGGATATATGAGAATCTGGATGCGGTTGCCGAAAAGTGGCGAAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATAGCGACGGCAATTTACAGAAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCAGTTACCAGGCGGTATGGTa                                                                                                                                                             <  2:239468/250‑1 (MQ=255)                                                                                      tCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTTGCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAat                                                                                                                                                 <  1:229467/250‑1 (MQ=255)                                                                                                   gggATATATGAGAATCTGGATGCTGTTGCCGAAAAGTGGCGGAAAAAATTAGTAACCCATAAAGAGAGGGCGTTTCTGTTCCTGGATAGTGACCGCTTACAATCCGATTTGCATATCGACTGCAATTTACAGAAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGCCGGTATGGTAAAGCTCTACAATATCCTCAAGCGtt                                                                                                                                    <  2:67948/250‑1 (MQ=255)                                                                                                                    ggATGCGGTTGCAGAAAAGTGGCGCAAAAAAGTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCGCATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAAGATCCTCAAGCGTTCCTTCACGCGGGTTGcc                                                                                                                   <  2:230428/250‑1 (MQ=255)                                                                                                                                            cAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTATAAACATCATCCAGTgcc                                                                                          <  2:224070/251‑1 (MQ=255)                                                                                                                                            cAAAAAATTAGAAACCCAAAAAGAGAGGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATTTCGACGGCAATTTACAGCAATTGCGGTTGGTAGGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgcc                                                                                          <  1:52293/251‑1 (MQ=255)                                                                                                                                                                     aTGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCtt                                                                 <  1:83032/251‑1 (MQ=255)                                                                                                                                                                                          atTGCCCGCTTACAAACCGAGTTGCATATCGACGGCAATTTACAGCAAGTGCGGTTGGTATGGTAACTGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGTCGGTATGGTAAAGCTCTACAATATCCTCAAGCTGTGCTTCACGCGGGTGGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATctc                                            <  2:37203/251‑3 (MQ=255)                                                                                                                                                                                                    tACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCGTGCGTACACCAACATCACGCAAAtgtg                                   <  2:136370/250‑1 (MQ=255)                                                                                                                                                                                                     aCAAACCGATTTGCATATCGACGGCAAGTTACAGCAATGGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGCTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTCCAAACATCCTCCAGTGCCGCCTGCGCGAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAATGTGGc                                 <  1:229835/251‑1 (MQ=255)                                                                                                                                                                                                                     atCGACGGCAATTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAATGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAATGTGGCGGAATACCGACATCa                  <  1:194200/250‑1 (MQ=255)                                                                                                                                                                                                                       cgccGGCAATTTCCAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGAAGACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAAGGTGGCGTAATACCGACATCACgg               <  2:22575/248‑1 (MQ=255)                                                                                                                                                                                                                               aTTTACAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAATGTGGCGGAATACCGACATCACGGTTGAGAGc       <  1:19214/251‑1 (MQ=255)                                                                                                                                                                                                                                    cAGCAATTGCGGTTGGTATGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAATGTGGCGGAATACCGACATCACGGTTGAGAGCAAaca  <  2:153090/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                         GGGGACTTGCGGGGATCAGCAGTTCAAAGGTACCGGGTGTTGCGGGAATCGGACCAAAAAGCGCCACGCAGCTGCTGGTCGAGTTTCAGAGTCTGGAAGGGATATATGAGAATCTGGATGCGGTTGCCGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCGTACACCAACATCACGCAAATGTGGCGGAATACCGACATCACGGTTGAGAGCAAACA  >  NC_000913/2931560‑2932037

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