Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB122	3,626,095	T→G	47.8%	G170G (GGT→GGG)	hcf136_3 →	Ycf48‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB122	3,626,095	0	T	G	47.8%	14.4 / 4.3	26	G170G (GGT→GGG)	hcf136_3	Ycf48‑like protein
Reads supporting (aligned to +/- strand):  ref base T (9/4);  new base G (2/11);  total (11/15)
Fisher's exact test for biased strand distribution p-value = 1.54e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGC  >  ACB122/3625972‑3626223                                                                                                                            |                                                                                                                                 ggcagCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGt                                                                                                                       >  6:454085/1‑135 (MQ=255)  gcagCCGCGGCTGCCGCAGACGAGCGGCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTg                                                                                                                      <  8:390417/135‑1 (MQ=255)    agCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGGGGCGACGGGGTGGc                                                                                                                    >  3:380579/1‑135 (MQ=255)         aGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGGTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATgc                                                                                                               >  4:331519/1‑135 (MQ=255)          ggCTGCCGCAGACGAGCGTCGCCGGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGGCCTGAATTTCACCGCTGCCCGCAACGGCCTGGTGGTGGGGGGCGACGGGGGGGGCCTgcg                                                                                                              >  5:342185/1‑135 (MQ=255)                                                          cgcgggCGCGCACAAGCGTTTGTGGGCCCTGAATTGCCCCGATCCCCGCAACGGCGTGGGCGGGGGGGCCTCCGGGGTGCCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                              <  8:224449/131‑1 (MQ=255)                                                          cgccGGCGCGGACACGCGGTTGTTGCCCCTCAATTTCCCCGAGCCCCCCACCGGCCGGGGCGGCGGGGCCTCGGGGGGGCCCAGGCGCACCCACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                              <  7:467193/132‑1 (MQ=255)                                                          cGAGGGCGCGGACACGCCGTTGTGGGCCCTGAATTCCCCCGATCCCCCCACCGGCGTGGCCGGCGGGGCCTACGGGGGGCCCAGGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGcc                                                              <  3:357732/135‑1 (MQ=255)                                                               gcgcGGACAAGCGGTTGTTGGCCCTGAATTGCACCGCTGCCCGCACGGGCGTGGGGGGGGGGGCCTCGGGGGTGGACTGCCCCACCCACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACg                                                         <  3:179300/135‑1 (MQ=255)                                                                  cgcacacGCGGTTGTTGGCCCTTAATTACCCCTATGCCCGCACGGGCCGGGGCGGCGGGGCCTCGGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGc                                                      <  6:206940/132‑1 (MQ=255)                                                                    caCAAGCGTTTGTGGGCCCTTAATTTCCCCTAGCCCCGCAACGGCGTGGGGGGCGGGGCCTAGGGGGGGGCCAGGCGCACCTACGACGGCGGGCGCAGCTGGCAGGCGGAGATCGGCCGACTGCCGAACGGGCGt                                                    <  8:204/134‑1 (MQ=255)                                                                       acgcggttgtgGGCCCTGAATTCCCCCGAGGCCCGCACGGGCGGGGTCGGCGGGGCCTAGGGGGGGCCCAGGCCCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGa                                                 <  4:269900/130‑1 (MQ=255)                                                                         gcggttgtgGGCCCTGATTTTCCCGGTTGCCCGCACCGGCCGGGGCGGCGGGGCCTCCGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACt                                               <  4:121353/132‑1 (MQ=255)                                                                              ttttgGCCCTTAATTTCCCGGAGGCCCGCACCGGCGGGGGGGTGGGGCCCTCGGGGGGGGCCTGGCCCCCCGCGGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCAcc                                          <  6:416246/133‑1 (MQ=255)                                                                              tgttgGCCGTGAATTGCACCGATGCCCGCACCGGCGGGGGCGGGGGTGCCTCCGGGGGGGCCTGGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCAcc                                          <  3:142559/135‑1 (MQ=255)                                                                                           tttCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTAAGGGGTGGCCATGCGCACCGACGATGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGgc                             >  8:425045/1‑135 (MQ=255)                                                                                                 cGAAGCCCGCAACGGCCTGGTCGACCGTGCATAACAGGTGCACATGCGCACAGACGACCGCGCGCGCAGCTGGCACGCGGAGAACGCGCGACAGCCCCAACGGCGTCGACAGCACCAGTCAGCCCTGccgcggcc                       >  8:214923/1‑132 (MQ=255)                                                                                                     ccccGCACCGGCCGGGGCGGCGGGGCCTCGGGGGGGCCCAGGCCCACCGCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGc                   <  6:142232/134‑1 (MQ=255)                                                                                                      cccGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGCCTGCCGAACGGGCGTGGCCTGCACCTGTATGCCCGGGCGCGccggtctc                  >  1:323503/1‑129 (MQ=255)                                                                                                      cccGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCa                  >  7:129095/1‑135 (MQ=255)                                                                                                       ccGCACCGGCCTGGGCGGCGGGGCCTACGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGc                   <  7:388804/133‑1 (MQ=255)                                                                                                                ccTCGTCGTCGGTGCCTACGGCGTGGCCATGCGCACCCACGACGGCGGGCGCAGCAGGCAGTCCGAAATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTCTATTCCCTGGCGCGCCACGGCAACCAGgtgt         >  5:365990/1‑134 (MQ=255)                                                                                                                   ggTCGTCGGTGCCTACGGGGTGGCCCTGCGCACCGACGCCGGCGGGCGCAGCTGGCAGGCGGAGATCGGCCGCCTGCCGACCGGGCGTGGCCTGCACCTGTATGCCCTGGCGCGCCCAGGCCACCCGGTGGACAt     >  3:68186/1‑135 (MQ=255)                                                                                                                     gcggcgGTGCCTCCGGGGGGGCCCTGCGCCCCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCg   <  4:140073/131‑1 (MQ=255)                                                                                                                      cgtcgGTGCCTACGGGGTGGCCATTCGCACCGCCGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGc  >  3:124002/1‑135 (MQ=255)                                                                                                                      cggcgGTGCCTTCGGGGGGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGc  <  4:141822/132‑1 (MQ=255)                                                                                                                            |                                                                                                                                 GGCAGCCGAGGCTGCCGCAGACGAGCGTCGCCTGGCGGCGGCGCAACGGCTGCTGGCCGACGGCGCGGACAAGCCGTTGTTGGCCCTGAATTTCACCGATGCCCGCAACGGCCTGGTCGTCGGTGCCTACGGGGTGGCCATGCGCACCGACGACGGCGGGCGCAGCTGGCAGTCGGAGATCGGCCGACTGCCGAACGGGCGTGGACTGCACCTGTATGCCCTGGCGCGCCAAGGCAACCAGGTGTACATCGC  >  ACB122/3625972‑3626223

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑