Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,725,236	A→G	V28A (GTT→GCT)	glyQ ←	glycine tRNA synthetase, alpha subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,725,236	0	A	G	100.0%	171.1 / NA	46	V28A (GTT→GCT)	glyQ	glycine tRNA synthetase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (22/24);  total (22/24)

AAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGA  >  NC_000913/3725089‑3725379                                                                                                                                                    |                                                                                                                                                aaaCGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCaa                                                                                                                                              <  1:379664/151‑1 (MQ=255) aaaCGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCa                                                                                                                                               >  1:23329/1‑150 (MQ=255) aaaCGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCa                                                                                                                                               <  1:441613/150‑1 (MQ=255)   aCGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATg                                                                                                                                            >  2:647897/1‑151 (MQ=255)   aCGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATg                                                                                                                                            >  1:564436/1‑151 (MQ=255)        tGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGc                                                                                                                                        >  1:538278/1‑150 (MQ=255)         ggggTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGc                                                                                                                                        >  2:270892/1‑149 (MQ=255)         ggggTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAg                                                                                                                                      <  2:501078/151‑1 (MQ=255)               ttCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTgg                                                                                                                                >  1:612513/1‑151 (MQ=255)                tCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGc                                                                                                                               >  2:387459/1‑151 (MQ=255)                              aTCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATcc                                                                                                                 >  2:419693/1‑151 (MQ=255)                                cggtcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTg                                                                                                               >  1:88545/1‑151 (MQ=255)                                  gtcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGAGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGt                                                                                                              >  1:499728/1‑150 (MQ=255)                                   tcggGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTa                                                                                                             >  1:625912/1‑150 (MQ=255)                                      ggCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGt                                                                                                         <  2:432928/151‑1 (MQ=255)                                         gacgaGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAg                                                                                                      <  2:178463/151‑1 (MQ=255)                                           cgaGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGa                                                                                                    <  1:204886/151‑1 (MQ=255)                                                  gCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGcc                                                                                             <  2:477755/151‑1 (MQ=255)                                                      aaCATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTgg                                                                                         <  1:554297/151‑1 (MQ=255)                                                        cATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGaa                                                                                       <  1:422220/151‑1 (MQ=255)                                                         aTAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAg                                                                                      <  2:491078/151‑1 (MQ=255)                                                               agcCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTg                                                                                >  2:221405/1‑151 (MQ=255)                                                                gccgccATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGACCAAGCCCTGGAAGGTCCTgg                                                                               <  1:644932/151‑1 (MQ=255)                                                                   gccATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTAt                                                                            >  1:127285/1‑151 (MQ=255)                                                                     cATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATc                                                                           >  1:598364/1‑150 (MQ=255)                                                                        cGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCaaa                                                                        <  2:161695/150‑1 (MQ=255)                                                                        cGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAAc                                                                       <  1:381216/151‑1 (MQ=255)                                                                              tGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGc                                                                 >  2:386047/1‑151 (MQ=255)                                                                                    cAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCAtatta                                                           <  1:482614/151‑1 (MQ=255)                                                                                           cgcAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTg                                                    <  2:278806/151‑1 (MQ=255)                                                                                            gcAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGc                                                   >  2:269903/1‑151 (MQ=255)                                                                                             cAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCt                                                  <  2:179684/151‑1 (MQ=255)                                                                                                   gggcATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATa                                             <  1:348217/147‑1 (MQ=255)                                                                                                         tGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGtgtg                                       <  2:258733/150‑1 (MQ=255)                                                                                                            gTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGAtt                                   <  2:496320/151‑1 (MQ=255)                                                                                                             tGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGAttt                                  >  2:487935/1‑151 (MQ=255)                                                                                                                  ggTTCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAaagg                             <  2:302871/151‑1 (MQ=255)                                                                                                                     tCCCGCGCCGACTTCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAaagga                            >  1:351512/1‑149 (MQ=255)                                                                                                                                 ttCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAc              >  2:72397/1‑151 (MQ=255)                                                                                                                                 ttCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAc              >  1:118020/1‑151 (MQ=255)                                                                                                                                 ttCCATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAc              >  1:565770/1‑151 (MQ=255)                                                                                                                                   ccATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATAccc            <  2:516950/151‑1 (MQ=255)                                                                                                                                    cATGTCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCg           <  1:56313/151‑1 (MQ=255)                                                                                                                                        tCCAATGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCt         <  1:565584/149‑1 (MQ=255)                                                                                                                                           aaTGGTTGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCGGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtggatg    <  2:329335/151‑1 (MQ=255)                                                                                                                                                 tGAGCAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCtggatgga  >  1:541930/1‑147 (MQ=255)                                                                                                                                                    |                                                                                                                                                AAACGGTTGGGGTTTTCGCCGTAGCGACCATCGGTCGGGCGACGAGAAGGCTGAACATAAGCAGCCGCCATCGGTTCTGGCCCCAGCTCGCGCAGACAGGTCATTGGGTGAGAGGTTCCCGCGCCGACTTCCATGTCCAATGGTTGAACAATGGTGCAGCCCTGGCGAGCCCAGTAATCCTGTAAGGTCAGGATCAAGCCCTGGAAGGTCCTGGTATCAAACTTTTGCATATTATTTCGTGCTGGATACGTGTGGATTTAAAGGAAGGGATCAGTATACCCGCTGGATGGA  >  NC_000913/3725089‑3725379

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