Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,429,941	C→T	G28D (GGT→GAT)	yrdB ←	DUF1488 domain‑containing protein YrdB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,429,941	0	C	T	94.1%	37.2 / NA	17	G28D (GGT→GAT)	yrdB	DUF1488 domain‑containing protein YrdB
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (4/12);  minor base G (0/1);  total (4/13)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGCTATTCCGCGAATAATTACGTGAAATGGAAGGACGAATATCTGGATCAGGGTAACCAGACCCAGCCTTGATCATCTTCACTTTGTTCCTGAATTAAGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATACCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCGGAAACTGGATGGCCTGATTCACGCGGACAATTCCTCCTGCAATTGCTTTATAACTGGTTCTACGTCAGGCAGAACACCGTGCCAGAGAAGAAAGGCATGAGCCGCCTGTGCCACCAGCATTCCTAAACCATCAGCATTACGCTTTGAGCCTCGCTG  >  NC_000913/3429719‑3430157                                                                                                                                                                                                                               |                                                                                                                                                                                                                         cGCTATTCCGCGAATAATTACGTGAAATGGAAGGACGAATATCTGGATCAGGGTAACCAGACCCAGCCTTGATCATCTTCACTTTGTTCCTGAATTAAGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATAGCATTCACGAGAGCgg                                                                                                                                                                                                           <  1:260384/238‑1 (MQ=255)        ccGCGAATAATTACGTGAAATGGAAGGACGAATATCTGGATCAGGGTAACCAGACCCAGCCTTGATCATCTTCACTTTGTTCCTGAATTAAGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACAtttttt                                                                                                                                                                                        >  1:178106/1‑250 (MQ=255)                          naTGGAAGGACGAATATCTGGATCAGGGTAACCAGACCCAGCCTTGATCATCTTCACTTTGTTCCTGAATTAAGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCTTcc                                                                                                                                                                     >  1:194132/2‑251 (MQ=255)                                 gtaCGAATATCTGGATCAGGGTAACCAGACCCAGCCTTGATCATCTTCACTTTGTTCCTGAATTAAGTTTTCCGCTTATTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGCCTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGCAAACATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCg                                                                                                                                                              >  1:310813/3‑251 (MQ=255)                                              gATCAGGGTAACCAGACCCAGCCTTGATCATCTTCACTTTGTTCCTGAATTAAGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCGGGCGTATCTCCAGTAAAGCGAGACGCCAGACTCTCGCCAGAGAGCGCGCATGTCAGTTTGATATCATTCACGAGAGCGGGAAAGCAGACACATTTTTGATTCTCGTCCCACTCTTCCCTGTCCGGAAACTGGATggn                                                                                                                                                 <  1:60102/251‑2 (MQ=255)                                                                           atcTTCACTTTGTTCCTGAATTAAGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCGGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATac                                                                                                                                                                                                 <  1:14946/174‑1 (MQ=255)                                                                               tCACTTTGTTCCTGAATTAAGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCGGAAACTGGATGGCCTGATTCACGCGGACAATTCCTCCTGCAATTg                                                                                                                 <  1:74755/250‑1 (MQ=255)                                                                                                tAAGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTGTATTCGCGTCCCACTCTTCCCTGTCCGGAAACTGGATGGCCTGATTCACGCGGACAATTCCTCCTGCAATTGCTTTATAACTGGTTc                                                                                                  <  1:238933/248‑1 (MQ=255)                                                                                                 aaGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCGTCCCTGTCCGGAAACTGGATGGCCTGATTCACGCGGACAATGCCTGCTGCAATTTCTGTATAACTGTTTCTACg                                                                                              <  1:112281/251‑1 (MQ=255)                                                                                                       tCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTTTATtct                                                                                                                                                                                   <  1:118087/160‑1 (MQ=255)                                                                                                                     cAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGACAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCGGAAACTGGATGGCCTGATTCACGCGGACAATTCCTCCTGCAATTGCTTTATAACTGGTTCTACGTCAGGCAGAACACCGTGcc                                                                           >  1:280857/1‑250 (MQ=255)                                                                                                                                       aCGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGGAAAGCGATACGCCAGACTCTCGCCAGAGAGCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAGCATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCGGAAGCTGGATGGCCTGATTCACGTGGACAATTCCTCCTGCAATTGCTTTATAACTGGTTCTACGTCAGGCAGAACACCGTGCCAGAGAAGAAAGGCATGAgc                                                        <  1:233121/251‑1 (MQ=255)                                                                                                                                                   tAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCGGAAACTGGATGGCCTGATTCACGCGGACAATtcctc                                                                                                                          <  1:207735/173‑1 (MQ=255)                                                                                                                                                                                   nnnnnnnnnnaCNNNNNCCANAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCGGAAACTGGATGGCCTGATTCACGCGGACAATTCCTCCTGCAATTGCttt                                                                                                             <  2:194132/144‑1 (MQ=255)                                                                                                                                                                                       cGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCGGAAACGGGATGGCCTGATTCACGCGGACAATTCCTCCTGCAATTGCTTGATAACTGGTTCTACGTCAGGCAGAACACCGTGCCaga                                                                        <  1:188580/187‑1 (MQ=255)                                                                                                                                                                                                ctcGCCAGAGATCGTGCATGTCAGTTGCATATCATTCACGAGAGCGTGAAAACATACACATTTTTTATTCTCGGCCCACTCTTCCCTGTCCGGAAACTGGATGGCCCTATTCACGCGGATAATGCTTGCTGCAATGGCTTTATAACTGGTTGTACGTCTTGCAGAACACCGTGGCAGAGAAGAAAGGCAGGAGCCGCCTGTGCCACCAGCATTCCTAAACCATCAGCATTACGCTTTGAGCCTCgctg  <  1:160604/248‑1 (MQ=255)                                                                                                                                                                                                             gcgcAAGTCAGTTGCATATCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCGGAAACGGGATGGCCTGATTCACGCGGACAATTCCTCCTGCAATTGCTTTATAACTGGTTCTACGTCAGGCAGAac                                                                                   <  1:321934/154‑1 (MQ=255)                                                                                                                                                                                                                               |                                                                                                                                                                                                                         CGCTATTCCGCGAATAATTACGTGAAATGGAAGGACGAATATCTGGATCAGGGTAACCAGACCCAGCCTTGATCATCTTCACTTTGTTCCTGAATTAAGTTTTCCGCTTCTTCTTCCAGGTCCCAGCGATGCTGACGAAAACTCGCTAACCACTGTTCTGGCGTATCTCCAGTAAAGCGATACGCCAGACTCTCGCCAGAGATCGCGCATGTCAGTTGCATACCATTCACGAGAGCGGGAAAACATACACATTTTTTATTCTCGTCCCACTCTTCCCTGTCCGGAAACTGGATGGCCTGATTCACGCGGACAATTCCTCCTGCAATTGCTTTATAACTGGTTCTACGTCAGGCAGAACACCGTGCCAGAGAAGAAAGGCATGAGCCGCCTGTGCCACCAGCATTCCTAAACCATCAGCATTACGCTTTGAGCCTCGCTG  >  NC_000913/3429719‑3430157

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