Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,923,941	C→T	R58Q (CGA→CAA)	yqcC ←	DUF446 domain‑containing protein YqcC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,923,941	0	C	T	90.9%	53.3 / ‑1.6	22	R58Q (CGA→CAA)	yqcC	DUF446 domain‑containing protein YqcC
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (4/16);  total (5/17)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-01

CGCGATCCAGCCAGCTGCGGTGAACCAGCCAGCCGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATAAAAAAAGGGAGCACTGTATT  >  NC_000913/2923697‑2924179                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               cgcgATCCAGCCAGCTGCGGTGAACCAGCCAGCCGGAGGGTTTATTTACCGCCACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTCGCTCGCCTCATCCGCAAAAAGAGCATCCCGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCTTAATAGGGTGCACCGGCAGAAGCGCCCGGTCACGGCTGTTTGTTATCTCACAAATCGTGCATTTGCggtg                                                                                                                                                                                                                                          >  2:87581/1‑249 (MQ=255)                                            atttaCCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCAt                                                                                                                                                                                               <  1:13933/251‑1 (MQ=255)                                                                cATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTa                                                                                                                                                                            <  1:84834/250‑1 (MQ=255)                                                                      gATAGAGTATTTCCAGCATTAGCGCGCATCATCCTCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCATAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAAtt                                                                                                                                                                     <  2:125508/251‑1 (MQ=255)                                                                      gATAGAGTAGTTCCAGCATTAGCGCGCATCATCCGCAAAAAGAGCAGCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTTCGGTGCCAGCGCCATTGCATAATAGGGGGCAACGGCAAAAGCGCCCGGGAACTGCTGTGTGTTAGCTAACAAAGCGTGCATTTGCGGGATCATCACCCACGGCAGCGTCTCCAGTGGGTCCAGGGTTGCCATTAAGAACGGTTGAGGACTATTaaa                                                                                                                                                                       <  2:18688/249‑1 (MQ=255)                                                                       aTAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTg                                                                                                                                                                    <  1:62396/251‑1 (MQ=255)                                                                                              cGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTc                                                                                                                                             <  2:188720/251‑1 (MQ=255)                                                                                                   catcCGCAAAAAGAGCATCCAGTGTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTGTCGCCa                                                                                                                                         <  1:207584/250‑1 (MQ=255)                                                                                                                              ccAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCATCCACTCCAGTGGTTCCATGGTGTCCATAAAGATCGGTTGGGTACTATTAAATTGATGCGTCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCa                                                                                                             <  2:164110/251‑1 (MQ=255)                                                                                                                                                              tGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATg                                                                             >  2:72023/1‑251 (MQ=255)                                                                                                                                                                    gTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGt                                                                       >  2:199887/1‑251 (MQ=255)                                                                                                                                                                      cGCCAGCGCCATTGCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCa                                                                     <  2:213494/251‑1 (MQ=255)                                                                                                                                                                                       aataGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGATCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTGCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATaa                                                    <  1:36655/251‑1 (MQ=255)                                                                                                                                                                                             gTGCAACGGCAAAAGCGACCGGTAACAGCTCTTTGGTATCTAACAAATCGTGCATGTGCGGGATCAACACGCTCTGCAACCACTGCAGTGGTTCCATGGTGTCCATAAAGAACAGGTGGGTACTATTTAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTTAGGTGCACGCAGTAACGCTTCAAGCGCCTGCAAGTGGAGACGAACGCGGTCAGGAGTGGTCATAAGTCAGCCTCAATAAGaaaaa                                              <  2:31522/251‑1 (MQ=255)                                                                                                                                                                                               gCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCCCTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAAc                                             <  2:55507/250‑1 (MQ=255)                                                                                                                                                                                                    ggCAAAAGCGCCCGGTAACGGCGGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGGTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAACCTGCcgc                                       <  2:146392/251‑1 (MQ=255)                                                                                                                                                                                                             gcCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAg                              <  2:198639/251‑1 (MQ=255)                                                                                                                                                                                                               ccggtcaaggCTGGTTGTTCTCTAACACATCGTTCATTTGCGGGATCAGCACCCACTGCAACCAATCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGCTGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGCCCTAAAGCAACCTCAATAAGAAACACTGCCGCGAAGGGTAGc                             >  2:131924/9‑250 (MQ=255)                                                                                                                                                                                                                    aaCGGCTGTTTGTTATCTAACAAATCGTGCATGTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATGGATGCGGCTGGGGTTCAGCGTTTCGCCAGTGCTGATGGTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCAGAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAAt                       <  1:126578/251‑1 (MQ=255)                                                                                                                                                                                                                     aCGGCTGGTTGTGAGCTAACAAATCGTGCATTTGCGGGATCAGCACCCACTGCAACAACTCCAGGGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCGGGGGTGCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATa                      <  1:87581/251‑1 (MQ=255)                                                                                                                                                                                                                               gttATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATAAAAAAAGGGa            <  1:72023/251‑1 (MQ=255)                                                                                                                                                                                                                                         aaaTCGTGCATTTGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACTCAGGAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATACAAAAAGAAAGCACTGTAtt  >  1:72635/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                               CGCGATCCAGCCAGCTGCGGTGAACCAGCCAGCCGGAGGGTTTATTTACCGCAACCAGCCATTCATCCTGATAGAGTATTTCCAGCATTAGCTCGCATCATCCGCAAAAAGAGCATCCAGTTTTTCCAGCTCAGCCAGAATAAGCGCGCGTTGCGGATGGTCCGTCGCCAGCGCCATTTCATAATAGGGTGCAACGGCAAAAGCGCCCGGTAACGGCTGTTTGTTATCTAACAAATCGTGCATTCGCGGGATCAGCACCCACTGCAACCACTCCAGTGGTTCCATGGTGTCCATAAAGAACGGTTGGGTACTATTAAATTGATGCGGCTGGGGTTCATCGTTTCGCCAGTGCTGATGTTCACGCAGTAACGCTTCAAGCGCCTGCAACTGGAGACGAACGCGGTCATGAGTGGTCATAAAGCAACCTCAATAAGAAAAACTGCCGCGAAGGATAGCAGGAATAAAAAAAGGGAGCACTGTATT  >  NC_000913/2923697‑2924179

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