Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	146,055	C→T	7.0%	P60L (CCT→CTT)	adk →	adenylate kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	146,055	0	C	T	7.0%	184.9 / 3.0	57	P60L (CCT→CTT)	adk	adenylate kinase
Reads supporting (aligned to +/- strand):  ref base C (40/13);  new base T (2/2);  total (42/15)
Fisher's exact test for biased strand distribution p-value = 2.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.25e-01

CGAACGAATTGTTGAGGATTATGGGATTCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTGAGGAATACGGC  >  BS168_reference_genome/145927‑146188                                                                                                                                 |                                                                                                                                      cgaacgaaTTGTTGAGGATTATGGGATTCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGt                                                                                                                              >  1:218700/1‑138 (MQ=255)     cgaaTTGTTGAGGATTATGGGATTCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACaa                                                                                                                         >  5:193942/1‑139 (MQ=255)         ttgttgAGGATTATGGGATTCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCg                                                                                                                      >  6:276762/1‑138 (MQ=255)          tgttgCGGATTATGGCATTCCTCATAACTCAACAGTATATTTGTTCGGGGCTCCTATAAAACAACAACCACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCgg                                                                                                                     <  1:54119/138‑1 (MQ=255)                      tGGGATTCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTaaaa                                                                                                           >  4:79304/1‑136 (MQ=255)                           ttCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACt                                                                                                   >  1:185117/1‑139 (MQ=255)                           ttCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACt                                                                                                   >  3:389219/1‑139 (MQ=255)                           ttCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACt                                                                                                   >  7:26794/1‑139 (MQ=255)                           ttCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACt                                                                                                   >  1:101616/1‑139 (MQ=255)                                        aaCAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAGACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGCCGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAgatgat                                                                                      >  8:127536/1‑139 (MQ=255)                                        aaCAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGGGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAgatgat                                                                                      >  3:100897/1‑139 (MQ=255)                                        aaCAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAgatgat                                                                                      >  1:322358/1‑139 (MQ=255)                                        aaCAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAgatgat                                                                                      >  5:182777/1‑139 (MQ=255)                                        aaCAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAgatgat                                                                                      >  2:195431/1‑139 (MQ=255)                                        aaCAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTGGGCAAAgatgat                                                                                      >  3:299075/1‑139 (MQ=255)                                        aaCAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACACGGACTCGAAGCAAAATCTTATATTGATAACGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGGAAAAGAGAGACTTGGCAAAgatgat                                                                                      >  2:266926/1‑139 (MQ=255)                                             gagaTATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATtgtg                                                                                  >  8:68927/1‑138 (MQ=255)                                             gagaTATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGCCAAAGATGATtgtg                                                                                  >  1:292181/1‑138 (MQ=255)                                               gaTATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGa                                                                                 >  1:30327/1‑137 (MQ=255)                                                    gTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTCCCTGACGAAGTCACAATCGGAATAGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAgga                                                                          >  2:310154/1‑138 (MQ=255)                                                     ttCCGTCCAGCAATCAAACAACACCCACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGtt                                                                         <  7:314114/139‑1 (MQ=255)                                                     tcccGGGCTGCAATCAAACAACAACCACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGtt                                                                         <  1:20092/137‑1 (MQ=255)                                                          tgctgcTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTttctt                                                                    >  6:116187/1‑139 (MQ=255)                                                          tgctgcTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGAGAAGGGAGAGCTCGTACCTGCCGAAGTCCCAATCGGAATTGTAAAAGAGAGCCTTGGCAAAGATGATTGTGAAAGAGGTTttctt                                                                    >  2:29507/1‑139 (MQ=255)                                                             tgctATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTg                                                                 >  5:36971/1‑139 (MQ=255)                                                               ctATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGa                                                               >  7:88807/1‑139 (MQ=255)                                                                    aaaGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGAt                                                          >  1:365839/1‑139 (MQ=255)                                                                         agaaACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCcg                                                     >  2:241110/1‑139 (MQ=255)                                                                                  acTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTc                                            >  1:253111/1‑139 (MQ=255)                                                                                  acTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTc                                            >  2:102227/1‑139 (MQ=255)                                                                                  acTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTAGGCAAAGATGATTGTGAATGAGGTTTTCTTCTGGACGGATTTCCGCTAACAGTc                                            >  6:32947/1‑139 (MQ=255)                                                                                      ggACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCt                                         >  2:236778/1‑138 (MQ=255)                                                                                      ggACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCt                                         <  6:36971/138‑1 (MQ=255)                                                                                      ggACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTc                                        <  8:88807/139‑1 (MQ=255)                                                                                             aaGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGa                                 <  2:218700/139‑1 (MQ=255)                                                                                                        tataTTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAg                               >  8:374282/1‑130 (MQ=255)                                                                                                        tataTTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAg                               <  7:374282/130‑1 (MQ=255)                                                                                                          tatTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTgaaga                     <  2:101616/138‑1 (MQ=255)                                                                                                               aTAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGCACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCt               >  6:149454/1‑139 (MQ=255)                                                                                                               aTAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACa                                               >  1:351352/1‑107 (MQ=255)                                                                                                               aTAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACa                                               <  2:351352/107‑1 (MQ=255)                                                                                                               aTAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCt               >  1:19291/1‑139 (MQ=255)                                                                                                               aTAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCt               >  1:348435/1‑139 (MQ=255)                                                                                                               aTAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCt               >  7:307586/1‑139 (MQ=255)                                                                                                                  cgggAGTGCTCGTACCTGTCGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCATAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTg             <  1:310154/137‑1 (MQ=255)                                                                                                                  aGGGCGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTg             <  4:100897/138‑1 (MQ=255)                                                                                                                  aGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTg             <  2:185117/138‑1 (MQ=255)                                                                                                                  aGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTg             <  5:276762/138‑1 (MQ=255)                                                                                                                         gCTCGTACTTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGaa                                                 <  3:362800/95‑1 (MQ=255)                                                                                                                         gCTCGTACTTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGaa                                                 <  3:368935/95‑1 (MQ=255)                                                                                                                         gCTCGTACTTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGaa                                                 >  4:362800/1‑95 (MQ=255)                                                                                                                         gCTCGTACTTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGaa                                                 >  4:368935/1‑95 (MQ=255)                                                                                                                         gCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGCACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTGAGGAATAc     >  5:285532/1‑139 (MQ=255)                                                                                                                          cTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTGAGGAATACg    >  5:177373/1‑139 (MQ=255)                                                                                                                              tACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGGAGCACTTGCAGAAATTCTTGAGGAATACGGc  >  2:260915/1‑137 (MQ=255)                                                                                                                              tACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTGAGGAATACGGc  >  6:87877/1‑137 (MQ=255)                                                                                                                              tACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTGAGGAATACGGc  >  7:241064/1‑137 (MQ=255)                                                                                                                                 |                                                                                                                                      CGAACGAATTGTTGAGGATTATGGGATTCCTCATATCTCAACAGGAGATATGTTCCGTGCTGCTATGAAAGAAGAAACACCACTCGGACTCGAAGCAAAATCTTATATTGATAAGGGAGAGCTCGTACCTGACGAAGTCACAATCGGAATTGTAAAAGAGAGACTTGGCAAAGATGATTGTGAAAGAGGTTTTCTTCTGGACGGATTTCCGCGAACAGTCGCTCAAGCCGAAGCTCTTGAAGAAATTCTTGAGGAATACGGC  >  BS168_reference_genome/145927‑146188

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