Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	2,298,184	C→A	9.0%	G112C (GGT→TGT)	ypjQ ←	putative phosphatidylglycerophosphatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	2,298,184	0	C	A	9.0%	142.9 / 3.5	52	G112C (GGT→TGT)	ypjQ	putative phosphatidylglycerophosphatase
Reads supporting (aligned to +/- strand):  ref base C (17/30);  new base A (3/2);  total (20/32)
Fisher's exact test for biased strand distribution p-value = 3.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

TGCAATTCTGCTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTCAACAAGGTGCTGCAGCGGTTCGGGGAGAAGTTTCTGTTCTGCAAGCTGATCGAGTGCAAGG  >  BS168_reference_genome/2298047‑2298317                                                                                                                                          |                                                                                                                                      tgaaaTTCTGCTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAAcc                                                                                                                                       >  5:78208/4‑138 (MQ=255)   cAATTCTGCTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGa                                                                                                                                     <  7:123941/138‑1 (MQ=255)    aaTTCTGCTCGCCGCTGCTGCAGCAAGATCTGCCACTATATCATCGAAAATCGTGTCAATACCGTCTGTACTTTCATCAAGTTGCTTAATAATTCCAATCTTCTCTTTAACCAAATAACCGAAATTGGTCAAACCGATc                                                                                                                                   <  5:184285/139‑1 (MQ=255)    aaTTCTGCTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATc                                                                                                                                   <  4:71919/139‑1 (MQ=255)    aaTTCTGCTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATc                                                                                                                                   <  1:161919/139‑1 (MQ=255)    aaTTCTGCTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATc                                                                                                                                   <  1:166282/139‑1 (MQ=255)       tCTGCTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGAc                                                                                                                                >  6:253868/1‑139 (MQ=255)        cTGCTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGAc                                                                                                                                <  2:142725/138‑1 (MQ=255)           cTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAACTAACCGAAAGTGGTCAACCCGATCGACCCGt                                                                                                                            >  6:104848/1‑139 (MQ=255)           cTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGATCCg                                                                                                                             >  8:34306/1‑138 (MQ=255)             cgccgcTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACAGATCGACCCGTa                                                                                                                           >  1:302934/1‑138 (MQ=255)             cgccgcTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACAGATCGACCCGTa                                                                                                                           >  1:298977/1‑138 (MQ=255)                cgcTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGCTCGACCCGTAcaca                                                                                                                       >  5:152420/1‑139 (MQ=255)                     ctgcAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTa                                                                                                                    <  4:102110/137‑1 (MQ=255)                              gagCTGCCACAATATCATCCAAAAAGGTGCGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAg                                                                                                         <  7:267426/139‑1 (MQ=255)                               agCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACGTTCATCAAGTTCCTTAAAAATTCCAATCTTCTCTTTATCCAATTTACCGAAATTGTTCAATCCGATCGACACGTACCCATTAACGATTGAAAGc                                                                                                        >  3:263464/1‑139 (MQ=255)                                    ccACAATATCATCCAAAAAGGTTTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATATCCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCggg                                                                                                     <  4:141895/137‑1 (MQ=255)                                       caATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAAt                                                                                                <  6:78208/139‑1 (MQ=255)                                              atcCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATc                                                                                          >  3:87035/1‑138 (MQ=255)                                              atcCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAAATTCATc                                                                                          >  7:189963/1‑138 (MQ=255)                                               tcCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTa                                                                                        <  8:126790/139‑1 (MQ=255)                                               tcCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTa                                                                                        <  3:105316/139‑1 (MQ=255)                                                 cAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATg                                                                                      <  4:87035/139‑1 (MQ=255)                                                 cAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATg                                                                                      <  7:242993/139‑1 (MQ=255)                                                 cAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAAATTCATCTATg                                                                                      <  8:189963/139‑1 (MQ=255)                                                 cAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAACTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATg                                                                                      <  6:138820/139‑1 (MQ=255)                                                    aaaGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCg                                                                                   >  4:218260/1‑139 (MQ=255)                                                        gtgtGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTGTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTa                                                                                 >  6:257752/1‑137 (MQ=255)                                                                aCCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACAGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCa                                                                       <  2:298977/139‑1 (MQ=255)                                                                aCCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACAGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCa                                                                       <  2:302934/139‑1 (MQ=255)                                                                   gTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGTTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCt                                                                                         >  5:291797/1‑118 (MQ=255)                                                                   gTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGa                                                                                                                                     <  7:47703/74‑1 (MQ=255)                                                                   gTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGa                                                                                                                                     >  8:47703/1‑74 (MQ=255)                                                                   gTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCt                                                                                         >  7:141244/1‑118 (MQ=255)                                                                   gTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCt                                                                                         <  8:141244/118‑1 (MQ=255)                                                                   gTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCt                                                                                         <  6:291797/118‑1 (MQ=255)                                                                   gTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATcc                                                                    >  3:312049/1‑139 (MQ=255)                                                                     cTGGACTTTCATCAAGTTCCTTAATAAATCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACAGATCGACCCGTACACATTAACGATTCAAATCGGGATAATTTCATCTATGCCGTAAAGCTGTTCATCCGt                                                                  >  5:110096/1‑139 (MQ=255)                                                                        gACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGt                                                                  <  5:253868/136‑1 (MQ=255)                                                                        gACGTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGt                                                                  <  5:201957/136‑1 (MQ=255)                                                                             tcatcaAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGATCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTcaaca                                                           <  7:34306/138‑1 (MQ=255)                                                                             tcatcaAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTcaac                                                            >  4:78007/1‑137 (MQ=255)                                                                             tcatcaAGTTCATTAATAATTACAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTcaaca                                                           <  6:110096/138‑1 (MQ=255)                                                                             tcatcaAGTACCTTAATAATTCCAATCTCCTCTTAATCCAAATAACCCAAAATGGTCAAACCGATCGACCCGTACACATTAACGAAAGAAAGCGGCATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTcaaca                                                           <  4:263464/138‑1 (MQ=255)                                                                                       ccTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTTGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGAGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTGTCAACAAGGtcatg                                                   >  5:74474/1‑132 (MQ=255)                                                                                        cTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTCAACAAGGTGCTGCAGc                                               <  4:246139/139‑1 (MQ=255)                                                                                                           tctTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTCAACAAGGTGCTGCAGCGGTTCGGGGAGAAGTttct                            >  7:214945/1‑139 (MQ=255)                                                                                                                aTCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTcaaca                                                           <  2:162280/103‑1 (MQ=255)                                                                                                                aTCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTcaaca                                                           >  1:162280/1‑103 (MQ=255)                                                                                                                   cAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTCAACAAGGTGCTGCAGCGGTTCGGGGAGAAGTTTCTGTTCTGCa                    <  6:140019/139‑1 (MQ=255)                                                                                                                           cGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTCAACAAGGTGCTGCAGCGGTTCGGGGAGAAGTTTCTGTTCTGCAAGCTGATc            <  4:312049/139‑1 (MQ=255)                                                                                                                           cGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGAGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTCAACAAGGTGCTGCAGCGGTTCGGGGAGAAGTTTCTGTTCTGCAAGCTGATc            <  6:152420/139‑1 (MQ=255)                                                                                                                                      aaaCCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTACATCTATGCCGTAAAGCGGTTCATCCGGTTCAACAAGGTGCTGCAGCGGTTCGGGGAGAAGTTTCTGTTCTGCAAGCTGATCGAGTGCAAgg  <  3:145369/138‑1 (MQ=255)                                                                                                                                          |                                                                                                                                      TGCAATTCTGCTCGCCGCTGCTGCAGCAAGAGCTGCCACAATATCATCCAAAAAGGTGTGAATACCGTCTGGACTTTCATCAAGTTCCTTAATAATTCCAATCTTCTCTTTATCCAAATAACCGAAATTGGTCAAACCGATCGACCCGTACACATTAACGATTGAAAGCGGGATAATTTCATCTATGCCGTAAAGCGGTTCATCCGTTTCAACAAGGTGCTGCAGCGGTTCGGGGAGAAGTTTCTGTTCTGCAAGCTGATCGAGTGCAAGG  >  BS168_reference_genome/2298047‑2298317

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: ‑