Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	REL606	4,479,649	C→T	noncoding (17/85 nt)	leuX →	tRNA‑Leu

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	REL606	4,479,649	0	C	T	98.5%	261.3 / ‑3.6	67	noncoding (17/85 nt)	leuX	tRNA‑Leu
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (38/28);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.93e-01

TGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAA  >  REL606/4479577‑4479720                                                                         |                                                                        tGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGt                                                                      >  1:1360173/1‑76 (MQ=255) tGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGt                                                                      >  1:258592/1‑76 (MQ=255)   cATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAg                                                                    <  1:366973/76‑1 (MQ=255)    aTGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGa                                                                   >  2:1321004/1‑76 (MQ=255)     tggtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGAc                                                                  >  2:868028/1‑76 (MQ=255)     tggtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGAc                                                                  <  1:1515623/76‑1 (MQ=255)       gtggCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGc                                                                >  1:710909/1‑76 (MQ=255)           cGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGtt                                                            >  2:1519286/1‑76 (MQ=255)           cGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGtt                                                            >  1:790427/1‑76 (MQ=255)             tGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGa                                                          >  2:408839/1‑76 (MQ=255)               cGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGAtt                                                        >  2:780427/1‑76 (MQ=255)                  cAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCaa                                                     <  2:1384081/76‑1 (MQ=255)                   aGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCaaa                                                    >  2:1380885/1‑76 (MQ=255)                    ggTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCaaaa                                                   >  1:387838/1‑76 (MQ=255)                       ataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCa                                                >  2:551442/1‑76 (MQ=255)                       ataATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCa                                                <  1:1608618/76‑1 (MQ=255)                        taATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCaa                                               >  2:457249/1‑76 (MQ=255)                         aaTCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAAc                                              >  2:1111857/1‑76 (MQ=255)                           tCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCg                                            >  1:1033078/1‑76 (MQ=255)                            ccACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGt                                           <  2:1443642/76‑1 (MQ=255)                              acaACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAg                                         <  1:385519/76‑1 (MQ=255)                               caACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGa                                        <  2:870367/76‑1 (MQ=255)                                aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGaa                                       <  1:979748/76‑1 (MQ=255)                                aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGaa                                       <  1:1555176/76‑1 (MQ=255)                                aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGaa                                       >  2:815720/1‑76 (MQ=255)                                aaCGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGaa                                       >  2:1243964/1‑76 (MQ=255)                                   gTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATa                                    <  2:282707/76‑1 (MQ=255)                                   gTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATa                                    <  1:294139/76‑1 (MQ=255)                                      ttCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGt                                 >  1:183759/1‑76 (MQ=255)                                       tCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTg                                <  2:1632699/76‑1 (MQ=255)                                       tCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTg                                >  2:1532091/1‑76 (MQ=255)                                       tCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTg                                <  2:998868/76‑1 (MQ=255)                                           cATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCgg                            <  2:1145608/76‑1 (MQ=255)                                           cATACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCgg                            >  1:1674571/1‑76 (MQ=255)                                             tACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGtt                          <  1:130333/76‑1 (MQ=255)                                             tACCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGtt                          >  1:95680/1‑76 (MQ=255)                                              aCCTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTc                         >  1:1574990/1‑76 (MQ=255)                                               ccTCTTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCg                        <  1:78507/76‑1 (MQ=255)                                                ctctTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGa                       >  1:974387/1‑76 (MQ=255)                                                 tctTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAg                      >  2:1357732/1‑76 (MQ=255)                                                  ctTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGt                     >  2:1426034/1‑76 (MQ=255)                                                  ctTCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGt                     >  2:1426048/1‑76 (MQ=255)                                                    tCAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTcc                   <  2:90122/76‑1 (MQ=255)                                                     cAGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCg                  <  2:164364/76‑1 (MQ=255)                                                      aGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCgg                 >  1:1348652/1‑76 (MQ=255)                                                      aGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCgg                 >  1:1321857/1‑76 (MQ=255)                                                      aGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCgg                 <  2:1325836/76‑1 (MQ=255)                                                      aGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCgg                 >  1:362368/1‑76 (MQ=255)                                                      aGTGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCgg                 <  2:419448/76‑1 (MQ=255)                                                        tGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCgggc               >  1:786660/1‑74 (MQ=255)                                                        tGCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGcc               >  1:771209/1‑76 (MQ=255)                                                         gCCGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCt              <  2:1434205/76‑1 (MQ=255)                                                          ccGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCtt             <  2:121510/76‑1 (MQ=255)                                                          ccGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCtt             >  2:1023276/1‑76 (MQ=255)                                                          ccGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCtt             <  2:609026/76‑1 (MQ=255)                                                          ccGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCtt             >  2:1416574/1‑76 (MQ=255)                                                           cGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGGCTTc            >  2:899570/1‑76 (MQ=255)                                                           cGAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTc            >  2:892697/1‑76 (MQ=255)                                                            gAAGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCg           >  2:314494/1‑76 (MQ=255)                                                             aaGTGGCGAAGTTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCgg          <  1:687390/76‑1 (MQ=255)                                                             aaGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCgg          <  2:492950/76‑1 (MQ=255)                                                              aGTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGc         <  2:1348651/76‑1 (MQ=255)                                                               gTGGCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCa        <  2:672397/76‑1 (MQ=255)                                                                tGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCCTCGGCAc       >  2:587759/1‑76 (MQ=255)                                                                 ggCGAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCAcc      >  1:1182396/1‑76 (MQ=255)                                                                    gAAATTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCaaa   <  1:777712/76‑1 (MQ=255)                                                                     aaaTTGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCaaaa  >  2:697730/1‑76 (MQ=255)                                                                         |                                                                        TGCATGGTGGCGTGCGACAGGTATAATCCACAACGTTTTCCGCATACCTCTTCAGTGCCGAAGTGGCGAAATCGGTAGACGCAGTTGATTCAAAATCAACCGTAGAAATACGTGCCGGTTCGAGTCCGGCCTTCGGCACCAAAA  >  REL606/4479577‑4479720

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