Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,017,470	C→T	P471S (CCG→TCG)	ygfK →	putative oxidoreductase, Fe‑S subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,017,470	0	C	T	96.3%	89.2 / ‑5.3	27	P471S (CCG→TCG)	ygfK	putative oxidoreductase, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (18/8);  minor base G (0/1);  total (18/9)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

ACCTGGCGGCTATCTGCGCTTAAGTGCCTGCATGCGCGAGCTGGAAGGCTCCGACGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTCCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGCCCGCCATTACCGGTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATC  >  NC_000913/3017232‑3017713                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    aCCTGGCGGCTATCTGCGCTTAAGTGCATGCATGCGCGAGCTGGAAGGCTCCGACGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATcc                                                                                                                                                                                                                                         >  2:149631/1‑251 (MQ=255)     ggcggcTATCTGCGCTTAAGTGCCTGCATGCGCGAGCTGGAAGGCTCCGACGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCt                                                                                                                                                                                                                                     <  1:299760/251‑1 (MQ=255)        ggcTATCTGCGCTTAAGTGCCTGCATGCGCGAGCTGGAAGGCTCCGACGCCTGGGGCCTTGACCATGTTGACGTCGAACGACTGAACAGCCTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCt                                                                                                                                                                                                                                  >  1:157874/1‑251 (MQ=255)                  gcTTAAGTGCCTGCATGCGCGAGCTGGAAGGCTCCGACGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAacac                                                                                                                                                                                                                        >  2:242308/1‑251 (MQ=255)                     tAAGTGCCTGCATGCGCGAGCTGGAAGGCTCCGACGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTAACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGc                                                                                                                                                                                                                     >  1:253199/1‑251 (MQ=255)                     tAAGTGCCTGCATGCGCGAGCTGGAAGGCTCCGACGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTGGCCCCCTGTGTTACTGCCGGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACAcc                                                                                                                                                                                                                       <  1:331164/249‑1 (MQ=255)                       agggctggcaggcgcGAGCGGGAAGGCTCCGACGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGAGGCGTTACCCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTTCGCGATCAAGCAAGATATGTCGGAATACATCCGTCTGCTTGGCGAACCCCGCTa                                                                                                                                                                                                                   <  2:34903/240‑1 (MQ=255)                                     ggAGCTGGAAGGCACCGACGCCTGGGGACGTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCAGGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTgg                                                                                                                                                                                                     <  2:124084/250‑1 (MQ=255)                                         cTGGAAGGCTCCGACGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGCCGCGCTGGAACt                                                                                                                                                                                                 >  2:266829/1‑251 (MQ=255)                                          ggggaGGCTCCGATGCCTGGGGACTTGACCATGTGGACGACGAACTACTGAACAGACTGGCAGCAGATGCGTTAACCAGGGAATACCCCCAGAAACCCTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAATACCTGCCGCTGACCGACTGCTACGTTGCCCACTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGACGACGCGCTGGAACTc                                                                                                                                                                                                <  2:273064/247‑1 (MQ=255)                                                      atgCCTGGGAACTTGACCATGGTGACGCCAAACGACTGAACAATCTGCCAGCAGAGGCGTTAACCAGGGAACACACCCAGAAACAGGGGAAGCCCGAAGAGCAGATGGAAGAGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTCACTGCCTGCGCGATCAAGCCAGAAATTGCGGAATACAGCCGTCTGCTTGGCGAACAACGCTATGCCGACGCGCTGGAACTCATCGACcaccgc                                                                                                                                                                                    <  2:242600/249‑1 (MQ=255)                                                      aCGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACcaacgc                                                                                                                                                                                    >  1:80294/1‑251 (MQ=255)                                                         ccTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACGGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATGTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACcaacgcaac                                                                                                                                                                                 <  2:16149/251‑1 (MQ=255)                                                             gggACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGctc                                                                                                                                                                             <  1:294773/251‑1 (MQ=255)                                                                aCTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGc                                                                                                                                                                          >  1:192379/1‑251 (MQ=255)                                                                 cTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGcc                                                                                                                                                                         >  2:238789/1‑251 (MQ=255)                                                                   tGACCATGTGGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCAGGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGCCCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATGTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCGACCAACGCAACTCTCTGccc                                                                                                                                                                        <  2:289865/250‑1 (MQ=255)                                                                                        aCTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGCCCGCCATTACCGGTCATATTTGCg                                                                                                                                                  >  1:294977/1‑251 (MQ=255)                                                                                                  cTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGCATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGCCCGCCATTACCGGTCATATTTGCGATCACCAGTg                                                                                                                                        >  2:276449/1‑251 (MQ=255)                                                                                                               cGTTAACCATGGAATACACCCAGACACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGCCCGCCATTACCGGTCATATTTGCGATCACCAGTGCCAATACAACTGt                                                                                                                           >  2:177224/1‑251 (MQ=255)                                                                                                                                     gAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGCTCTTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGCACTCATCTACCAACGAACCGCTCTGCCCGCCATTACCGGTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGACAGGATTAAGACAGTg                                                                                                     >  2:123082/1‑251 (MQ=255)                                                                                                                                           cTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGCCCGCCATTACCGGTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTg                                                                                                >  1:168998/1‑250 (MQ=255)                                                                                                                                                                  aaGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCCAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGCCCGCCATTACCGGTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAgac                                                                        >  1:350544/1‑249 (MQ=255)                                                                                                                                                                  aaGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGCCCGCCATTACCGGTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTc                                                                        >  1:217620/1‑251 (MQ=255)                                                                                                                                                                                            gACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGACCTCATCTACCAAAGCAACGCTCTGCCCGACATTACCGGTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCAACCTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAat                                               >  1:287902/1‑250 (MQ=255)                                                                                                                                                                                                     gttgccgcatgaggtACTGCCTGCGCTATCAAGCAAGATATTTCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGACGACGCGCTGGAACTCATCTCCCAACGCAACGCTATGCACACCATTACCGGTCATATTTGAGATCACCAGTGCCAATACAACCGTACCCGCCTGGATTCCGACAGTGCGCTGAATATCCGCGAACTGAAACAAGTCGCGCTGGAAAAAGGATGGGATGAATATAAGCCACg                                     >  2:321813/15‑251 (MQ=255)                                                                                                                                                                                                                                        agATATTTCTGAATACATCCGTCTGCTTGGCGGACACCGCTATGCCGAAGCGCTGGAACTCAGCTACCAACACAACGCTCTGCCCTCCATTACCGGTCATATTTGCGATCACCAGTGCCAATACAACTGTCCCCGCCTGGATTCCGACAGTGCGCTGAATATCCGCGAACTGGGAAAAGACGCGCTGGAAGAAGTATGGGGTGAATATAACCAACGCTGGAACAAAACAGCCGGTTATGGTTCACGCCATc  >  2:253932/1‑251 (MQ=255)                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                    ACCTGGCGGCTATCTGCGCTTAAGTGCCTGCATGCGCGAGCTGGAAGGCTCCGACGCCTGGGGACTTGACCATGTTGACGTCGAACGACTGAACAGACTGGCAGCAGATGCGTTAACCATGGAATACACCCAGAAACACTGGAAGCCAGAAGAGCGTATTGAAGTGGCAGAAGACCTGCCGCTGACCGACTGCTACGTTGCCCCCTGTGTTACTGCCTGCGCTATCAAGCAAGATATTCCGGAATACATCCGTCTGCTTGGCGAACACCGCTATGCCGACGCGCTGGAACTCATCTACCAACGCAACGCTCTGCCCGCCATTACCGGTCATATTTGCGATCACCAGTGCCAATACAACTGTACCCGCCTGGATTACGACAGTGCGCTGAATATCCGCGAACTGAAAAAAGTCGCGCTGGAAAAAGGTTGGGATGAATATAAGCAACGCTGGCACAAACCAGCCGGTTCTGGTTCACGCCATC  >  NC_000913/3017232‑3017713

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