Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	32	22,712	G→T	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	32	22,712	0	G	T	98.4%	212.1 / ‑3.7	62	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (21/40);  total (21/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.26e-01

ACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCAT  >  32/22614‑22804                                                                                                   |                                                                                             acacAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTa                                                                                            >  2:111849/1‑101 (MQ=255) acacAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTa                                                                                            <  2:1246048/101‑1 (MQ=255)   acaGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAg                                                                                          >  2:1257343/1‑101 (MQ=255)    caGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAgg                                                                                         >  2:44870/1‑101 (MQ=255)     aGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGa                                                                                        <  1:514832/101‑1 (MQ=255)       atatCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAAt                                                                                      >  1:814283/1‑101 (MQ=255)       atatCGCAAGGCAAATGGCGACCAAGCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAAt                                                                                      <  2:627482/101‑1 (MQ=255)        tatCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGTCCAAAAAGCGGCACCAGGAATTAAGGAATg                                                                                     <  1:242379/101‑1 (MQ=255)        tatCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATg                                                                                     <  2:1024266/101‑1 (MQ=255)         atCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGt                                                                                    <  1:866028/101‑1 (MQ=255)           cGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGtct                                                                                  <  2:100259/101‑1 (MQ=255)           cGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGtct                                                                                  <  1:1209008/101‑1 (MQ=255)           cGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGtct                                                                                  <  2:715921/101‑1 (MQ=255)             cAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTctct                                                                                <  1:864702/101‑1 (MQ=255)             cAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTctct                                                                                >  1:241669/1‑101 (MQ=255)              aaGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTg                                                                               <  1:492753/101‑1 (MQ=255)              aaGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTg                                                                               >  1:418897/1‑101 (MQ=255)               aGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTgc                                                                              >  2:341772/1‑101 (MQ=255)                  cAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTgccgc                                                                           <  1:1335408/101‑1 (MQ=255)                   aaaTGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTgccgcc                                                                          <  1:577229/101‑1 (MQ=255)                      tGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGt                                                                       <  2:424239/101‑1 (MQ=255)                          gACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAg                                                                   <  2:226351/101‑1 (MQ=255)                           aCCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAgg                                                                  <  1:21543/101‑1 (MQ=255)                           aCCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAgg                                                                  <  1:797028/101‑1 (MQ=255)                           aCCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAACTAAGGAATGTCTCTGCAGCCGGTAAAgg                                                                  >  1:1610255/1‑101 (MQ=255)                             cAATCAGTTTGCGGTAGGCAAATACACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCg                                                                <  2:256281/101‑1 (MQ=255)                              aaTCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGa                                                               >  2:1421009/1‑101 (MQ=255)                                 cAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACAt                                                            <  1:207746/101‑1 (MQ=255)                                 cAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACAt                                                            <  1:406045/101‑1 (MQ=255)                                  aGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACAtt                                                           <  1:1634811/101‑1 (MQ=255)                                  aGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACAtt                                                           <  2:239670/101‑1 (MQ=255)                                        cGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCa                                                     >  2:966660/1‑101 (MQ=255)                                          gTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGc                                                   >  1:1532480/1‑101 (MQ=255)                                           tAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGcc                                                  <  1:986080/101‑1 (MQ=255)                                              gCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgc                                               >  2:919703/1‑101 (MQ=255)                                               cAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgcc                                              >  2:203250/1‑101 (MQ=255)                                                aaaTGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgcca                                             <  1:1257348/101‑1 (MQ=255)                                                aaaTGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgcca                                             >  1:230813/1‑101 (MQ=255)                                                  aTGCACCTACAATCGCGCGAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgccagc                                           <  2:1639811/101‑1 (MQ=255)                                                    gCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgccagcca                                         >  2:61134/1‑101 (MQ=255)                                                    gCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCgccagcca                                         <  2:1004145/101‑1 (MQ=255)                                                     cACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGccagccag                                        <  2:418892/101‑1 (MQ=255)                                                       ccTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGc                                      >  1:28021/1‑101 (MQ=255)                                                        cTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCa                                     <  1:111849/101‑1 (MQ=255)                                                         tACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCaa                                    <  1:991710/101‑1 (MQ=255)                                                          aCAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCATGAATTAAGGAATGTCTTTGCCGCCGGTAAAGGCGACATTGTCACAGCCCGCCAGCCCGGCaaa                                   >  2:609549/1‑101 (MQ=255)                                                          aCAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCaaa                                   <  1:692357/101‑1 (MQ=255)                                                                 cgcCAACGCCAGGGCCAAAAAGCGGCACCAGGAAGTAAGGAATGTCTAAGCAGCCGGTAAAGGCGACATTGTCCCAGCCCGCCAGCCAGGCAAAGACtttg                            >  2:1411756/1‑100 (MQ=255)                                                                 cgcCAAAGATGGGGCGAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCCAACAGGCAAAGACTTTc                            <  2:1610249/101‑1 (MQ=255)                                                                   ccAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCgg                          >  2:610156/1‑101 (MQ=255)                                                                   ccAACGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCgg                          >  1:331556/1‑101 (MQ=255)                                                                     aaCGATAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAATGACTTTCGGAc                        <  2:531948/101‑1 (MQ=255)                                                                          tAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAg                   <  1:1280989/101‑1 (MQ=255)                                                                          tAGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAg                   <  2:28021/101‑1 (MQ=255)                                                                           aGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGt                  <  2:1285955/101‑1 (MQ=255)                                                                           aGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGt                  <  2:241669/101‑1 (MQ=255)                                                                           aGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGt                  <  2:413134/101‑1 (MQ=255)                                                                           aGGGCCAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGt                  >  2:1295096/1‑101 (MQ=255)                                                                               ccAAAAAGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCAcg              <  2:814278/101‑1 (MQ=255)                                                                                 aaaaaGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCAcgcg            <  1:1282150/101‑1 (MQ=255)                                                                                    aaGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTg         >  1:775347/1‑101 (MQ=255)                                                                                    aaGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTg         <  2:1676045/101‑1 (MQ=255)                                                                                     aGCGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTgg        <  1:203250/101‑1 (MQ=255)                                                                                      gCGGCACCAGGAATGAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTggg       <  2:1049241/101‑1 (MQ=255)                                                                                       cGGCACCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGt      <  2:1579373/101‑1 (MQ=255)                                                                                           aCCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCAt  >  1:176014/1‑101 (MQ=255)                                                                                           aCCAGGAATTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCAt  >  2:588687/1‑101 (MQ=255)                                                                                                   |                                                                                             ACACAGATATCGCAAGGCAAATGGCGACCAATCAGTTTGCGGTAGGCAAATGCACCTACAATCGCGCCAACGATAGGGCCAAAAAGCGGCACCAGGAAGTAAGGAATGTCTCTGCCGCCGGTAAAGGCGACATTGCCCCAGCCCGCCAGCCAGGCAAAGACTTTCGGACCGAAGTCACGCGCTGGGTTCAT  >  32/22614‑22804

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 23 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑