Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	506,071	C→A	88.4%	P307Q (CCG→CAG)	ydbM →	putative acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	506,071	0	C	A	88.4%	110.4 / 9.9	43	P307Q (CCG→CAG)	ydbM	putative acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (20/18);  total (23/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.88e-01

AGTATTCTCCAAACAGCTTGAATGGACCGATCAAGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCCGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCAAAAGCCTGGAACGAACCAACCCGCTTCAGCG  >  BS168_reference_genome/505939‑506206                                                                                                                                     |                                                                                                                                        aGTATTCTCCAAACAGCTTGAATGGACCGATCAAGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTgcg                                                                                                                                   <  8:121224/139‑1 (MQ=255)  gTATTCACCAAACAGCTTGAAGGTACCCTTCTACATTGTCCCGGCCGACCAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTgcg                                                                                                                                   <  6:31162/138‑1 (MQ=255)     ttCTCCAAACAGCTTGAATGGACCGATCAAGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGt                                                                                                                               >  8:22341/1‑139 (MQ=255)           aaaCAGCTTGAATGGACCGATCAAGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTc                                                                                                                              <  5:2700/134‑1 (MQ=255)           aaaCAGCTTGAATGGACCGATCAAGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTc                                                                                                                              >  6:2700/1‑134 (MQ=255)                  ttGAATGGACCGATCAAGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGAc                                                                                                                    >  6:171720/1‑137 (MQ=255)                           ccGATCAAGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCg                                                                                                         >  4:90957/1‑139 (MQ=255)                              aTCAAGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTg                                                                                                      >  6:146588/1‑139 (MQ=255)                                 aaGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTgcgg                                                                                                   <  5:146588/139‑1 (MQ=255)                                    aaTGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGa                                                                                                <  7:22341/139‑1 (MQ=255)                                        tCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACa                                                                                            >  6:153579/1‑139 (MQ=255)                                        tCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTACCATATCGCTCAATTGTATCACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGAGaacc                                                                                             >  5:17355/1‑136 (MQ=255)                                                  aCAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGAGTCCACACTTGACCAGCGAGATCGGTGCGGCGAAACATATCGTcac                                                                                   >  6:192825/1‑138 (MQ=255)                                                      cAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcg                                                                              >  4:97577/1‑139 (MQ=255)                                                      cAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcg                                                                              <  3:90957/139‑1 (MQ=255)                                                           aagaGGAGATATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCTCTCAATTGTATGACGATCAGGTGCGGCGTCCACTCTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcgctgcg                                                                         <  6:58863/136‑1 (MQ=255)                                                           aaCAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcgctgcg                                                                         <  5:14824/139‑1 (MQ=255)                                                                gCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCCGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACAGATCGTCACAAACGCTGCGctgtc                                                                    >  1:90906/1‑139 (MQ=255)                                                                           cTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATTGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATa                                                         <  3:121177/139‑1 (MQ=255)                                                                                   gCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCCGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGAt                                                  >  5:87197/1‑138 (MQ=255)                                                                                    cTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCCGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACata                                                                                          <  5:223000/97‑1 (MQ=255)                                                                                    cTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCCGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACata                                                                                          >  6:223000/1‑97 (MQ=255)                                                                                    cTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGAt                                                  >  5:219406/1‑137 (MQ=255)                                                                                    cTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGAt                                                  >  7:86584/1‑137 (MQ=255)                                                                                            cAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCCGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTc                                        <  6:87197/139‑1 (MQ=255)                                                                                              gACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcgctgc                                                                          >  4:6134/1‑103 (MQ=255)                                                                                              gACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcgctgc                                                                          <  3:6134/103‑1 (MQ=255)                                                                                              gACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGAGAGCTCCGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCgg                                      >  7:149253/1‑139 (MQ=255)                                                                                              gACACTTTCTTATCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACATGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCgg                                      >  7:216281/1‑139 (MQ=255)                                                                                                   tttctttTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcgctgcg                                                                         >  8:173253/1‑99 (MQ=255)                                                                                                   tttctttTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcgctgcg                                                                         <  7:173253/99‑1 (MQ=255)                                                                                                    ttctgtTCCATATCGCTCAATTGTATGACGATCAGGTGCGGCGTCCAGACTTGACCCGGGAGCTCGGTGCGGCGAAACATTTCGTCACAAACGCTGCGCTGTCTGTTGTTGATATAGCGATGAGAATTGTCGGGGCCa                                 <  5:192825/138‑1 (MQ=255)                                                                                                     tctttTCCATATAGCTCAATTGTTTGACGTTCAGGTGCGGCGTCCACACTTGTCAAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCaa                                <  6:17355/138‑1 (MQ=255)                                                                                                                   cTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCCGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCAAAAGCCTGgaacgaac                 <  5:171720/139‑1 (MQ=255)                                                                                                                   cTCAATTGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCAAAAGCCTGgaacgaac                 >  8:12108/1‑139 (MQ=255)                                                                                                                   atCAATTGGATGACGAGCAGGTGCGGCGTCGACACTTGACCAGCGAGCTCGGTGCGTCGAAACATATGGTCACAAACGCTGAGCTGTCTGTTGTTGATAAATCGATGAGACATGTCGGGGCTAAAATCCTGgaacgaac                 <  2:185161/138‑1 (MQ=255)                                                                                                                        ttGTATGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCAAAAGCCTGGAACGaaccaacc            <  7:12108/139‑1 (MQ=255)                                                                                                                             tGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCAAAAGCCTGGAACGAACCAACCCGCtt       >  7:120160/1‑139 (MQ=255)                                                                                                                             tGACGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCAAAAGCCTGGAACGAACCAACCAGCtt       >  4:137397/1‑139 (MQ=255)                                                                                                                                cGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcgctgcg                                                                         >  7:216167/1‑70 (MQ=255)                                                                                                                                cGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAAcgctgcg                                                                         <  8:216167/70‑1 (MQ=255)                                                                                                                                cGATCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCAAAAGCCTGGAACGAACCAACCCGCTTCAg    <  6:219406/139‑1 (MQ=255)                                                                                                                                  aTCAGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCAAAAGCCTGGAACGAACCAACCCGCTTCAgcg  >  7:149616/1‑139 (MQ=255)                                                                                                                                     |                                                                                                                                        AGTATTCTCCAAACAGCTTGAATGGACCGATCAAGAATGTCCCGGCCGTACAGCAAAGAACAGGCGAAATGGAGCTTGAACTGCTCAATGCCAGACACTTTCTTTTCCATATCGCTCAATTGTATGACGATCCGGTGCGGCGTCCACACTTGACCAGCGAGCTCGGTGCGGCGAAACATATCGTCACAAACGCTGCGCTGTCTGTTGTTGATAAAGCGATGAGAATTGTCGGGGCCAAAAGCCTGGAACGAACCAACCCGCTTCAGCG  >  BS168_reference_genome/505939‑506206

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