Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,275,902	A→T	E207D (GAA→GAT)	garD →	GarD

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,275,902	0	A	T	95.7%	69.6 / ‑5.3	23	E207D (GAA→GAT)	garD	GarD
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (7/15);  minor base G (1/0);  total (8/15)
Fisher's exact test for biased strand distribution p-value = 3.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.65e-01

CGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGAAGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTT  >  NC_000913/3275659‑3276133                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         ggATGGCAGCGTGTGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATgg                                                                                                                                                                                                                                   <  2:252086/250‑1 (MQ=255)             gggCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGTTGATTGGCCTgg                                                                                                                                                                                                                       <  2:298477/251‑1 (MQ=255)                           ggCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTg                                                                                                                                                                                                        <  2:221053/250‑1 (MQ=255)                               cGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCa                                                                                                                                                                                                      <  2:160039/249‑1 (MQ=255)                                ggTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGCACTCTGTACCAAAAATCATTCAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGAGGGGCTGAATCCTTTGTACGGTTGTGGCGTGGCGATTAACCCACCGGCGGCAG‑GTGTACCTATCCCTACCACTCACCATATTTCGCTGAATCCTAAATTTGGCGGCGCTGTAATGGGGATTGGCCTGGGTTGT‑GGAAAGTTGCAGcc                                                                                                                                                                                                   >  1:175023/1‑251 (MQ=255)                                                             aGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGGTTAACGCACCGGCGGCAGGTTG‑ACCTTTCCCTTCCCTTTACAATATTTCTCTGAATTCTAACTTTTGCGGGGTGGTTATGGTGGTGGGTCTGGGTTGTTGAAAAGTTGC‑GCCTGAGCGCCTGCTGAATGGAAAGgatggtt                                                                                                                                                                      >  2:244906/1‑248 (MQ=255)                                                                   ggtggACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCaa                                                                                                                                                                 <  2:245016/250‑1 (MQ=255)                                                                      ggACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTa                                                                                                                                                             <  2:159840/251‑1 (MQ=255)                                                                                   aaaaTCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGGTTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAg                                                                                                                                                <  1:145378/251‑1 (MQ=255)                                                                                        cATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAAAATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCa                                                                                                                                           <  2:210745/251‑1 (MQ=255)                                                                                           tGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCCG‑TTGTACCTATCCGTACCATTCCCAATATTTCGCTGAATCCTAACTTTGGCGGCGATGGAATGGTGATTGGCCTGGGGTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGc                                                                                                                                         >  1:253733/1‑250 (MQ=255)                                                                                                gcgATCTGCGACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTGCCTATCCGTACCATTCACAGTATTTCGCTGAATCCTAGCTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTc                                                                                                                                   <  1:125868/251‑1 (MQ=255)                                                                                                  gATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGGTTGGCCTGGGGTTT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGGACGGATGATGTGCAAGCTTTTCCAGTAGAAAGCGCCAGCAttttc                                                                                                                                   >  1:163760/1‑246 (MQ=255)                                                                                                            cGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGGAAGCTATTCCAGTAGAAAGCGCCAGCATTGGCAGTTGGCAGGAt                                                                                                                       <  1:246204/251‑1 (MQ=255)                                                                                                                        aCGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCAACGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCCGGGTTTT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCCGTTTGCAGGATGAAAAGCATGTc                                                                                                           >  1:267527/1‑251 (MQ=255)                                                                                                                                     ggtggGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTcc                                                                                               <  2:9396/250‑1 (MQ=255)                                                                                                                                                                                                gTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGc                                     >  1:177847/1‑250 (MQ=255)                                                                                                                                                                                                 tACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAACAACTGAATCAACGGc                                     >  1:126258/1‑249 (MQ=255)                                                                                                                                                                                                  cccTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAAGCAACGGCAGc                                  <  2:169885/250‑1 (MQ=255)                                                                                                                                                                                                  aCCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAg                                   <  2:255342/250‑1 (MQ=255)                                                                                                                                                                                                                        atTTCGCTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAg             <  2:251661/250‑1 (MQ=255)                                                                                                                                                                                                                              cTGAATCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGt      >  1:298304/1‑251 (MQ=255)                                                                                                                                                                                                                                  aTCCTAACTTTGGCGGCGATGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGtt  <  2:93128/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                         CGATGGCAGCGTGGGCACCAAAAACCTGCTCGGTATCACCACCAGCGTCCACTGTGTGGCAGGCGTGGTGGACTATGTAGTAAAAATCATTGAACGCGATCTGCTACCGAAATACCCGAACGTCGATGGCGTGGTGGGGCTGAATCATTTGTACGGTTGTGGCGTGGCGATTAACGCACCGGCGGCAG‑TTGTACCTATCCGTACCATTCACAATATTTCGCTGAATCCTAACTTTGGCGGCGAAGTAATGGTGATTGGCCTGGGTTGT‑GAAAAGTTGCAGCCTGAGCGCCTGCTGACTGGAACGGATGATGTGCAAGCTATTCCAGTAGAAAGCGCCAGCATTGTCAGTTTGCAGGATGAAAAGCATGTCGGTTTTCAGTCCATGGTCGAGGATATTTTGCAGATCGCCGAACGCCATCTACAAAAACTGAATCAACGGCAGCGAGAAACCTGCCCGGCTTCAGAACTGGTCGTT  >  NC_000913/3275659‑3276133

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