Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	BS168_reference_genome	4,155,001	T→G	21.9%	intergenic (‑30/+52)	trnY‑Asp ← / ← trnY‑Glu	tRNA‑Asp/tRNA‑Glu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	BS168_reference_genome	4,155,001	0	T	G	21.9%	66.3 / 5.2	43	intergenic (‑30/+52)	trnY‑Asp/trnY‑Glu	tRNA‑Asp/tRNA‑Glu
Reads supporting (aligned to +/- strand):  ref base T (12/21);  new base G (8/2);  total (20/23)
Fisher's exact test for biased strand distribution p-value = 2.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGTTTTATTTTTCTCTACTAAGTAAAGTAAAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCGGTGTCTTAACCGCTTGACCAACGGGCCTTCACGAAATGGT  >  BS168_reference_genome/4154863‑4155137                                                                                                                                           |                                                                                                                                         tGTTTTATTTTTCTCTACTAAGTAAAGTAAAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCgtgt                                                                                                                                          <  3:422896/139‑1 (MQ=255) tGTTTTATTTTTCTCTACTAAGTAAAGTAAAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACACGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCgtgt                                                                                                                                          >  7:126143/1‑139 (MQ=255)     ttATTTTTCTCTAATAAGTAAAGTAAAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAg                                                                                                                                       <  6:29705/138‑1 (MQ=255)           ttCTCTACTAAGTAAAGTAAAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTca                                                                                                                                 <  4:8711/138‑1 (MQ=255)             ctctACTAAGTAAAGTAAAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTcaca                                                                                                                               <  7:130860/138‑1 (MQ=255)               ctagtaATTAAAGTAAATTGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGGGACAGGCAGGCATTCTAACCAACTGAACTACCGCACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGa                                                                                                                             <  8:368440/138‑1 (MQ=21)                  ctaAGTAAAGTAAAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGAc                                                                                                                            <  3:93523/136‑1 (MQ=255)                    aagtaaagtaAAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATa                                                                                                                       <  5:287172/139‑1 (MQ=255)                     agtaaagtaaAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGCACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGGCCGCTCTCACAGACAAAt                                                                                                                        >  8:94612/1‑137 (MQ=255)                      gtaaagtaacAGGGCGGGCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAg                                                                                                                     <  6:102505/139‑1 (MQ=255)                            taaatgggCGGGCCGGCCGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCTTTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGGCGGGGCAGCTCTCACAGACAAATAAGGAAttt                                                                                                               <  4:338551/133‑1 (MQ=21)                            taaaATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCCGCTCTCACAGACAAATAAGGAAtt                                                                                                                >  7:281777/1‑138 (MQ=255)                                 tGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTGAGCGGTCACAGGCAAAGAAGGTATTTAGCa                                                                                                           >  3:79012/1‑138 (MQ=21)                                 tGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACACACCAATAAGGAATTTATCa                                                                                                           >  4:125457/1‑138 (MQ=255)                                 tGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCACCGCTCACCCACAAATAATGAATTTAGCa                                                                                                           >  2:50336/1‑138 (MQ=21)                                  ggCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAAt                                                                                                         <  5:73691/139‑1 (MQ=255)                                      gTCCGGACGGGACTAGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCGAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCCGCTCTCAAAGACAAATAAGGAATTTAGCAATAACg                                                                                                     >  7:234179/1‑139 (MQ=255)                                      gTCCGGACCGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTGAACTCTCACACACAAAGAAAGTATTTAGGAAAAAca                                                                                                     >  6:403840/1‑138 (MQ=18)                                              gggACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTa                                                                                             <  2:82846/139‑1 (MQ=255)                                                 aCTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTGGGGACGGCCCGACAAATAAAGAATATATCAAAAACCCTAAAGTaaa                                                                                           >  1:389923/1‑138 (MQ=18)                                                     gAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAAt                                                                                      <  1:175833/139‑1 (MQ=255)                                                     gAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAAt                                                                                      <  6:184144/139‑1 (MQ=255)                                                     gAACACGCGACCTCCTGCGTGTCAGGCAGGCATTGTAACCAACTGTAATAGCGGTCCTTACATTAATTGAGCGGAGATCTGTCGTGTCAGGTCTAACAGACAAATAAGGAATTTAGCAATAAGGGTAAAGTAAATCAAt                                                                                      <  7:314229/139‑1 (MQ=255)                                                          cgcgACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAAt                                                                                                         >  2:9829/1‑115 (MQ=255)                                                          cgcgACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAAt                                                                                                         <  1:9829/115‑1 (MQ=255)                                                          cgcgACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGGCTGGGGGGCCACACAAATAAAGAATTTATAAATAAACCGAAATAAAAAAAATGAcca                                                                                 >  2:174192/1‑138 (MQ=12)                                                          cgcgACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGGCAGCGCTCCCAGACAAAAAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGAccc                                                                                 >  7:97583/1‑139 (MQ=255)                                                               cctcctGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGGCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACg                                                                            <  3:276470/139‑1 (MQ=255)                                                                ctcctGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACgg                                                                           <  3:125457/139‑1 (MQ=255)                                                                     gCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGGTCTCACAGACAAAAAATGAATTTAGCAATAACGCTAAATTAAATCAATGACCCGTAAggggt                                                                        >  4:357632/1‑135 (MQ=38)                                                                           caggcaggCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTGTGCAGTCACAGACCAATAAAGAATTTATAAAAAACACTAAAGAAAATAAAAAACCCGAAAGGGAGTAGGAccc                                                                >  3:10325/1‑139 (MQ=19)                                                                              gcaggcaTTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGGCAGCGCGCACACACAAATAAGTAATTTAGCAATAAAGCTAAAGTAAATAAATGAACCCTACGGGGTGCGAACCCggg                                                             >  1:88687/1‑137 (MQ=255)                                                                                aggcaTTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGGGAGCTCGCACAGACAAATAAAGAATTTTTGAATAACCCAAAAATAAATCAAATACCCCTAAGGGTGTCGTACCCGGGtt                                                           >  6:387636/1‑139 (MQ=255)                                                                                          aCCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGGCAGCGCTCACAGACAAATAAGGAATTTATCAATAAAGCTAAAGTAAATAAATGACCCGTACGGGCTTCGAACCCGGGTTAACGCCGTg                                                  >  6:107981/1‑138 (MQ=255)                                                                                                tGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGCCAAATAAGGAATTTATCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCCTGTTTACGCCGGGGAAGGGc                                           >  4:112739/1‑139 (MQ=255)                                                                                                tGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCGCGCACACACAAAAAAGGAATTTAGCAATAACGTTAAAATAAATAAATGACCCGTACGGGATTCGCACCCCTGTTAACGCCGTGGAAGGGc                                           >  4:345179/1‑139 (MQ=37)                                                                                                  caCTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGc                                           <  8:218564/136‑1 (MQ=255)                                                                                                  aaCTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGc                                           <  1:174192/137‑1 (MQ=255)                                                                                                  aaCTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGc                                           <  7:3249/137‑1 (MQ=255)                                                                                                  aaCTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCgg                                         <  8:64782/139‑1 (MQ=255)                                                                                                  aaCTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGGCCGCGGGCACAGACAAATAAAGAATTTATAAATAACGCTAAAGTAAATAAAAGACCCCTACGGGAGTCGCACCCGTGTTAACCCCGGGGAAGGGCgg                                         >  3:242841/1‑139 (MQ=255)                                                                                                  aaCTACCGGACCTTACAGTAATTGCCCGGAGAGCTGGCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCg                                          <  3:112739/138‑1 (MQ=255)                                                                                                                 cAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCGGTGTCTTAACCGCTTg                          <  2:2247/139‑1 (MQ=255)                                                                                                                             ggAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCGGTGTCTTAACCGCTTGACCAACGGGCCt              <  3:75848/139‑1 (MQ=255)                                                                                                                              gagaGCTGTCGTGGCCGCGCTCACACACAAAAAACGAATTTAGCAATAACGCGAAAGTAAATCAATGAACCGTACGGGGTTCGCAACCGTGTTACCGCCGTGAAAGGGCCGGGGGTTTTCCGCATTGCCGACGGGCCtt             >  2:312973/1‑139 (MQ=16)                                                                                                                                         tgtCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACTCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCGGTGTCTTAACCGCTTGACCAACGGGCCTTCACGAAATGGt  <  6:286480/139‑1 (MQ=255)                                                                                                                                         tgtCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCGGTGTCTTAACCGCTTGACCAACGGGCCTTCACGAAATGGt  <  4:220928/139‑1 (MQ=255)                                                                                                                                         tgtCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCGGTGTCTTAACCGCTTGACCAACGGGCCTTCACGAAATGGt  <  3:357632/139‑1 (MQ=255)                                                                                                                                         tgtCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCGGTGTCTTAACCGCTTGACCAACGGGCCTTCACGAAATGGt  <  1:64296/139‑1 (MQ=255)                                                                                                                                         tgtCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCGGTGTCTTAACCGCTTGACCAACGGGCCTTCACGAAATGGt  <  5:327743/139‑1 (MQ=255)                                                                                                                                           |                                                                                                                                         TGTTTTATTTTTCTCTACTAAGTAAAGTAAAATGGCGGTCCGGACGGGACTCGAACCCGCGACCTCCTGCGTGACAGGCAGGCATTCTAACCAACTGAACTACCGGACCTTACAGTAATTGCCCGGAGAGCTGTCGTGTCAGCTCTCACAGACAAATAAGGAATTTAGCAATAACGCTAAAGTAAATCAATGACCCGTACGGGATTCGAACCCGTGTTACCGCCGTGAAAGGGCGGTGTCTTAACCGCTTGACCAACGGGCCTTCACGAAATGGT  >  BS168_reference_genome/4154863‑4155137

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