Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	1,561,301	C→T	8.3%	P279L (CCG→CTG)	ctaC →	cytochrome caa3 oxidase (subunit II)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	1,561,301	0	C	T	8.3%	146.0 / 3.2	48	P279L (CCG→CTG)	ctaC	cytochrome caa3 oxidase (subunit II)
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base T (2/2);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.93e-01

TGTCTGCAAAAGAATTTCAGGGCTGGACAAAAGAAATGAAGAACTATAAGTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAGGCTTGGCTGAAGGACCCTGACAGCATAAAGC  >  BS168_reference_genome/1561169‑1561438                                                                                                                                     |                                                                                                                                          tGTCTGCAAAAGAATTTCAGGGCTGGACAAAAGAAATGAAGAACTATAAGTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGa                                                                                                                                     >  7:96825/1‑139 (MQ=255) tGTCTGCAAAAGAATTTCAGGGCTGGACAAAAGAAATGAAGAACTATAAGTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGa                                                                                                                                     >  7:259618/1‑139 (MQ=255)   tCTGCAAAAGAATTTCGGGGCTGGACAAAAGAAATGAAGAACTATAAGTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGAt                                                                                                                                    <  4:102907/138‑1 (MQ=255)              aTTTCAGGGCTGGACAAAAGAAATGAAGAACTATAAGTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGa                                                                                                                         <  2:88932/138‑1 (MQ=255)                          gACAAAAGAAATGCAGAACTATAAGTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGaa                                                                                                              <  1:130369/137‑1 (MQ=255)                           aCAAAAGAAATGAAGAACTATAAGTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACgg                                                                                                           <  8:96825/139‑1 (MQ=255)                           aCAAAAGAAATGAAGAACTATAAGTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACgg                                                                                                           <  8:259618/139‑1 (MQ=255)                                                  gTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGt                                                                                    >  1:104661/1‑139 (MQ=255)                                                       cagcagAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATaa                                                                                                                                  >  2:156543/1‑88 (MQ=255)                                                       cagcagAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATaa                                                                                                                                  <  1:156543/88‑1 (MQ=255)                                                       cagcagAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAg                                                                               >  1:206974/1‑139 (MQ=255)                                                         gcagAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGa                                                                             >  4:305906/1‑139 (MQ=255)                                                          cagAGAGCGATTTGGCTAAGCTGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCAATGCGGTAGTGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGAc                                                                            <  6:43779/139‑1 (MQ=255)                                                               agCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGaa                                                                                                              >  3:26911/1‑100 (MQ=255)                                                               agCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGaa                                                                                                              <  4:26911/100‑1 (MQ=255)                                                               agCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAg                                                                       >  8:42455/1‑139 (MQ=255)                                                                   aTTTGTCTATGCAGTGTGAATAACTGTTTAAAGAGAAAAACTGCCTGAGCGGCCATGCGGTAGAGCCGAATGATAAGCGTTCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGc                                                                   <  6:203356/139‑1 (MQ=255)                                                                   aTTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGc                                                                   <  2:255971/139‑1 (MQ=255)                                                                   aTTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGc                                                                   <  2:206974/139‑1 (MQ=255)                                                                    tttGGCTAAGCAGGGTGAAGATCTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCTGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTgg                                                                    >  4:269116/1‑137 (MQ=255)                                                                    tttGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCTGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTgg                                                                    >  4:273649/1‑137 (MQ=255)                                                                         cTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGt                                                             >  8:17554/1‑139 (MQ=255)                                                                                 ggTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAgg                                                     >  7:164020/1‑139 (MQ=255)                                                                                 ggTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAgg                                                     >  1:181311/1‑139 (MQ=255)                                                                                  gTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGc                                                    >  7:139631/1‑139 (MQ=255)                                                                                   tGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGc                                                    >  7:128819/1‑138 (MQ=255)                                                                                     aagtattGTTTATAGAGAAAAACTGCCTGAGCTTCCATGCGGTAGAGCCGATTGATAAGCGTGCAGTAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGTGTGAAAGAGGCcaa                                                 <  1:34336/133‑1 (MQ=255)                                                                                      agaaCTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCcaa                                                 >  6:133464/1‑138 (MQ=255)                                                                                      agaaCTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCcaa                                                 >  5:141560/1‑138 (MQ=255)                                                                                             tttAAAGAGAAAAACTGCCTTTGCTTCCATGCGGTAGAGACGAATGATAATCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAATAGGCCAACAAAGaa                                          <  6:42415/138‑1 (MQ=255)                                                                                             tttAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGaa                                          >  8:63/1‑138 (MQ=255)                                                                                             tttAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGaa                                          >  5:42415/1‑138 (MQ=255)                                                                                             tttAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGaa                                          >  2:198908/1‑138 (MQ=255)                                                                                             tttAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGaa                                          >  6:46773/1‑138 (MQ=255)                                                                                             tttAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGAAAAAATTGGCAGGGGTAAAAAAGGCCAACAAAGaa                                          >  2:146649/1‑138 (MQ=255)                                                                                              ttAAATAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGGGGCAGGGGTGAAAGAGGCCAACAAAGaaa                                         <  7:60060/138‑1 (MQ=255)                                                                                                  agagAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAAtgtg                                    <  8:164020/139‑1 (MQ=255)                                                                                                  agagAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAAtgt                                     >  2:42497/1‑138 (MQ=255)                                                                                                      aaaaaCTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAgg                                <  7:17554/139‑1 (MQ=255)                                                                                                       aaaaCTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGCAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAGGc                               >  5:187185/1‑139 (MQ=255)                                                                                                         aaCTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAGGCtt                             <  7:146115/139‑1 (MQ=255)                                                                                                               cTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAGGCTTGGCTGa                       <  4:170577/139‑1 (MQ=255)                                                                                                                         aTGCGGTAGAGCTGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAATTGTGAAGGCTTGGCTGAAGGACCCTGa             <  3:273649/139‑1 (MQ=255)                                                                                                                         aTGCGGTAGAGCTGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAATTGTGAAGGCTTGGCTGAAGGACCCTGa             <  3:269116/139‑1 (MQ=255)                                                                                                                         aTGCGGTAGAGCCGAATGATAATCCTGCAGAAGCAGCTAGAACGGCACCTAACTGGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAGGCTTGGCTGAAGGACCCTGa             <  6:187185/139‑1 (MQ=255)                                                                                                                         aTGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAGGCTTGGCTGAAGGACCCTGa             <  8:128819/139‑1 (MQ=255)                                                                                                                         aTGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAGGCTTGGCTGAAGGACCCTGa             <  5:271460/139‑1 (MQ=255)                                                                                                                         aTGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAGGCTTGGCTGAAGGACCCTGa             <  5:46773/139‑1 (MQ=255)                                                                                                                                    ccGAATGATAAGCGTGCCGACGCCGCTAGAACGGCCCCTACCGTGGCGACCGTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGATGCCACCAAAGAACATGTGAAGGCTGGGCTGAAGGACCCTGCCAGCAtcacgc  >  2:165696/1‑134 (MQ=255)                                                                                                                                     |                                                                                                                                          TGTCTGCAAAAGAATTTCAGGGCTGGACAAAAGAAATGAAGAACTATAAGTCTACAGCAGAGAGCGATTTGGCTAAGCAGGGTGAAGAACTGTTTAAAGAGAAAAACTGCCTGAGCTGCCATGCGGTAGAGCCGAATGATAAGCGTGCAGAAGCCGCTAGAACGGCACCTAACTTGGCGACCTTCGGTGAAAGGACAAAAGTGGCAGGGGTGAAAGAGGCCAACAAAGAAAATGTGAAGGCTTGGCTGAAGGACCCTGACAGCATAAAGC  >  BS168_reference_genome/1561169‑1561438

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