Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	1,246,787	C→T	59.1%	G70G (GGC→GGT)	znuC →	zinc transporter subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	1,246,787	0	C	T	59.1%	34.7 / 79.3	66	G70G (GGC→GGT)	znuC	zinc transporter subunit
Reads supporting (aligned to +/- strand):  ref base C (16/11);  new base T (16/23);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 2.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.28e-01

TACGGGTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTT  >  Exported/1246717‑1246852                                                                       |                                                                  tACGGGTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGc                                                                   <  1:1469516/71‑1 (MQ=255)  aCGGGTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGc                                                                   >  1:229549/1‑70 (MQ=255)   cGGGTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTa                                                                 >  1:842117/1‑71 (MQ=255)    gggTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGCTAt                                                                >  1:2421605/1‑71 (MQ=255)     ggTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATg                                                               <  1:3397748/71‑1 (MQ=255)     ggTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGCTATg                                                               <  1:964335/71‑1 (MQ=255)     ggTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGCTATg                                                               <  1:592642/71‑1 (MQ=255)      gTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGt                                                              <  1:394029/71‑1 (MQ=255)      gTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGt                                                              <  1:389566/71‑1 (MQ=255)      gTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGt                                                              <  1:2213201/71‑1 (MQ=255)      gTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGt                                                              <  1:2104519/71‑1 (MQ=255)            cTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCg                                                          >  1:1773825/1‑69 (MQ=255)                  cTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCt                                                  <  1:3330847/71‑1 (MQ=255)                    ggTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTg                                                 >  1:613794/1‑70 (MQ=255)                    ggTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTg                                                 >  1:1714440/1‑70 (MQ=255)                    ggTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTg                                                 >  1:2851991/1‑70 (MQ=255)                     gTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCTGTa                                               >  1:1205692/1‑71 (MQ=255)                     gTAACACCCGATCAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGt                                                <  1:1776191/70‑1 (MQ=255)                         cacCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTAtctc                                           <  1:3301369/71‑1 (MQ=255)                          acCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCTGTAtctc                                           >  1:3183471/1‑70 (MQ=255)                           cccGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGa                                         <  1:1069112/71‑1 (MQ=255)                            ccGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTAt                                              >  1:1064411/1‑65 (MQ=255)                             cGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGaca                                       >  1:187879/1‑71 (MQ=255)                              gATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGacac                                      >  1:307972/1‑71 (MQ=255)                                 gAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGAcaccac                                   >  1:1359465/1‑71 (MQ=255)                                  aaGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACg                                  <  1:2423336/71‑1 (MQ=255)                                           tCAAGCGCAGCGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACt                          <  1:1347061/70‑1 (MQ=255)                                           tCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGt                                 <  1:484814/63‑1 (MQ=255)                                           tCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGt                                 <  1:216463/63‑1 (MQ=255)                                           tCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGt                                 <  1:2726531/63‑1 (MQ=255)                                           tCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACt                          <  1:3482309/70‑1 (MQ=255)                                           tCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACt                          <  1:536669/70‑1 (MQ=255)                                           tCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACt                          <  1:1988327/70‑1 (MQ=255)                                           tCAAGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTg                         <  1:2454261/71‑1 (MQ=255)                                             aaGCGCAACGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGa                        <  1:3425672/70‑1 (MQ=255)                                              aGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGAcc                      >  1:2405257/1‑71 (MQ=255)                                              aGCGCAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGAc                       >  1:1679221/1‑70 (MQ=255)                                                cgcAACGGAAAACTGCTCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGt                    <  1:3426700/71‑1 (MQ=255)                                                 gcAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGAcc                      <  1:44488/68‑1 (MQ=255)                                                 gcAACGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGAcc                      <  1:123984/68‑1 (MQ=255)                                                 gcAACGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGt                    >  1:405703/1‑70 (MQ=255)                                                 gcAACGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTa                   <  1:3534100/71‑1 (MQ=255)                                                     cGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAAcc               >  1:3525787/1‑71 (MQ=255)                                                     cGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAAcc               >  1:715521/1‑71 (MQ=255)                                                     cGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAAcc               >  1:2069950/1‑71 (MQ=255)                                                     cGGAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCAAGTTGCCACTGACCGTAAAcc               >  1:3550296/1‑71 (MQ=255)                                                      ggAAAACTGCGCATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCg              >  1:754314/1‑71 (MQ=255)                                                         aaaCTGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGacac                                      <  1:2965837/44‑1 (MQ=255)                                                         aaaCTGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGacac                                      <  1:950605/44‑1 (MQ=255)                                                             tGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTtta       >  1:1912939/1‑71 (MQ=255)                                                                 cATCTGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTttacgct   <  1:1259374/71‑1 (MQ=255)                                                                 cATCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTttacgct   <  1:3593601/71‑1 (MQ=255)                                                                 cATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTttacgc    <  1:1501075/70‑1 (MQ=255)                                                                  aTCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTttacgct   <  1:2502800/70‑1 (MQ=255)                                                                  aTCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTaa                  >  1:1276348/1‑55 (MQ=255)                                                                  aTCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTaa                  >  1:2086450/1‑55 (MQ=255)                                                                  aTCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTaa                  >  1:3643215/1‑55 (MQ=255)                                                                  aTCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTaa                  >  1:267963/1‑55 (MQ=255)                                                                  aTCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTttacgctt  <  1:1134068/71‑1 (MQ=255)                                                                   tCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGcc                             <  1:721706/43‑1 (MQ=255)                                                                   tCGGTTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGcc                             <  1:1781777/43‑1 (MQ=255)                                                                     ggCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTttacgct   >  1:383020/1‑67 (MQ=255)                                                                     ggCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTttacgct   >  1:2453742/1‑67 (MQ=255)                                                                     ggCTATGTACAGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTttacgct   >  1:91584/1‑67 (MQ=255)                                                                      gCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTaaa                 >  1:127323/1‑52 (MQ=255)                                                                      gCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTaaa                 >  1:1881632/1‑52 (MQ=255)                                                                      gCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTaaa                 >  1:1848215/1‑52 (MQ=255)                                                                      gCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTaaa                 >  1:2948777/1‑52 (MQ=255)                                                                       cTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGAc                       >  1:940333/1‑45 (MQ=255)                                                                       |                                                                  TACGGGTAGTGCTCGGGCTGGTAACACCCGATGAAGGGGTTATCAAGCGCAACGGAAAACTGCGCATCGGCTATGTACCGCAGAAGCTGTATCTCGACACCACGTTGCCACTGACCGTAAACCGTTTTTTACGCTT  >  Exported/1246717‑1246852

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