Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,287,779	C→T	intergenic (‑157/+3)	tyrT ← / ← purU	tRNA‑Tyr/formyltetrahydrofolate deformylase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,287,779	0	C	T	96.0%	74.5 / NA	25	intergenic (‑157/+3)	tyrT/purU	tRNA‑Tyr/formyltetrahydrofolate deformylase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (13/11);  minor base A (1/0);  total (14/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.87e-01

GGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGTCATTCACATAGTGAGCGGTTGCGCCAATAATCTTCACACCACGTTCATAGGCCTGGTGATAAG  >  NC_000913/1287547‑1288023                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     gggAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGtt                                                                                                                                                                                                                                    >  2:165242/1‑251 (MQ=255)                     aaggcGATGACTGCAGATTTACAGTCTGCTCCCGTGTGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTGTACTGGCCTGCTCCCTTATCGGGAAGCTGGGCGCATCATATCAAATTACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGTATTAAAGGTAACCAGGTATTTTTGCGGGTGATTAAAGAAGAAGCGTTCGATTACCGTAACCAAAGAc                                                                                                                                                                                                                <  2:20143/247‑1 (MQ=255)                      aGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGAcg                                                                                                                                                                                                               >  2:333030/1‑251 (MQ=255)                      aGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAAACCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCATTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGAcg                                                                                                                                                                                                               >  1:222061/1‑251 (MQ=255)                             gACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCtgtg                                                                                                                                                                                                        <  2:77371/251‑1 (MQ=255)                                                             gCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACg                                                                                                                                                                        <  1:321115/251‑1 (MQ=255)                                                                                     aTGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcg                                                                                                                                                >  1:71292/1‑251 (MQ=255)                                                                                        cTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATAAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGCTTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgcgc                                                                                                                                             >  1:315352/1‑251 (MQ=255)                                                                                                    cTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTACAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGCCTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCCTATCTTc                                                                                                                                  >  2:110996/1‑250 (MQ=255)                                                                                                                                 catATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATacc                                                                                                    >  2:77317/1‑249 (MQ=255)                                                                                                                                 catATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAAc                                                                                                    <  1:110996/251‑1 (MQ=255)                                                                                                                                 catATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTGCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGCGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCAGATCTTCAGCTGTGTATGTATGATCGACATGAATAAc                                                                                                    <  2:387782/251‑1 (MQ=255)                                                                                                                                    atCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGATGCGTCATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGCCAGTACTTTGTATAGTGAACGACTTAAGACGTTTTTCTCGACGTAACGACCTGCGCGCATAATATCTTCAGCTGTGTAGGACTGATCGACATGCATAACgac                                                                                                 >  2:343203/1‑249 (MQ=255)                                                                                                                                        aaTGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGc                                                                                             <  1:198815/251‑1 (MQ=255)                                                                                                                                          tGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCat                                                                                           >  1:136552/1‑251 (MQ=255)                                                                                                                                             cgcgCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatgat                                                                                        <  1:333030/251‑1 (MQ=255)                                                                                                                                                             gtTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTcc                                                                        <  2:15545/251‑1 (MQ=255)                                                                                                                                                               tACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCGGTGTAGGTATGATCGACATGaa                                                                                                        <  1:162441/217‑1 (MQ=255)                                                                                                                                                               tACGTTGAGAAAAATGAACTGGTTGCGAAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTCTTTTTGCTGGTGCTTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGaa                                                                                                        >  2:162441/1‑217 (MQ=255)                                                                                                                                                                  gTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGAtt                                                                   <  2:382012/251‑1 (MQ=255)                                                                                                                                                                                      tGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGTCATTCACATAGTGAGCgg                                               >  1:111887/1‑251 (MQ=255)                                                                                                                                                                                                      gTAACGGATTAAAGGAAACCAGTTATTTTTGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGTCATTCACATAGTGAGCGGTTGCGCCAATAATCtt                               >  1:123283/1‑251 (MQ=255)                                                                                                                                                                                                                          aGTTATTTTTGCTGGAGATTAAAGAATAATCGTTCGATTACCGTACTCATAGACGCGCTGTGCCAGTACTATGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGAGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGTCATTCACATAGTGAGCGGTTGCGCCAATAATCTTCACACCACGATCATAGGCCt           >  2:127582/1‑251 (MQ=255)                                                                                                                                                                                                                                   tGCTGGTGATTAAAGAATAATCGTTCGATTACCGTAACCAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGTCATTCACATAGTGAGCGGTTGCGCCAATAATCTTCACACCACGTTCATAGGCCTGGTGATa    >  1:437075/1‑249 (MQ=255)                                                                                                                                                                                                                                    gCTGGTGATTAAAGAATAATCGTTCGATTACCGGAATAAAAGACACGCCGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCCTCATATCTTCAGATGTGTAGGTATGATCGACCTGAATAACGCCTTGCATGATGATTGGGCCTTCGTCCAGATTGTCATTCACATAGTGAGCGTTTGCGCCAATAATCTTCACACCACGTTCATAGGCCTGGTGATAAg  <  1:351332/250‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                     GGGAAGGATTCGAACCTTCGAAGTCGATGACGGCAGATTTACAGTCTGCTCCCTTTGGCCGCTCGGGAACCCCACCACGGGGTAATGCTTTTTACTGGCCTGCTCCCTTATCGGGAAGCGGGGCGCATCATATCAAATGACGCGCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCCAGATTGTCATTCACATAGTGAGCGGTTGCGCCAATAATCTTCACACCACGTTCATAGGCCTGGTGATAAG  >  NC_000913/1287547‑1288023

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