Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	Exported	2,415,850	G→A	22.2%	noncoding (29/77 nt)	proM ←	tRNA‑Pro

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	Exported	2,415,850	0	G	A	22.2%	117.6 / 24.1	54	noncoding (29/77 nt)	proM	tRNA‑Pro
Reads supporting (aligned to +/- strand):  ref base G (30/12);  new base A (6/6);  total (36/18)
Fisher's exact test for biased strand distribution p-value = 1.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

GAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCG  >  Exported/2415785‑2415920                                                                  |                                                                       gAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                                                                   >  1:3122587/1‑71 (MQ=255) gAAGCGGGGTTAAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                                                                   >  1:3793516/1‑71 (MQ=255)    gCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgcg                                                                   >  1:577056/1‑68 (MQ=255)        ttttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCag                                                            >  1:3108423/3‑70 (MQ=12)          ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCag                                                            <  1:2367851/69‑2 (MQ=12)          ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                          >  1:5014009/1‑71 (MQ=12)          ttCGAATTTGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGc                                                          >  1:468949/1‑71 (MQ=12)          ttCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTgc                                                                    <  1:4029939/61‑1 (MQ=255)            cAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                                >  1:3907592/1‑63 (MQ=255)            cAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTg                                                        >  1:4132185/1‑71 (MQ=255)                tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                                    <  1:1308051/70‑2 (MQ=16)                tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                                    <  1:4180925/70‑2 (MQ=16)                tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                                      >  1:2772012/2‑69 (MQ=16)                tttGGTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcg                                                      >  1:1922640/2‑69 (MQ=16)                aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                                >  1:4671372/1‑59 (MQ=255)                aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                                >  1:1358257/1‑59 (MQ=255)                aTTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTa                                                                >  1:2281846/1‑59 (MQ=255)                   ggTCGGTGATAGAGGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAc                                                               <  1:3823222/57‑1 (MQ=12)                       ggCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc                                             <  1:1147009/71‑1 (MQ=255)                        gCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGcc                                             >  1:1599202/1‑70 (MQ=255)                          gatagagGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTGCGCCAATCACCGaa                                          >  1:1185964/4‑71 (MQ=12)                          gagagaGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGa                                           >  1:3291001/1‑70 (MQ=255)                           atagagGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTAc                                                               <  1:3747175/47‑1 (MQ=9)                              gagGATTCGAACCTCCGACCCCTTCGTCCCGAACGAAGTGCGCTACCAGGCTgcgcc                                                    <  1:4244743/57‑2 (MQ=12)                                     cGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCtttt                               >  1:2900712/1‑71 (MQ=255)                                            ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCt                                                    >  1:561022/1‑43 (MQ=255)                                            ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTa                                                   >  1:2476629/1‑44 (MQ=255)                                            ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCg                                            <  1:3706236/51‑1 (MQ=255)                                            ccGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                        >  1:2233676/1‑71 (MQ=255)                                             cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                       <  1:1273409/71‑1 (MQ=255)                                             cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                       <  1:4269871/71‑1 (MQ=255)                                             cGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                       <  1:2270917/71‑1 (MQ=255)                                              gACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGaa                                          >  1:5520604/1‑51 (MQ=255)                                              gACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAtttt                       >  1:291555/1‑70 (MQ=255)                                                   aCTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATAc                                      >  1:3173002/1‑50 (MQ=255)                                                    cTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGaaa                                         <  1:4971219/46‑1 (MQ=255)                                                    cTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGaaa                                         <  1:2936474/46‑1 (MQ=255)                                                    cTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTc                 >  1:1032025/1‑70 (MQ=255)                                                      ggTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATa                                       >  1:1049354/1‑46 (MQ=255)                                                      ggTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAt                          >  1:2381708/1‑59 (MQ=255)                                                       gTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCaa               <  1:3630416/69‑1 (MQ=255)                                                          ccAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCtt          >  1:3243/1‑71 (MQ=255)                                                          ccAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCt           <  1:2149250/70‑1 (MQ=255)                                                          ccAAACCAGTTGCGCTACCAAGCTGAGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCt           <  1:2672216/70‑1 (MQ=255)                                                            aaaCCAGTTGCGCTACGAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt         >  1:2451534/1‑70 (MQ=255)                                                            aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt         >  1:2421618/1‑70 (MQ=255)                                                            aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCttt         >  1:5201225/1‑70 (MQ=255)                                                            aaaCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCtt          >  1:2426201/1‑69 (MQ=255)                                                              aCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGc                                   >  1:3526204/1‑42 (MQ=255)                                                              aCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTa                     >  1:270388/1‑56 (MQ=255)                                                              aCCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTaaa      <  1:1735515/71‑1 (MQ=255)                                                                 aGTTTCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                        >  1:5364138/1‑50 (MQ=255)                                                                 aGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                        >  1:255122/1‑50 (MQ=255)                                                                 aGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAAttt                        >  1:5595757/1‑50 (MQ=255)                                                                  gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg  <  1:1975606/71‑1 (MQ=255)                                                                  gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg  <  1:289974/71‑1 (MQ=255)                                                                  gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg  <  1:5412913/71‑1 (MQ=255)                                                                  gTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCg  <  1:1000378/71‑1 (MQ=255)                                                                  |                                                                       GAAGCGGGGTTCAAAATTGGTCGGCGAGAGAGGATTCGAACCTCCGACCCACTGGTCCCAAACCAGTTGCGCTACCAAGCTGCGCTACTCGCCGAAATACTGCTTTTTGAATTTTTAGTTCAATTCTTTAAAGTCG  >  Exported/2415785‑2415920

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