Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	CP009273	4,051,806	G→T	100%	V14L (GTA→TTA)	typA →	GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	CP009273	4,051,806	0	G	T	100.0%	88.0 / NA	29	V14L (GTA→TTA)	typA	GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base T (14/15);  total (14/15)
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGGCTGCACCACGATGTGAATGTTGCACCAATATAGTGCTTCAATGGAAACATTAAGCACCATGTTGGTGCAATGACCTTTGGATAACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGT  >  CP009273/4051544‑4052034                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     gggCTGCACCACGATGTGAATGTTGCACCAATATAGTGCTTCAATGGAAACATTAAGCACCATGTTGGTGCAATGACCTTTGGATAACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCTAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACAAttgaa                                                                                                                                                                                                                          >  1:146/1‑271 (MQ=255)          anagATGTTAATGTTGAACCAATATAGTGCTTCAATGGAAACATTAAGCACCATGTTGGTGCAATGACTTTTGGATAACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAacca                                                                                                                                                                                                                    <  2:146/269‑1 (MQ=255)            cGATGTGAATGTTGCACCAATATAGTGCTTCAATGGAAACATTAAGCACCATGTTGGTGCAATGACCTTTGGATAACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATAGTAAAACCACCCt                                                                                                                                                                                                                >  1:83406/1‑274 (MQ=255)                     ttttGCATCAATATTGTGCTTCAATGGATACATTAAGCACCATGTTGGTGTAATTACCTTTTGTTAACTCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAACCACCCTGGTAGACAAg                                                                                                                                                                                                      <  1:21307/273‑1 (MQ=255)                                       cTTCAATGGAAACATAAAGCACCATGTTGGTGCAATGACCTTTGGATAACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTa                                                                                                                                                                                                                                    >  1:365230/1‑227 (MQ=255)                                           aaTGGAAACATTAAGCACCATGTTGGTGCAATGACCTTTGGATAACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGAcc                                                                                                                                                                                                                                <  2:365218/227‑1 (MQ=255)                                                                                                   cTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGAcc                                                                                                                                                                                                                                >  1:246800/1‑171 (MQ=255)                                                                                                   atCCGTTAAAGCTATAACAAAGGGACTATTATATACTTTGTTGATGAGAATGTTTTTTATCTTGTTTTGTACTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGGCTGAGGCAAAGTTGTGATCGAAAGATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAACCACCCTGGTAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGaa                                                                                                                         <  1:132849/273‑1 (MQ=255)                                                                                                       gtgtatgCGATCACCAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGGTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGCGGCAAAGTTGTGAACGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATgg                                                                                                                                                                                                                            <  2:246791/165‑1 (MQ=255)                                                                                                                             tctGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGt                                                                                                                                                                                                                           >  2:185152/1‑150 (MQ=255)                                                                                                                             tctGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAACCACCCTGGTAGACAAg                                                                                                                                                                                                      >  1:314137/1‑171 (MQ=255)                                                                                                                               tGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAac                                                                                                                                                                                                                      >  1:211062/1‑153 (MQ=255)                                                                                                                                 cAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTaaaa                                                                                                                                                                                                                       <  1:185159/150‑1 (MQ=255)                                                                                                                                 cAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTGTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAACCACCCTGGTAGACAAgctgc                                                                                                                                                                                                  <  2:314126/171‑1 (MQ=255)                                                                                                                                   aTACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAaccacc                                                                                                                                                                                                                  <  2:211053/153‑1 (MQ=255)                                                                                                                                         gTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAACCACCCTGGTAGa                                                                                                                                                                                                          >  2:295885/1‑155 (MQ=255)                                                                                                                                             gCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGAACATGGTAAAACCACCCTGGTAGACAAg                                                                                                                                                                                                      <  1:295896/155‑1 (MQ=255)                                                                                                                                                                  gTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAACCACCCTGGTAGACAAGCTGCTCCAACAATCCGGTNCGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGAATTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTaac                                                            >  1:419715/1‑270 (MQ=255)                                                                                                                                                                             cGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGt                                                                                                                                                                                                                           >  2:343941/1‑102 (MQ=255)                                                                                                                                                                                 tAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTaaaa                                                                                                                                                                                                                       <  1:343953/102‑1 (MQ=255)                                                                                                                                                                                       gtgtACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAACCACCCTGGTAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATAGAATGAGTACCGTATCCACATCCTTCATACCCCGGTGcct                                    >  2:147393/1‑273 (MQ=255)                                                                                                                                                                                                  cgcgCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAacc                                                                                                                                                                                                                     >  1:156545/1‑87 (MQ=255)                                                                                                                                                                                                      cTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAACCAccc                                                                                                                                                                                                                 <  2:156542/87‑1 (MQ=255)                                                                                                                                                                                                                     gAGGCAAAGTGGTGATCTGAAAATTGCGTAATCTCACCATCATCGCGCACTTAGACCATGGGAAAACCACCCGGGTAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTa      <  1:147395/275‑1 (MQ=255)                                                                                                                                                                                                                         gaaaGTTGTTATCGAAAAATTGCGTAATATCTCCATCATCGCGCACTTAGACCATGGTAAAACCAGCCTGGTAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGGGTGATGGACTGCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCAAGCCGACTTCGGTGGTGAAGTTGTACGTGTAATGt  <  1:314734/274‑1 (MQ=255)                                                                                                                                                                                                                                 tgATCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTaaaa                                                                                                                                                                                                                       >  1:47048/1‑54 (MQ=255)                                                                                                                                                                                                                                    tCGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGAc                                                                                                                                                                                                                                 <  2:117924/41‑1 (MQ=255)                                                                                                                                                                                                                                     cGAAAAATTGCGTAATATCGCCATCATCGCGCACTTAGACCATGGTAAAaccac                                                                                                                                                                                                                   <  2:47048/54‑1 (MQ=255)                                                                                                                                                                                                                                                                    cACTTAGACCATGGTAAAACCACCCTGGTAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGt                                                                                                                                                 >  1:50702/1‑89 (MQ=255)                                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     GGGCTGCACCACGATGTGAATGTTGCACCAATATAGTGCTTCAATGGAAACATTAAGCACCATGTTGGTGCAATGACCTTTGGATAACCCTTTTTATGCTCCGTGAAAGCGATCACAAAGGGACTCTGCAATACTTGTTTGCGGAGGATGTTTGTGATCCTGTTTTGTAGTGCGATTAATCCGTGTACAATAACGCGCTATTTCTAATGCCTGAGGCAAAGTTGTGATCGAAAAATTGCGTAATATCGCCATCATCGCGCACGTAGACCATGGTAAAACCACCCTGGTAGACAAGCTGCTCCAACAATCCGGTACGTTCGACTCTCGTGCCGAAACCCAAGAGCGCGTGATGGACTCCAACGATTTGGAGAAAGAGCGTGGGATTACCATCCTCGCGAAAAACACCGCTATCAAATGGAATGATTACCGTATCAACATCGTTGATACCCCGGGGCACGCCGACTTCGGTGGTGAAGTTGAACGTGTAATGT  >  CP009273/4051544‑4052034

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