Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	10	42,414	G→A	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	10	42,414	0	G	A	98.5%	227.1 / ‑4.3	68	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (33/34);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

CGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTAC  >  10/42316‑42504                                                                                                   |                                                                                           cGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACt                                                                                          <  1:243425/101‑1 (MQ=255)   gTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTaa                                                                                        <  2:1272003/101‑1 (MQ=255)    ttATAGTGGCATAATCAGCCGCAGGCATGGGGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAAt                                                                                       >  1:1143792/1‑101 (MQ=255)    gtataGGGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAAt                                                                                       <  1:1255783/100‑1 (MQ=255)     tataGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAAtt                                                                                      <  2:949994/101‑1 (MQ=255)     tataGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAAtt                                                                                      >  1:461656/1‑101 (MQ=255)      ataGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTg                                                                                     >  1:1058916/1‑101 (MQ=255)        aGTGGCATAATCAGCCGCAGGCATGGAGCCTGAGGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTggg                                                                                   >  2:1563362/1‑101 (MQ=255)         gTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGt                                                                                  <  1:261090/101‑1 (MQ=255)         gTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGt                                                                                  <  1:932274/101‑1 (MQ=255)          tGGCATAATCAGCCGCAGGCATGGAGCCTGAAGGCAACACCCTACGGAAACAAAAGACAACGGGAGATTTTCATGACCTGGGAATATGCACTAATTGGGtt                                                                                 <  2:874821/101‑1 (MQ=255)             cATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgt                                                                              <  1:499512/101‑1 (MQ=255)             cATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgt                                                                              <  2:718159/101‑1 (MQ=255)              aTAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                                                             <  2:1585335/101‑1 (MQ=255)              aTAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                                                             >  2:1548391/1‑101 (MQ=255)              aTAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                                                             >  2:880579/1‑101 (MQ=255)              aTAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                                                             >  2:261221/1‑101 (MQ=255)              aTAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtc                                                                             >  2:1696879/1‑101 (MQ=255)               tAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcg                                                                            <  2:1143787/101‑1 (MQ=255)                aaTCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAgtcgt                                                                           >  1:1108563/1‑101 (MQ=255)                      ccGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGcat                                                                     >  2:1523637/1‑101 (MQ=255)                           ggCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGTAATATGCACTAATTGGGTTAGTCGTCGGCATCatta                                                                <  2:1011812/101‑1 (MQ=255)                            gCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCattat                                                               <  1:1258049/101‑1 (MQ=255)                            gCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCattat                                                               <  1:1553692/101‑1 (MQ=255)                              aTGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTg                                                             >  1:1549443/1‑101 (MQ=255)                              aTGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTg                                                             >  1:1560937/1‑101 (MQ=255)                                ggAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGt                                                           <  1:1416925/101‑1 (MQ=255)                                 gAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTg                                                          <  1:169185/101‑1 (MQ=255)                                 gAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTg                                                          <  1:1362153/101‑1 (MQ=255)                                   gCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCt                                                        <  2:1272811/101‑1 (MQ=255)                                    ccTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtg                                                       >  1:1063322/1‑101 (MQ=255)                                     cTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgt                                                      >  1:350586/1‑101 (MQ=255)                                      tGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCtgtg                                                     >  2:94271/1‑101 (MQ=255)                                           cccacccccTGCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCAt                                                <  2:682838/95‑1 (MQ=255)                                            cAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATg                                               >  1:345927/1‑101 (MQ=255)                                            cAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATg                                               <  1:158292/101‑1 (MQ=255)                                            cAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATg                                               >  1:471190/1‑101 (MQ=255)                                                acCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGtt                                           >  2:1571545/1‑101 (MQ=255)                                                     aCGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGt                                      >  1:484331/1‑101 (MQ=255)                                                      cGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTa                                     <  2:1329800/101‑1 (MQ=255)                                                       ggAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTaa                                    >  2:1295743/1‑101 (MQ=255)                                                       ggAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTaa                                    >  1:1394920/1‑101 (MQ=255)                                                       ggAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTaa                                    <  2:468024/101‑1 (MQ=255)                                                          aaCAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCg                                 >  2:39389/1‑101 (MQ=255)                                                          aaCAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCg                                 <  1:979145/101‑1 (MQ=255)                                                            cAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTAGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTa                               <  1:678479/101‑1 (MQ=255)                                                             aaaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTaa                              >  1:717996/1‑101 (MQ=255)                                                              aaaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTaaa                             >  2:870795/1‑101 (MQ=255)                                                               aaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAAc                            >  2:976584/1‑101 (MQ=255)                                                               aaGACAACGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAAc                            <  2:1549437/101‑1 (MQ=255)                                                                      cGGGAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCa                     <  1:880584/101‑1 (MQ=255)                                                                        ggAGATGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGc                   >  2:1289237/1‑101 (MQ=255)                                                                            aTGTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACa               <  2:1580723/101‑1 (MQ=255)                                                                              gTTCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAgg             >  2:1619454/1‑101 (MQ=255)                                                                               ttCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTTATCGTAAACTACGCCAGCAACAGGc            >  2:162206/1‑101 (MQ=255)                                                                                tCATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCg           <  1:1107648/101‑1 (MQ=255)                                                                                 cATGACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGt          >  2:1341240/1‑101 (MQ=255)                                                                                    gACCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGc       >  1:1390356/1‑101 (MQ=255)                                                                                     aCCTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCa      <  1:1547210/101‑1 (MQ=255)                                                                                      ccTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAg     <  1:1427565/101‑1 (MQ=255)                                                                                       cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGt    <  2:587488/101‑1 (MQ=255)                                                                                       cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGt    >  2:74747/1‑101 (MQ=255)                                                                                       cTGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGt    <  1:1563368/101‑1 (MQ=255)                                                                                        tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCTAGCAACAGGCGTTGCAGTa   >  1:263140/1‑101 (MQ=255)                                                                                        tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTa   <  1:261221/101‑1 (MQ=255)                                                                                        tGGGAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTa   <  2:1273193/101‑1 (MQ=255)                                                                                         gggAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTAc  <  2:509244/101‑1 (MQ=255)                                                                                         gggAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTAc  >  2:208705/1‑101 (MQ=255)                                                                                         gggAATATGCACTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTAc  <  2:583197/101‑1 (MQ=255)                                                                                                   |                                                                                           CGGTTATAGTGGCATAATCAGCCGCAGGCATGGAGCCTGAAGCCAACACCCTACGGAAACAAAAGACAACGGGAGATGTTCATGACCTGGGAATATGCGCTAATTGGGTTAGTCGTCGGCATCATTATTGGTGCTGTGGCCATGCGTTTTGGTAATCGTAAACTACGCCAGCAACAGGCGTTGCAGTAC  >  10/42316‑42504

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