Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	ACB143	4,729,213	G→C	10.8%	A252A (GCC→GCG)	sdhA_2 ←	Succinate dehydrogenase flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	ACB143	4,729,213	0	G	C	10.8%	144.4 / 8.1	45	A252A (GCC→GCG)	sdhA_2	Succinate dehydrogenase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (18/22);  new base C (2/3);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

ACGGGCAACCACGTCGCGGCCGGCCAGGTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGACGCGTAGATACGGCCCG  >  ACB143/4729087‑4729346                                                                                                                               |                                                                                                                                      aCGGGCAACCACGTCGCGGCCGGCCAGGTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTAAAccggcgccggc                                                                                                                               >  1:85733/1‑135 (MQ=255)     gCAACCACGTCGCGGCCGGCCAGGTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATAc                                                                                                                          <  6:207160/136‑1 (MQ=255)       aaCCACGTCGCGGCCGGCCAGGTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGCCCAGTACCCCGGCGCCGGCGATACCg                                                                                                                        >  8:119274/1‑136 (MQ=255)       aaCCACGTCGCGGCCGGCCAGGTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCCGCGCCCGCGATACCg                                                                                                                        >  7:13730/1‑136 (MQ=255)                gcgGCCGGCCAGGTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTg                                                                                                                <  8:49542/135‑1 (MQ=255)                     cggccAGGTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTg                                                                                                          <  3:94189/136‑1 (MQ=255)                       gccAGGTCTTTGGCGTTCGGCGCATAACGCTCCATAAAACGTTCGCCATGTTTGTTCAGCAGATAACCGCCTTCACCACGGCAACCTTCGGTGACCAGTACGCCCGCACCGGCAATGCCGGTCGGGTGGAACTGcc                                                                                                        <  2:252019/136‑1 (MQ=255)                          aGGTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACAGCcac                                                                                                      >  1:132677/1‑135 (MQ=255)                            gTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACAtt                                                                                                   <  6:78306/136‑1 (MQ=255)                            gTCTTTGGCGTTCGGCGCATAACGCTCCATAAAACGTTCGCCATGTTTGTTCAGCAGATAACCGCCTTCACCACGGCAACCTTCGGTGACCAGTACGCCCGCACCGGCAATGCCGGTCGGGTGGAACTGCCACAtt                                                                                                   <  4:251024/136‑1 (MQ=255)                                    cGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGACCTGCCACATTTCGATGTc                                                                                           >  5:172499/1‑136 (MQ=255)                                         ggCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCa                                                                                      >  5:87837/1‑136 (MQ=255)                                          gcgcGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCAc                                                                                     <  5:108053/136‑1 (MQ=255)                                                      ccATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGc                                                                         <  1:147162/136‑1 (MQ=255)                                                           aaGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGAAGTCCTGCACCGGCCCGCCTGCACGcc                                                                    >  2:297989/1‑135 (MQ=255)                                                            aGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCTGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAg                                                                   <  3:206428/136‑1 (MQ=255)                                                               gcTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGc                                                                >  1:110345/1‑136 (MQ=255)                                                                    ccGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCCGGGCCATGc                                                           >  5:57666/1‑136 (MQ=255)                                                                    ccGTGGGCGGTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTACTGAACCGGCACGCCTGCTCGCAGGGCCATGc                                                           >  7:278675/1‑136 (MQ=255)                                                                             ttcagcAGATAACCGCCTTCACCACGGCAACCTTCGGTGACCAGTACGCCCGCACCGGCAATGCCGGTCGGGTGGAACTGCCACATTTCCATATCCTGCACCGGTACGCCGGCACGGATAGCCATGCCGACACCGt                                                  <  4:189441/131‑1 (MQ=255)                                                                             ttGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGt                                                  <  5:48720/136‑1 (MQ=255)                                                                             ttGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGt                                                  <  1:89447/136‑1 (MQ=255)                                                                             ttGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGt                                                  <  7:112676/136‑1 (MQ=255)                                                                                   aGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCgg                                            <  7:61730/136‑1 (MQ=255)                                                                                        gCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGtt                                        <  7:185332/135‑1 (MQ=255)                                                                                        gCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGtt                                        >  5:487/1‑135 (MQ=255)                                                                                          cACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGAt                                     <  2:255447/136‑1 (MQ=255)                                                                                             cTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAg                                  >  2:171623/1‑136 (MQ=255)                                                                                               tCACCACGGCAACCTTCGGTGACCAGTACGCCCGCACCGGCAATGCCGGTCGGGTGGAACTGCCACATTTCCATATCCTGCACCGGTACGCCGGAACGGATAGCCATGCCGCCACCGTCGCCGGTGTTTAtgtggg                                >  6:26936/1‑131 (MQ=255)                                                                                                 aCCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGc                               <  8:145011/135‑1 (MQ=255)                                                                                                  ccGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGCACTGCCACATTTCGATGTCCTGCACCGGCCCGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGt                             >  6:9354/1‑136 (MQ=255)                                                                                                  ccGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCg                                                                                   >  5:170025/1‑82 (MQ=255)                                                                                                  ccGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCg                                                                                   <  6:170025/82‑1 (MQ=255)                                                                                                  ccGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCCCGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGt                             >  1:194478/1‑136 (MQ=255)                                                                                                  ccGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGt                             >  1:52213/1‑136 (MQ=255)                                                                                                       gCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTtggtg                        <  1:250351/136‑1 (MQ=255)                                                                                                       gCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTtggtg                        <  5:65841/136‑1 (MQ=255)                                                                                                         aaCCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGAt                                     <  7:67277/121‑1 (MQ=255)                                                                                                         aaCCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGAt                                     >  8:67277/1‑121 (MQ=255)                                                                                                         aaCCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTtggtggt                      <  5:217558/136‑1 (MQ=255)                                                                                                         aaCCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTtggtggt                      <  7:248908/136‑1 (MQ=255)                                                                                                             cTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGAc                  <  1:262785/136‑1 (MQ=255)                                                                                                                        gTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGACGCGTAGATACg       >  8:92079/1‑136 (MQ=255)                                                                                                                         tACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGACGCGTAGATACgg      >  3:46013/1‑136 (MQ=255)                                                                                                                         tACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGACGCGTAGATACgg      >  3:16478/1‑136 (MQ=255)                                                                                                                           cacCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGACGCGTAGATACGGcc    <  1:139547/136‑1 (MQ=255)                                                                                                                             ccggcgccggcgATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGACGCGTAGATACGGCCCg  <  3:134428/136‑1 (MQ=255)                                                                                                                             ccggcgccggcgATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGACGCGTAGATACGGCCCg  <  1:197934/136‑1 (MQ=255)                                                                                                                             ccggcgccggcgATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGACGCGTAGATACGGCCCg  <  1:127189/136‑1 (MQ=255)                                                                                                                               |                                                                                                                                      ACGGGCAACCACGTCGCGGCCGGCCAGGTCTTTGGCGTTCGGCGCGTAACGCTCCATGAAGCGCTCGCCGTGGGCGTTGATCAGGTAGCCACCTTCACCGCGGCAACCCTCGGTGACCAGTACACCGGCGCCGGCGATACCGGTCGGGTGGAACTGCCACATTTCGATGTCCTGCACCGGCACGCCTGCACGCAGGGCCATGCCGATACCGTCACCGGTGTTGATCAGGGCGTTGGTGGTGGACGCGTAGATACGGCCCG  >  ACB143/4729087‑4729346

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