Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,923,281	G→A	12.5%	N647N (AAC→AAT)	valS ←	valyl‑tRNA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,923,281	0	G	A	12.5%	179.0 / 14.5	64	N647N (AAC→AAT)	valS	valyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (23/33);  new base A (3/5);  total (26/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATG  >  Exported/2923211‑2923341                                                                       |                                                             aCAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCg                                                              >  1:2950589/1‑71 (MQ=255) aCAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCg                                                              >  1:2412348/1‑71 (MQ=255)  cAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGt                                                             <  1:494890/71‑1 (MQ=255)  cAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGt                                                             <  1:1846282/71‑1 (MQ=255)  cAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGt                                                             >  1:1802274/1‑71 (MQ=255)   aGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCg                                                              >  1:722794/1‑69 (MQ=255)    gCGTCGTTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTc                                                           <  1:1652077/71‑1 (MQ=255)    gCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTc                                                           <  1:2014190/71‑1 (MQ=255)    gCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTc                                                           <  1:3266870/71‑1 (MQ=255)    gCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTc                                                           <  1:492703/71‑1 (MQ=255)      gTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCca                                                         <  1:2976266/71‑1 (MQ=255)      gTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCca                                                         <  1:1374307/71‑1 (MQ=255)         aTTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAg                                                      >  1:1908120/1‑71 (MQ=255)          tttCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGc                                                     >  1:1802991/1‑71 (MQ=255)              gccgccGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAg                                                      <  1:591889/66‑1 (MQ=255)                cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCATCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                               <  1:3041354/71‑1 (MQ=255)                cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCt                                                    >  1:27520/1‑66 (MQ=255)                cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCt                                                    >  1:819965/1‑66 (MQ=255)                cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                               <  1:594717/71‑1 (MQ=255)                cgccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCAACAGCt                                                    >  1:2182499/1‑66 (MQ=255)                 gccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                               >  1:2374124/1‑70 (MQ=255)                  ccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                               >  1:1063362/1‑69 (MQ=255)                  ccgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                               >  1:1141252/1‑69 (MQ=255)                   cgTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTa                                               >  1:605590/1‑68 (MQ=255)                          gCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGt                                                 <  1:3396312/59‑1 (MQ=255)                          gCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGt                                                 <  1:2360525/59‑1 (MQ=255)                            cGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACgg                                   >  1:313003/1‑71 (MQ=255)                             gCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAg                                            >  1:2218269/1‑61 (MQ=255)                             gCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGt                                  <  1:876502/71‑1 (MQ=255)                             gCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGt                                  <  1:3054158/71‑1 (MQ=255)                              cAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGt                                  <  1:2744344/70‑1 (MQ=255)                                      aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         <  1:2506201/71‑1 (MQ=255)                                      aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         <  1:1476123/71‑1 (MQ=255)                                      aCCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         <  1:3081431/71‑1 (MQ=255)                                       ccTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAg                        <  1:1039188/71‑1 (MQ=255)                                         ttCGGTGTTCATCATCACGAAGCTGGTCGCATTCCACAGCTTGTTGCAGa                                           <  1:1666247/50‑1 (MQ=255)                                         ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                           <  1:1853491/50‑1 (MQ=255)                                         ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                           <  1:638743/50‑1 (MQ=255)                                         ttCGGTGTTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                           <  1:1915694/50‑1 (MQ=255)                                           cggtgtTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGaa                                          >  1:1847163/3‑49 (MQ=255)                                           cggtgtTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGAAg                                         >  1:252185/3‑50 (MQ=255)                                           cggtgtTCATCATCACGAAGCGGGTCGCATTCCACAGCTTGTTGCAGAAg                                         >  1:3149956/3‑50 (MQ=255)                                           cTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACg                     <  1:2735248/70‑1 (MQ=255)                                           cTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGc                    >  1:2921775/1‑71 (MQ=255)                                           cTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGc                    >  1:1226892/1‑71 (MQ=255)                                            tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         >  1:522003/1‑65 (MQ=255)                                            tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         >  1:433400/1‑65 (MQ=255)                                            tgtgTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCa                         >  1:2185184/1‑65 (MQ=255)                                              tgtTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGc                    <  1:1545896/68‑1 (MQ=255)                                              tgtTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCtt                  <  1:2107639/70‑1 (MQ=255)                                              tgtTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCtt                  <  1:3341670/70‑1 (MQ=255)                                              tgtTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCt                   <  1:3015916/69‑1 (MQ=255)                                                      aGCACAAAGCGGCTGGCGTTACACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATcc           >  1:89855/1‑69 (MQ=255)                                                        cacgAAGCGGGTCGCATTCCACAGCTTGTTGCAGa                                           <  1:3458943/35‑1 (MQ=16)                                                         acaAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGtt      <  1:1641966/71‑1 (MQ=255)                                                         acaAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGtt      <  1:2736336/71‑1 (MQ=255)                                                          caAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTa                                 <  1:217556/43‑1 (MQ=255)                                                           aaaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGa    >  1:2555996/1‑71 (MQ=255)                                                            aaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCatat             >  1:3187046/1‑61 (MQ=255)                                                            aaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGa    >  1:1464942/1‑70 (MQ=255)                                                            aaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGa    >  1:2678563/1‑70 (MQ=255)                                                            aaGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGAt   >  1:1549479/1‑71 (MQ=255)                                                             aGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATg  <  1:3418372/71‑1 (MQ=255)                                                                      cGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTg     <  1:2274417/59‑1 (MQ=255)                                                                      cGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTg     <  1:2355920/59‑1 (MQ=255)                                                                      cGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTg     <  1:1519350/59‑1 (MQ=255)                                                                      cGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTg     <  1:1496548/59‑1 (MQ=255)                                                                       |                                                             ACAGCGTCATTTCGCCGCCGTTGAAGCCGCAATCCTGACCTTCTGTGTTCATCAGCACAAAGCGGCTGGCGTTCCACAGCTTGTTACAGAAGTTACGGTAACCTTCCAGACGCTTCATATCCCAGTTGATG  >  Exported/2923211‑2923341

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