Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	292,377	C→T	S191N (AGC→AAC)	yagK ←	uncharacterized protein YagK

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	292,377	0	C	T	94.1%	41.1 / NA	17	S191N (AGC→AAC)	yagK	uncharacterized protein YagK
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (9/7);  minor base G (0/1);  total (9/8)
Fisher's exact test for biased strand distribution p-value = 4.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATG  >  NC_000913/292130‑292619                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   cGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGtttt                                                                                                                                                                                                                                                 >  1:310712/1‑251 (MQ=255)      gtTATATTTGAAACGGCAGGATGTGCAGGATGCTGTTTCTACCGCTTAAAATGCAGAAACGTTATTTAGTAGCGTACTGGCTGATGACTATCAGTAATCTATAGAAAATTGCGAAACGCAATGGCAGGTTAGTAATTTTTTGAAACAATACAAATAATGCCTTCTCACTGAGAATGCCAGTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACTCTTTGTTTGCAAc                                                                                                                                                                                                                                             <  2:409896/250‑1 (MQ=255)        tAGAATTGAAACAGCAGGAGTTGCAGGATGCTGTTGCTACCGCTGAAAATGCATAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATATTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGGTTTCAACtt                                                                                                                                                                                                                                           <  2:292785/250‑1 (MQ=255)               aaaacggcaggtatttcattaTGCGGTTGCTACCGCTGAAAAATCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACGATCAGTAATCTATAGAAAATTGCGTAACTAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTGCTCGCTGAGAAGGCCATTACCGACTTACAGTTCGTTTTGGCTACAGCCGAGGTTACGGTCTCCTTCGCCAAACACTTTGTTTTCAACTTTGGTCAgg                                                                                                                                                                                                                                   <  2:3326/231‑1 (MQ=255)                                                                                                          ataGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATgg                                                                                                                                         <  1:32049/250‑1 (MQ=255)                                                                                                            aGAAAATTACGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCATTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGGCAGTATTATTGGTGTCTAGCACATACTTAAAGTTTTCAGGAAAATGgacc                                                                                                                                      >  2:400815/1‑251 (MQ=255)                                                                                                                                        aGTAAATTTTTGAAACAATACAAATATGGCCCTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCAGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTTCAgc                                                                                                           <  2:289321/250‑2 (MQ=255)                                                                                                                                                                      cTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGatc                                                                             >  2:174435/1‑250 (MQ=255)                                                                                                                                                                           cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATcc                                                                         <  1:174435/249‑1 (MQ=255)                                                                                                                                                                           cACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCct                                                                       >  1:141102/1‑251 (MQ=255)                                                                                                                                                                            aCTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCCCATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATcctccctc                                                                      >  2:393888/1‑251 (MQ=255)                                                                                                                                                                                     gCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTCAGGTTAt                                                             >  1:143408/1‑251 (MQ=255)                                                                                                                                                                                               aCTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATc                                                   >  1:265839/1‑251 (MQ=255)                                                                                                                                                                                                 ttccAGTTCGCTTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGt                                                 >  2:305205/4‑251 (MQ=255)                                                                                                                                                                                                                  tACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGcaca                                                                                                                                                                 >  2:149013/1‑122 (MQ=255)                                                                                                                                                                                                                  tACAGCCGAAGTTACGGTCTCCTTCACCAAACACGTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGcaca                                                                                                                                                                 <  1:149013/122‑1 (MQ=255)                                                                                                                                                                                                                                         tCACCAAACACTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGGGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTTATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGatcaac           <  2:143408/249‑1 (MQ=255)                                                                                                                                                                                                                                                   cTTTGTTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAAAAGCGCCTGGTAGTTCAGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGATTTCAGGAAAATGGAACAGACCAGGGTGATCGTCTCTATCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAAtg  >  1:184870/1‑248 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   CGAGTGTTAGATTTGAAACGGCAGGATGTGCAGGATGCTGTTGCTACCGCTGAAAATGCAGAAACGTTATTTAGTAGCGGACTTGCTGATGACTATCAGTAATCTATAGAAAATTGCGTAACGAAACTGCATTTTAGTAATTTTTTGAAACAATACAAATAAGGCCTTCTCACTGAGAAGGCCATTACCGACTTACAGTTCGATTTGGCTACAGCCGAAGTTACGGTCTCCTTCACCAAACACTTTGCTTTCAACTTTGGTCAGGTAGTCCAGTCGGGTTAACAGCGCCTGGTAGTTCTGCTGGAAGTCAGTACTATTGGTGTCTAGCACATACTTACAGTTTTCAGGAAAATGGACCAGACCAGGGTGATCGTCTCTCTCCAGTCTCAGGGCGCTGTACCAGGCTCCTGTGATCATCCCTCTAAGGTTATCGTCCTGATCGTAATCCCCGAGGTGGTAGTAAGCGTCTTTGTTGAACAACAGGCAAATG  >  NC_000913/292130‑292619

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