Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,627,641	C→T	F109F (TTC→TTT)	nadR →	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,627,641	0	C	T	94.1%	45.3 / NA	17	F109F (TTC→TTT)	nadR	DNA‑binding transcriptional repressor/NMN adenylyltransferase NadR
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

AGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAT‑GGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATAT  >  NC_000913/4627401‑4627874                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           atggTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAgcg                                                                                                                                                                                                                                  <  1:207595/249‑3 (MQ=255)              cAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCgcagc                                                                                                                                                                                                                      >  1:235171/1‑250 (MQ=255)                                                cgcAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTAtta                                                                                                                                                                                   <  1:367895/251‑2 (MQ=255)                                                         tGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACtttt                                                                                                                                                                          >  1:357712/1‑251 (MQ=255)                                                               cGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATCTCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTCAata                                                                                                                                                                     >  2:167444/1‑250 (MQ=255)                                                                    tcccGTTTTTTGGGGCGTTAGTTTCCCCGGCAGAAGAAAACCATTTGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATAGCAGa                                                                                                                                                                <  2:357712/248‑1 (MQ=255)                                                                               ggggCTTGAGTTTCCCCGGCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCa                                                                                                                                                    <  1:149187/251‑1 (MQ=255)                                                                                                ggCAGAAGAAAACCATTGGTGTGGTATTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGa                                                                                                                                    >  1:163210/1‑250 (MQ=255)                                                                                                                         aTTCGGTAAGTTCTACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAT‑GGAGCCGTATCCGCACg                                                                                                         <  1:91710/251‑1 (MQ=255)                                                                                                                           tCGGTAAGTTCTAACCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTATAGCCAGGTTGACGAACTGCATATCATTGGGAGTTTTGACGATAACCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGACCCCCGTGCCGGATCGTCTGCGCTGGTTAGTACAGACTTATAAAAATCCGAAAAATATTCGCATTAATGCTTTCAACGAGGAGAG‑ATCGGAGCAGTATCCGCACGGc                                                                                                       >  2:315797/1‑251 (MQ=255)                                                                                                                                      tACCCGCTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAACTATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAT‑GGAGCCGTATCCGCACGGCTGGGAtgtgtg                                                                                            >  1:313221/1‑251 (MQ=255)                                                                                                                                                                               cTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCa                                                                                                                                        >  2:108786/1‑167 (MQ=255)                                                                                                                                                                               cTGTAGCCAGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCa                                                                                                                                        <  1:108786/167‑1 (MQ=255)                                                                                                                                                                               cTGTAGCCAGGATGACGTACTGAATATCATTATGGGTTTTGACGATACCCTCGATGGCGCGTTGCTTGAAGACAGCGCTATGTCGCAGCAGCCCCCCGTTCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCCGTACACTATTCGCATTCATGATTTCACCGACGAGGGCAT‑GGAGCCGTATCCGCACGGCTGGGATGCGTGGAGCAACGGCATCAAACAGTTTATGCCTggtagtgggattc                                                   >  2:231858/1‑239 (MQ=255)                                                                                                                                                                                       aGGTTGACGAACTGCATATCATTATGGGTTTTGACGATACCCGCGATCGCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAT‑GGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGgac                                           <  1:311302/251‑2 (MQ=255)                                                                                                                                                                                               gaacTGCATATCATTATGGGTTTTGACGATCCCCGCGACCCCGCGTTGTTTGAAGACAGCGCTATGTCGCAGCAGCCCCCCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAT‑GGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCAGgacctgatct                                    >  1:308236/4‑250 (MQ=255)                                                                                                                                                                                                                                cgcgatctcgcgTTGTTGGAAGACAGCGCTATGTCGCAGCAGCCCACCGTGCCGGATCGTCTGCGCTGGTTATTACAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAT‑GGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAGGGGATTCAGCCGGACCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTatat  <  1:74886/243‑1 (MQ=255)                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                           AGGGTATTTAAGCCAGTTACTGAATGCCAAAATCAAAAGCCCCAGCGCGCAAAAGCTGGAGGCGTTGCACCGTTTTTTGGGGCTTGAGTTTCCCCGGCAGAAGAAAACGATCGGTGTCGTATTCGGTAAGTTCTACCCACTGCATACCGGACATATCTACCTTATCCAGCGCGCCTGTAGCCAGGTTGACGAGCTGCATATCATTATGGGTTTTGACGATACCCGTGACCGCGCGTTGTTCGAAGACAGTGCCATGTCGCAGCAGCCGACCGTGCCGGATCGTCTGCGTTGGTTATTGCAAACTTTTAAATATCAGAAAAATATTCGCATTCATGCTTTCAACGAAGAGGGCAT‑GGAGCCGTATCCGCACGGCTGGGATGTGTGGAGCAACGGCATCAAAAAGTTTATGGCTGAAAAAGGGATCCAGCCGGATCTGATCTACACCTCGGAAGAAGCCGATGCGCCACAGTATAT  >  NC_000913/4627401‑4627874

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