Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	218,459	A→G	9.3%	L106P (CTC→CCC)	proS ←	prolyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	218,459	0	A	G	9.3%	122.6 / 4.4	43	L106P (CTC→CCC)	proS	prolyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base A (22/17);  new base G (2/2);  total (24/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAT  >  NC_000913/218311‑218595                                                                                                                                                     |                                                                                                                                         gAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGa                                                                                                                                          <  1:79123/149‑1 (MQ=255) gAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGa                                                                                                                                          <  2:427674/149‑1 (MQ=255)     tCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTAc                                                                                                                                      <  2:24631/149‑1 (MQ=255)        cgGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGaa                                                                                                                                   <  1:325482/149‑1 (MQ=255)           gAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACgg                                                                                                                                >  3:127686/1‑149 (MQ=255)                 aTGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGGGTACGAACGGACGCTc                                                                                                                          >  2:4332/1‑149 (MQ=255)                 aTGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGGGTACGAACGGACGCTc                                                                                                                          >  2:275834/1‑149 (MQ=255)                 aTGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGAAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTc                                                                                                                          <  2:342116/149‑1 (MQ=255)                         gAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGt                                                                                                                  <  3:59261/149‑1 (MQ=255)                                 gACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACaaaa                                                                                                          <  2:130374/149‑1 (MQ=255)                                    gCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACgc                                                                                                       >  3:87266/1‑149 (MQ=255)                                            tcgCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTc                                                                                               <  4:111607/149‑1 (MQ=255)                                               cgGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCgg                                                                                            >  1:298116/1‑149 (MQ=255)                                                ggAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGa                                                                                           >  2:534556/1‑149 (MQ=255)                                                     ttGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTa                                                                                      >  4:260024/1‑149 (MQ=255)                                                           tGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCTGGTCAGTGATAACTGCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCTGACCGTACTGTTc                                                                                <  4:491017/149‑1 (MQ=255)                                                            ggATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTcc                                                                               <  2:143927/149‑1 (MQ=255)                                                                  gATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGa                                                                         <  4:42533/149‑1 (MQ=255)                                                                     aGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCa                                                                      <  2:425037/149‑1 (MQ=255)                                                                     aGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCa                                                                      >  1:80155/1‑149 (MQ=255)                                                                     aGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCa                                                                      >  3:176160/1‑149 (MQ=255)                                                                                    gCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACAGTTCCCAACGACCACTCTCTTGCCACaaa                                                       <  1:464742/149‑1 (MQ=255)                                                                                      tGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATc                                                     >  4:374470/1‑149 (MQ=255)                                                                                      tGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATc                                                     <  4:310841/149‑1 (MQ=255)                                                                                       gTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCg                                                    <  4:76647/149‑1 (MQ=255)                                                                                          tGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCggct                                                 <  1:141431/149‑1 (MQ=255)                                                                                              agagCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCggctggct                                             >  3:414130/1‑149 (MQ=255)                                                                                                 gCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGaa                                          >  3:270515/1‑149 (MQ=255)                                                                                                  cTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAac                                         <  2:365158/149‑1 (MQ=255)                                                                                                          gTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGTCCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCa                                 <  4:327244/149‑1 (MQ=255)                                                                                                           ttACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCAt                                >  4:307448/1‑149 (MQ=255)                                                                                                           ttACGAAGCGGGTCGGTGAGAACTTCTTCATGAGGTGGGTCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACGGTACTGTTACCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGCCAt                                >  4:44712/1‑149 (MQ=255)                                                                                                            tACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATc                               <  3:374470/149‑1 (MQ=255)                                                                                                                 aTCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGacac                          >  4:470737/1‑149 (MQ=255)                                                                                                                 aTCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGacac                          >  4:99748/1‑149 (MQ=255)                                                                                                                  tCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACAcc                         >  2:551174/1‑149 (MQ=255)                                                                                                                     ggTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCTCTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCg                      >  4:384167/1‑149 (MQ=255)                                                                                                                       tCAGTGATAACTTCTTCATGAGTTGGGCCGGGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGAt                    <  1:4332/149‑1 (MQ=255)                                                                                                                       tCAGTGATAACTTCTTCATGAGTTGGGCCGGGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGAt                    <  1:275834/149‑1 (MQ=255)                                                                                                                       tCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGAt                    >  2:318555/1‑149 (MQ=255)                                                                                                                            gATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCAc               >  1:469060/1‑149 (MQ=255)                                                                                                                               aaCTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCgg            >  3:246432/1‑149 (MQ=255)                                                                                                                                  ttcttcATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCgt         >  1:509252/1‑149 (MQ=255)                                                                                                                                  ttcttcATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCgt         >  1:372154/1‑149 (MQ=255)                                                                                                                                         tGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAt  >  4:228246/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                         GAATTCGCGGGAACGCATGACGCCGAAACGCGGACGCACTTCGTCGCGGAACTTGGTCTGGATCTGATAGAAGTTCAGCGGCAGCTGTTTGTAAGAGCTAAGCTCGTTACGAATCAGGTCAGTGATAACTTCTTCATGAGTTGGGCCGAGTACGAACGGACGCTCGCCACGGTCAACAAAACGCAGCAGTTCCGGACCGTACTGTTCCCAACGACCACTCTCTTGCCACAAATCGGCTGGCTGAACCACCGGCATCGACACCTCGATCGCACCGGCGTTGTTCAT  >  NC_000913/218311‑218595

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