Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	35	23,402	A→C	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	35	23,402	0	A	C	98.3%	211.6 / ‑4.4	60	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base C (34/25);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.95e-01

CGGCGTTTCCTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTAGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGACGCGCACGCCGCGACGAA  >  35/23303‑23499                                                                                                    |                                                                                                  cGGCGTTTCCTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTAg                                                                                                  <  2:82202/101‑1 (MQ=255)  ggCGTTTCCTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGa                                                                                                 >  1:1202152/1‑101 (MQ=255)   gCGTTTCCTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGAt                                                                                                >  2:676062/1‑101 (MQ=255)    cGTTTCCTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGAtt                                                                                               >  1:1206535/1‑101 (MQ=255)     gTTTCCTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGAttt                                                                                              >  1:206030/1‑101 (MQ=255)      tttCCTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTa                                                                                             >  2:504135/1‑101 (MQ=255)       ttCCTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTNATTTCATTTTTCAGNANTGCCAGCAGGGTCGATTTAc                                                                                            >  2:226626/1‑101 (MQ=255)         ccTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACca                                                                                          <  2:1452415/101‑1 (MQ=255)         ccTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACca                                                                                          >  2:1568976/1‑101 (MQ=255)           tGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACcaca                                                                                        <  1:770549/101‑1 (MQ=255)            gATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAg                                                                                       >  1:671705/1‑101 (MQ=255)               tCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCg                                                                                    <  1:104326/101‑1 (MQ=255)                cACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGt                                                                                   >  1:644147/1‑101 (MQ=255)                  cccACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGttt                                                                                 <  1:714846/101‑1 (MQ=255)                     aCGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTa                                                                              <  2:928962/101‑1 (MQ=255)                      cGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTAc                                                                             <  2:671062/101‑1 (MQ=255)                       gCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTAcc                                                                            >  1:1080947/1‑101 (MQ=255)                        ccAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTAccc                                                                           >  1:623731/1‑101 (MQ=255)                         cAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACcca                                                                          <  1:676067/101‑1 (MQ=255)                         cAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACcca                                                                          <  1:76587/101‑1 (MQ=255)                          aGTTGCCAGCTTCCCGGAAAGGGGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACccac                                                                         >  1:488248/1‑101 (MQ=255)                           gTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACccacc                                                                        >  2:32790/1‑101 (MQ=255)                                 aGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCg                                                                  >  2:999472/1‑101 (MQ=255)                                   cTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTCCCACAGCCGTTTTTACCCACCAGGCCGAc                                                                <  2:844983/101‑1 (MQ=255)                                      cccGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACttt                                                             <  2:742617/101‑1 (MQ=255)                                       ccGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTc                                                            >  1:919661/1‑101 (MQ=255)                                       ccGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTc                                                            >  2:1032067/1‑101 (MQ=255)                                        cGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCt                                                           >  1:1006437/1‑101 (MQ=255)                                         ggAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTg                                                          <  2:877217/101‑1 (MQ=255)                                            aaGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGccc                                                       <  1:1113749/101‑1 (MQ=255)                                              ggTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCgg                                                     >  2:1296246/1‑101 (MQ=255)                                                 gtAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGt                                                  >  1:517863/1‑101 (MQ=255)                                                   aGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTg                                                <  2:1297146/101‑1 (MQ=255)                                                     cTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGAt                                              <  2:1182629/101‑1 (MQ=255)                                                      tGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGAtg                                             >  2:438376/1‑101 (MQ=255)                                                       gccgccGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGAtgg                                            <  1:32790/101‑1 (MQ=255)                                                          gccgTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGAtggtgg                                         <  2:209493/101‑1 (MQ=255)                                                           ccgTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATggtggc                                        >  1:705128/1‑101 (MQ=255)                                                             gTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATggtggcgg                                      <  2:488243/101‑1 (MQ=255)                                                                 gcgcTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATggtggcggtggc                                  >  2:876285/1‑101 (MQ=255)                                                                  cgcTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCa                                 >  1:1223499/1‑101 (MQ=255)                                                                    cTGATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCAtt                               >  2:1719240/1‑101 (MQ=255)                                                                      gATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTAt                             <  1:1296252/101‑1 (MQ=255)                                                                      gATTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTAt                             >  1:745493/1‑101 (MQ=255)                                                                       aTTTCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATc                            <  1:1250659/101‑1 (MQ=255)                                                                        tttCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATcc                           >  1:1228823/1‑101 (MQ=255)                                                                        tttCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATcc                           >  1:490984/1‑101 (MQ=255)                                                                        tttCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATcc                           <  1:1568982/101‑1 (MQ=255)                                                                         ttCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCca                          <  1:914104/101‑1 (MQ=255)                                                                         ttCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCca                          >  2:1079873/1‑101 (MQ=255)                                                                         ttCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCca                          <  2:671700/101‑1 (MQ=255)                                                                          tCATTTTTCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCcag                         <  2:123909/101‑1 (MQ=255)                                                                                ttCAGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGAc                   <  1:1037920/101‑1 (MQ=255)                                                                                   aGCAATGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGAcgcg                >  1:628542/1‑101 (MQ=255)                                                                                       aTGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGAcacgcacg            >  1:1078503/1‑101 (MQ=255)                                                                                       aTGCCAGCAGGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGACGCGCAcg            >  2:500182/1‑101 (MQ=255)                                                                                           cagcagGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGACGCGCAcgccgc        >  1:233636/1‑101 (MQ=255)                                                                                            agcagGGTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGACGCGCACGCcgcg       <  2:1223494/101‑1 (MQ=255)                                                                                                gggTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGACGCGCACGCCGcgacga   >  2:299956/1‑101 (MQ=255)                                                                                                 ggTCGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGACGCGCACGCCGcgacgca  >  1:39038/1‑99 (MQ=255)                                                                                                    |                                                                                                  CGGCGTTTCCTGATTCACCCACGCCAGTTGCCAGCTTCCCGGAAAGGTGTAGCTGCCGCCGTCGGCGCTGATTTCATTTTTCAGCAATGCCAGCAGGGTAGATTTACCACAGCCGTTTTTACCCACCAGGCCGACTTTCTGCCCGGGGTTGATGGTGGCGGTGGCATTATCCAGCAGGACGCGCACGCCGCGACGAA  >  35/23303‑23499

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