Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	3,442,493	T→A	6.8%	E346V (GAA→GTA)	copA ←	copper transporter ATPase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	3,442,493	0	T	A	6.8%	178.8 / 4.0	59	E346V (GAA→GTA)	copA	copper transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base A (2/2);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.82e-01

ATATGAGAAAGCGCAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCTGTCCGTCTCTTACAACAGTA  >  BS168_reference_genome/3442358‑3442627                                                                                                                                        |                                                                                                                                       atatGAGAAAGCGCAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGt                                                                                                                                     >  4:90687/1‑139 (MQ=255) ataaGAGAAAGCGCAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGt                                                                                                                                     >  4:44447/1‑139 (MQ=255)  tatGAGAAAGCGCAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTc                                                                                                                                    <  7:234935/139‑1 (MQ=255)  taagagaAAGCGCAGTATCCTTTCCGTCGTTGACATCTTTAATTTTCAGTAATCCGTTTTCGTTAACGGTTTAACGTGTCACGCTGTCACCAGGTTTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTc                                                                                                                                    <  1:62291/136‑1 (MQ=255)    tGAGTATGCGCAGTATCCTTTCCGACGTTGACAGCTTTAATTATGAGAAAGCCGTTAGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGAATTTTCTCAACACGAAGTCTTACTGCTGTAATCATCGATTCGTCaa                                                                                                                                  <  4:179877/139‑1 (MQ=255)      agaAAGCGCAGTATCCTTTCAGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAAc                                                                                                                                 <  4:3505/138‑1 (MQ=255)       gaAAGCGCAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCg                                                                                                                               <  7:376596/139‑1 (MQ=255)            ctcAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGaa                                                                                                                          <  3:54453/137‑1 (MQ=255)            cgcAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGaa                                                                                                                          >  4:216360/1‑139 (MQ=255)            cgcAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGaa                                                                                                                          >  2:78813/1‑139 (MQ=255)            cgcAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGaa                                                                                                                          <  8:221001/139‑1 (MQ=255)                  aTCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCt                                                                                                                    >  3:147240/1‑139 (MQ=255)                      tttCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCCAGCGGTCACCAGGAGTTTTGTCAACAGGAAGACTTTCGCCTGTAATCATCGATTCGTCAACCCGTGAAAGGCCttgcc                                                                                                                >  2:172969/1‑136 (MQ=255)                        tCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAAc                                                                                                               <  2:29562/138‑1 (MQ=255)                             cGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCTACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGaa                                                                                                                          <  6:292657/122‑1 (MQ=255)                             cGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGaa                                                                                                                          >  5:292657/1‑122 (MQ=255)                             cGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTtc                                                                                                         >  4:353463/1‑139 (MQ=255)                                 gACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCg                                                                                                     >  3:168463/1‑139 (MQ=255)                                 gACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCg                                                                                                     >  2:369688/1‑139 (MQ=255)                                         tAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAg                                                                                              >  4:111368/1‑138 (MQ=255)                                                cAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAAttcgt                                                                                      >  6:243971/1‑139 (MQ=255)                                                    aaGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAATACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGATCGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGc                                                                                  <  3:44447/139‑1 (MQ=255)                                                    aaGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGc                                                                                  <  3:90687/139‑1 (MQ=255)                                                         gTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGt                                                                             <  4:35475/139‑1 (MQ=255)                                                         gTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGt                                                                             <  4:10526/139‑1 (MQ=255)                                                         gTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGt                                                                             <  8:255648/139‑1 (MQ=255)                                                         gTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGt                                                                             <  1:369688/139‑1 (MQ=255)                                                               gTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACg                                                                       <  3:190712/139‑1 (MQ=255)                                                                     ggTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACa                                                                 <  3:353463/139‑1 (MQ=255)                                                                           aCCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCa                                                           <  7:164101/139‑1 (MQ=255)                                                                               gTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGCCAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTaa                                                       >  6:138458/1‑139 (MQ=255)                                                                                  aCGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACa                                                     <  3:73627/138‑1 (MQ=255)                                                                                    gCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGaa                                                                                                                          <  7:313911/67‑1 (MQ=255)                                                                                    gCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGaa                                                                                                                          >  8:313911/1‑67 (MQ=255)                                                                                              ggATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGa                                        <  3:216360/139‑1 (MQ=255)                                                                                                    ttGTCTACAGGAAGCCGTTCTCCTGTAATCATCGATTCGTCAACCTCTGAACGGCCTTCAAGCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTTAACAATGTCATTAACAAGCACTTCGTCGATCgtg                                   <  8:29526/138‑3 (MQ=255)                                                                                                    ttGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGga                                  >  1:271538/1‑139 (MQ=255)                                                                                                     tGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGga                                  <  7:122273/138‑1 (MQ=255)                                                                                                     tGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACGTCTCCGTCGACAGGAATTCGGTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGga                                  <  5:62324/138‑1 (MQ=255)                                                                                                         aaCAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGAtt                             >  8:277730/1‑139 (MQ=255)                                                                                                                      tctcCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCTGTCCGt                >  5:229473/1‑139 (MQ=255)                                                                                                                      tctcCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCTGTCCGt                >  5:27340/1‑139 (MQ=255)                                                                                                                          cTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCTGTCCGtc               <  3:111368/136‑1 (MQ=255)                                                                                                                          cTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGGCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCTGTCCGTctct            <  8:97115/139‑1 (MQ=255)                                                                                                                           tGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAAAAATTTCATTAACAAGCACTTCGTAGATCTTGATGATTTGCTCCAGTCCGTctct            >  6:21352/1‑138 (MQ=255)                                                                                                                             tAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCt                      <  4:346661/126‑1 (MQ=255)                                                                                                                             tAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCt                      >  3:346661/1‑126 (MQ=255)                                                                                                                              aaTCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGt                                                                             >  3:212243/1‑70 (MQ=255)                                                                                                                              aaTCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGt                                                                             <  4:212243/70‑1 (MQ=255)                                                                                                                                  atcGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGgatgat                              >  3:115867/1‑113 (MQ=255)                                                                                                                                  atcGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGgatgat                              <  4:115867/113‑1 (MQ=255)                                                                                                                                  atcGATACGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGAtt                             >  6:388887/1‑114 (MQ=255)                                                                                                                                  atcGATACGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGAtt                             <  5:388885/114‑1 (MQ=255)                                                                                                                                  atcGATACGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGAtt                             >  6:388885/1‑114 (MQ=255)                                                                                                                                  atcGATACGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGAtt                             <  5:388887/114‑1 (MQ=255)                                                                                                                                    cGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGga                                  >  2:257397/1‑107 (MQ=255)                                                                                                                                    cGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGga                                  <  1:257397/107‑1 (MQ=255)                                                                                                                                    cGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCTGTCCGTCTCTTACAACAg    >  1:397723/1‑137 (MQ=255)                                                                                                                                    cGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCTGTCCGTCTCTTACAACAGTa  >  2:95105/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                       ATATGAGAAAGCGCAGTATCCTTTCCGACGTTGACAGCTTTAATTTTCAGAAAGCCGTTTGCGTTAACGGTTGAACCTGTCACGCTGTCACCAGGATTTTTGTCAACAGGAAGACTTTCTCCTGTAATCATCGATTCGTCAACCGCTGAACGGCCTTCAACCACTTCTCCGTCGACAGGAATTCGTTCGCCAGGTTTGACGTAAACAATGTCATTAACAAGCACTTCGTCGATCGGGATGATTTGCTCCTGTCCGTCTCTTACAACAGTA  >  BS168_reference_genome/3442358‑3442627

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