Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	927,416	Δ1 bp	coding (17/705 nt)	aat ←	leucyl/phenylalanyl‑tRNA‑‑protein transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	927,416	0	A	.	97.1%	139.1 / ‑4.1	34	coding (17/705 nt)	aat	leucyl/phenylalanyl‑tRNA‑‑protein transferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (15/18);  total (16/18)
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

GGGCGAGCGTTTATGAAATCGCTTCATACTACGGCTGATATGCAGTGATTCTGGCCATAGCACCGCGCGGGGATCGGGCGACCACCAGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAG  >  NC_000913/927175‑927647                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        gggCGAGCGTTTATGAAATAGCTTCATACTACGGCTGATATGCAGTGATTCTGGCCATAGCACCGCGCGGGGATCGGGCGCCCACCAGAGGATGGGGTCACATGGAGAAAACCACGGAAAAATAACCCGATGGAAAGCCATTAACAGGCGCGCAGGGCTAAAATCGCCCCCAAGTGCCAGGAGGCCGTTAGGCTCACGTAAGGATCCTTCCGGAGACGGGCAGGCTTATGAATGGCGA‑AAAGCTGAAca                                                                                                                                                                                                                                 >  2:167772/1‑248 (MQ=255)                               acGGCTGATATGCAGTGATTCTGGCCATAGCACCGCGCGGGGATCGGGCGACCACCCGAGGATGGGGTCGCCTGGCGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCg                                                                                                                                                                                                 >  1:196233/1‑251 (MQ=255)                               acGGCTGATATGCAGTGATTCTGGCCATAGCACCGCGCGGGGATCGGGCGACCACCAGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCg                                                                                                                                                                                                 <  2:201669/251‑1 (MQ=255)                                      atatGCAGTGATTCTGGCCATAGCACCGCGCGGGGATCGGGCGACCACCAGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAACCTCCACGCAAGTCGGATCGtt                                                                                                                                                                                          >  2:270505/1‑251 (MQ=255)                                      atatGCAATGATTCTGGCCATAGCACCGCGCGGGGATAGGGCGACCACCAGAGGATGGGGTCGCCTGGAGCAAACAACGGAAAAATAACACGCTGGAAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCAAGAAGGCCGATAGGCTCACGTAATGCGCCTTCAGGGGACGAGAAGGCAATTGACTGGCGAGAA‑GCAGAACCAGGCGCATGACCGCAAAACTACCCGCCAGTCGGAt                                                                                                                                                                                              >  2:83991/1‑247 (MQ=255)                                        atGCAGTGATTCTGGCCATAGCACCGCGCGGGGATCGGGCGACCACCAGATGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCa                                                                                                                                                                                        <  1:281072/251‑1 (MQ=255)                                          gccgTGATTCTGGCCATAGCACCGCGCGGGGATCGGGCGACCACCAGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAGCCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCaat                                                                                                                                                                                      <  1:46718/248‑1 (MQ=255)                                                                cgcgggGGGATCGGGCGACCACCAGAGGATGGGGTCGCCGGGAGAAAGCCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGAGCGCCCCCAAGTGCCAGCAGGCCGGTAGGCGCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGTACGGCGAGAA‑GCTGAACCAGGCGCATGTCCGCAAAACTCCACGCAAGTCGGAGCGTTCAATAATAGCTTACAAACCCTGCTTg                                                                                                                                                                <  2:84537/251‑1 (MQ=255)                                                                     ggggATCGGGCGCCCACCAGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCCGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGACGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTg                                                                                                                                                           >  2:140996/1‑251 (MQ=255)                                                                            gggCGACCACCAGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGCAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTCTTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCCACAATAGCTTACAAACCCTGCTTGAACTggcacta                                                                                                                                                     >  2:50854/1‑245 (MQ=255)                                                                                gACCACCAGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGcc                                                                                                                                                <  1:287362/251‑1 (MQ=255)                                                                                 accaccAGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGccc                                                                                                                                               >  2:234855/1‑251 (MQ=255)                                                                                    accaGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATTCCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTTCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGTCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTg                                                                                                                                            <  1:103923/251‑1 (MQ=255)                                                                                              ggggTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGc                                                                                                                                  <  2:252127/251‑1 (MQ=255)                                                                                                        ggAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGt                                                                                                                        <  1:291066/251‑1 (MQ=255)                                                                                                         gagaAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCATGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCGGCTTGAACTGTTAAGAACGCCCCTGTCTGGCAAGCAGTTCTGCGTg                                                                                                                       <  1:196838/251‑1 (MQ=255)                                                                                                                                 gCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGccc                                                                                               >  1:200217/1‑251 (MQ=255)                                                                                                                                  cTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCg                                                                                              >  1:73025/1‑251 (MQ=255)                                                                                                                                   ggctAAGCCATTAACACGCGCCCAGGGCTAAGCTCGCCCCCAAGGGCCAGCAGGCCGTTAGGCTCACTTAATGCGCCTTCCGGGGTAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGTCCGCACAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAATCCCTGCTTCATCTGCTAATTACTCCCCTGTCTGGCTAGCAGTTCTGCGTGAGTACTTGGCTCAATAATTTGCCCGt                                                                                             <  2:50443/248‑1 (MQ=255)                                                                                                                                               aaCAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACta                                                                                 <  2:145983/251‑1 (MQ=255)                                                                                                                                               aaCAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACt                                                                                  <  1:140996/250‑1 (MQ=255)                                                                                                                                                   ggCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACtattat                                                                             >  1:118520/1‑251 (MQ=255)                                                                                                                                                             gCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCAAGTAATGCGCCTGCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGaaa                                                                   <  1:292083/251‑1 (MQ=255)                                                                                                                                                                  gATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACg                                                                 >  1:277119/1‑248 (MQ=255)                                                                                                                                                                                   aacgggCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATTGCGAGAA‑GATGAACCAGGCGCATGACCGCAAAACTCCACGCAATTCGGATCGTTCAATAATAGCTTACAATCCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTGCTGCGGGAGTTCCTTGTTCAATAATTTGCCCGTTGTCCATCACTGTTATGTGTTGGAAACTAGATAGTCGGCGAATTCGAt                                             <  2:23491/247‑1 (MQ=255)                                                                                                                                                                                    gcagGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATg                                            <  1:119381/251‑1 (MQ=255)                                                                                                                                                                                                     cGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAC‑GCTGACCCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGCTCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCg                           >  1:186609/1‑251 (MQ=255)                                                                                                                                                                                                          tGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATc                                                                                     <  2:193377/188‑1 (MQ=255)                                                                                                                                                                                                          tGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATc                                                                                     >  1:193377/1‑188 (MQ=255)                                                                                                                                                                                                              ccTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCAc                  >  2:68307/1‑251 (MQ=255)                                                                                                                                                                                                              ccTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCAc                  >  1:105426/1‑251 (MQ=255)                                                                                                                                                                                                                         agggaaggCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCt       >  1:267679/1‑251 (MQ=255)                                                                                                                                                                                                                             aaggCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTgcaa   <  1:234855/251‑1 (MQ=255)                                                                                                                                                                                                                              aggCTATTGAATGGCGAGAA‑GCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTTCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTTACCATTAACACCGTTTTCTCACGCATCATTTCTgcaag  <  1:147301/251‑1 (MQ=255)                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        GGGCGAGCGTTTATGAAATCGCTTCATACTACGGCTGATATGCAGTGATTCTGGCCATAGCACCGCGCGGGGATCGGGCGACCACCAGAGGATGGGGTCGCCTGGAGAAAACCACGGAAAAATACCACGCTGGTAAGCCATTAACAGGCGCGCAGGGCTAAGATCGCCCCCAAGTGCCAGCAGGCCGTTAGGCTCACGTAATGCGCCTTCCGGGGAAGGGAAGGCTATTGAATGGCGAGAAAGCTGAACCAGGCGCATGACCGCAAAACTCCACGCAAGTCGGATCGTTCAATAATAGCTTACAAACCCTGCTTGAACTGGTAATAACGCCCCTGTCTGGCAAGCAGTTCTGCGTGAGTACCTTGCTCAATAATTTGCCCGTTGTCCATCACTATTATTTGTTGGAAACGAGAGAGTCCGCGAAGTCGATGGGTGACCATTAACACCGTTTTCTCACGCATCATTTCTGCAAG  >  NC_000913/927175‑927647

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