Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,229,894	C→T	E164E (GAG→GAA)	dusC ←	tRNA‑dihydrouridine(16) synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,229,894	0	C	T	96.7%	95.8 / NA	30	E164E (GAG→GAA)	dusC	tRNA‑dihydrouridine(16) synthase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (11/18);  minor base A (1/0);  total (12/18)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGACGGGCAGCCGCAATTGAGATCCACGCCCCAG  >  NC_000913/2229653‑2230119                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  gCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAgcg                                                                                                                                                                                                                          >  2:142293/1‑251 (MQ=255)       gcgGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGcc                                                                                                                                                                                                                    >  1:273082/1‑251 (MQ=255)                tataTTGTACCACCCGGCTCAGGTTGGGAATATGGAGCGCCCCGCGACCAAGCATCAATGCGTCGCATCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTGTGCGATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCGTTCTTCCGCCCATGCACCCCCAGTTCAGTAGTGCCAGCTTGTTGAAc                                                                                                                                                                                                            <  2:198531/250‑1 (MQ=255)                   atTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCGCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATAGCGCCATTCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCa                                                                                                                                                                                                        <  2:273082/251‑1 (MQ=255)                         ccaccCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCg                                                                                                                                                                                                  <  1:225628/251‑1 (MQ=255)                                        tGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATttcttc                                                                                                                                                                                   >  1:92569/1‑251 (MQ=255)                                         gggAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTtcttct                                                                                                                                                                                  >  2:36436/1‑251 (MQ=255)                                                 tGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCt                                                                                                                                                                           <  1:2893/250‑1 (MQ=255)                                                    gcgcCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTcc                                                                                                                                                                       <  1:142293/251‑1 (MQ=255)                                                         ccGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGTATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTATCCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAATCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTcccagcc                                                                                                                                                                  <  1:192926/251‑1 (MQ=255)                                                          cgcgACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTcccagccc                                                                                                                                                                 >  2:240494/1‑251 (MQ=255)                                                                   tcatcaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACt                                                                                                                                                        >  1:128396/1‑251 (MQ=255)                                                                   tcatcaCTGCGTCGCAGCCACTGAGCGCCATGCCGTGTTGCGCGCTCTGCCGGTACCAGATTTCACCAGTAGCAATCACCGGAATATTCAGCCGCTAGCGAATATCACCGATCGCCGGCCAGTCAATATGCTCCTCGCGGTAACCCTGCGCGTTCATCCGCCCAGGCACCGCCAGATCCGTAACGCCAGCCTGTTGAACCGCCTCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCACACGCACt                                                                                                                                                        >  2:242121/1‑251 (MQ=255)                                                                      ccaCTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTc                                                                                                                                                     <  2:238477/250‑1 (MQ=255)                                                                        ccTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCCCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCAc                                                                                                                                                   <  2:1218/250‑1 (MQ=255)                                                                           gCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACg                                                                                                                                                  <  1:198268/249‑1 (MQ=255)                                                                               cGCAGCCGCTGATCGCCATGCATGGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACg                                                                                                                                            <  2:92569/251‑1 (MQ=255)                                                                                  aGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGc                                                                                                                                         >  1:49337/1‑251 (MQ=255)                                                                                      gCTGATCGCCATGCATTGTTGCGCGATCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGTCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGACAGCCTGTTGAACCGCATCGGCGACTTCAAATTCCTGCTCACCGCTGTCCCAGCCCAGACGCACGGTCACGCTGACGGGCAAAt                                                                                                                                     <  1:242121/251‑1 (MQ=255)                                                                                                                        tCCCAGATTTCCCAGTTGGCAATCACCGGAATATTCATCCGCTGCCGCATATAGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCCCTTTCACGCTGACGGGCAAATGCGCCGGTACCGCTTCACGCATCGCTTTTGCAcc                                                                                                   >  2:157653/1‑251 (MQ=255)                                                                                                                           cAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTg                                                                                                <  1:216444/251‑1 (MQ=255)                                                                                                                                          gCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAgg                                                                                 >  1:201307/1‑251 (MQ=255)                                                                                                                                                            aGCCGCTGGCGAATATCGCCAATGGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGCACAGCTTCACGCATCGCTTTTGCACCCTTGTAGATGTGTTCAGGATCTTTGAGTAACGTctc                                                               >  2:200394/1‑249 (MQ=255)                                                                                                                                                              ccGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGccc                                                             <  1:208768/251‑1 (MQ=255)                                                                                                                                                                                 aTCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAAc                                          <  1:36436/251‑1 (MQ=255)                                                                                                                                                                                           aGTCAATAGGCTCCGCGCGGTAACCCCGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCTCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGACTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGAc                                <  2:188103/251‑1 (MQ=255)                                                                                                                                                                                                    gCTCCGCGCGGTAACCCTGAGCTGTCGCCCGCCCATGCACCACCAGTGCGGTAGCGCCAGCCGGTTGAAGCGCATCGGCGATTTCAAATGTCTTCTCACCGCTGCCCCAGCCCAGACGCAGTTGCGCGCTGACGGGCAAAGGCGCGGGTAAAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGAGCTTTGAGTAACGTCGCACCGCCGCCGCTACCGTTAACCGTTTGCGACGGGCAGcc                        <  2:52925/250‑1 (MQ=255)                                                                                                                                                                                                       ccGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGTTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGCCGGGCAGCCGCa                     >  2:128693/1‑250 (MQ=255)                                                                                                                                                                                                                        ccttacgcccGCCCACGCACCACCAGTTCCGTCGCGCCAGCCTGTTGACCCTACTCGGCGATTTAAAATTACTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTATATGGGTTCAGGATTTTTGAGTAACGTCTCCCCGCCGCCGCTACCTTTTACCGTTTTGGACGGGCAGCCGCAATTGAGATCCACGCCCCa   <  1:89781/243‑1 (MQ=255)                                                                                                                                                                                                                          ttccgtcacccagtgtaccaccAGTTCCGTGGCGCCAGCCTGGTGAACCGAATCGCGGATGTCAAATTGCTTCTCGCCGCTGGCCCAGCCCAGACGCCCTTTCACGCTGACGGGCAAATGCGCCGGTAAAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCTCTCCGCCGCCGCTACCGTTAACCGTTTTCGACGGGCAGCCGCAATTGAGATCCACGCCCCAg  <  2:128396/234‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                  GCATTCGCGGTTCGTTATATTTTACCACCCGGCTCAGGTTGGGAATATTGAGCGCCCCGCGACCAATCATCACTGCGTCGCAGCCGCTGATCGCCATGCATTGTTGCGCGCTCTGCCAGTCCCAGATTTCACCGTTGGCAATCACCGGAATATTCAGCCGCTGGCGAATATCGCCAATCGCCTGCCAGTCAATATGCTCCGCGCGGTAACCCTGCTCTTTCGTCCGCCCATGCACCACCAGCTCCGTAGCGCCAGCCTGTTGAACCGCATCGGCGATTTCAAATTTCTTCTCACCGCTGTCCCAGCCCAGACGCACTTTCACGCTGACGGGCAAATGCGCCGGTACAGCTTCACGCATCGCTTTTGCACCCTGGTAGATGAGTTCAGGATCTTTGAGTAACGTCGCCCCGCCGCCGCTACCGTTAACCGTTTTCGACGGGCAGCCGCAATTGAGATCCACGCCCCAG  >  NC_000913/2229653‑2230119

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