Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	19	18,000	(T)8→7	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	19	17,993	0	T	.	96.9%	290.5 / 1.5	65	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base . (37/26);  total (37/28)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAATGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTATTTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTAAT  >  19/17893‑18089                                                                                                     |                                                                                                 gTAATGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAt                                                                                                  <  1:1026504/101‑1 (MQ=255) gTAATGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAt                                                                                                  <  2:325610/101‑1 (MQ=255) gTAATGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAt                                                                                                  <  2:825243/101‑1 (MQ=255)  tAATGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtt                                                                                                 <  2:278277/101‑1 (MQ=255)   aaTGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAttt                                                                                                >  2:929556/1‑101 (MQ=255)   aaTGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAttt                                                                                                <  1:677084/101‑1 (MQ=255)    aTGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtttt                                                                                               >  2:1640604/1‑101 (MQ=255)    aTGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtttt                                                                                               >  1:614044/1‑101 (MQ=255)      gACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtttttt                                                                                             >  2:1088656/1‑101 (MQ=255)      gACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtttttt                                                                                             <  2:1296721/101‑1 (MQ=255)      gACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtttttt                                                                                             >  2:1076455/1‑101 (MQ=255)        cGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtttttttc                                                                                           >  1:1246645/1‑100 (MQ=255)         gATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtttttttca                                                                                          >  1:1233334/1‑99 (MQ=255)          aTTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtttttttcat                                                                                         <  1:254634/101‑4 (MQ=255)           tttCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTAtttttttcata                                                                                        >  2:276114/1‑97 (MQ=255)                gCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATg                                                                                  >  1:1113249/1‑101 (MQ=255)                  aCTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGAt                                                                                <  1:1512531/101‑1 (MQ=255)                   cTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATc                                                                               <  2:578321/101‑1 (MQ=255)                    tGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCa                                                                              <  1:919781/101‑1 (MQ=255)                    tGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCa                                                                              <  2:792132/101‑1 (MQ=255)                     gaaagaaaTCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAg                                                                             >  2:383845/1‑101 (MQ=255)                      aaagaaaTCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAgt                                                                            >  2:678520/1‑101 (MQ=255)                      aaagaaaTCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTCA‑TTTTTTTCATAACATGATCAgt                                                                            >  2:1496547/1‑101 (MQ=255)                               gAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGAttt                                                                   <  1:545331/101‑1 (MQ=255)                               gAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGAttt                                                                   >  1:92919/1‑101 (MQ=255)                                aaaTGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGAtttt                                                                  >  1:629087/1‑101 (MQ=255)                                aaaTGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGAtttt                                                                  >  2:145637/1‑101 (MQ=255)                                aaaTGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGAtttt                                                                  >  2:1513209/1‑101 (MQ=255)                                 aaTGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGAttttt                                                                 <  1:284006/101‑1 (MQ=255)                                 aaTGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGAttttt                                                                 <  2:589708/101‑1 (MQ=255)                                  aTGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTa                                                                >  2:574687/1‑101 (MQ=255)                                  aTGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTa                                                                <  1:929799/101‑1 (MQ=255)                                   tGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTAc                                                               <  2:1411096/101‑1 (MQ=255)                                       gTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGaa                                                           >  2:1309273/1‑101 (MQ=255)                                       gTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGaa                                                           >  2:422845/1‑101 (MQ=255)                                           tGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGa                                                       <  2:67342/101‑1 (MQ=255)                                            gTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTCA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGaa                                                      <  1:1496553/101‑1 (MQ=255)                                              cAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGaaaa                                                    <  2:42224/101‑1 (MQ=255)                                               aGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGGTAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAAc                                                   >  2:774753/1‑101 (MQ=255)                                               aGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAAc                                                   <  2:620698/101‑1 (MQ=255)                                               aGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAAc                                                   >  2:1608108/1‑101 (MQ=255)                                               aGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAAc                                                   <  2:1246640/101‑1 (MQ=255)                                                gATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACg                                                  >  1:403977/1‑101 (MQ=255)                                                   atTCCGCTTGTGTGCCGTTTTAATGACAAAAGCAAATAAAAAAGTTGTT‑TTTTTTTTCATAACATGAT‑AGTGTCAGATTTTTACCGAATGGAAAACgatga                                              <  2:617739/101‑1 (MQ=255)                                                   atTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACgatg                                               <  2:1106940/101‑1 (MQ=255)                                                     tACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACgatgat                                             >  1:256409/1‑101 (MQ=255)                                                       cGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGAttt                                           <  2:502291/101‑1 (MQ=255)                                                       cGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGAttt                                           <  1:476484/101‑1 (MQ=255)                                                          cTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGAtttttt                                        >  2:963338/1‑101 (MQ=255)                                                          cTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGAtttttt                                        <  2:392895/101‑1 (MQ=255)                                                           tgtgtgCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGAttttttt                                       >  1:282885/1‑101 (MQ=255)                                                             tgtgCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTAt                                     <  1:1598370/101‑1 (MQ=255)                                                               tgCCGTGTTAATTACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCa                                   <  1:287250/101‑1 (MQ=255)                                                               tgCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCa                                   >  1:716827/1‑101 (MQ=255)                                                               tgCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCa                                   <  1:300746/101‑1 (MQ=255)                                                                 ccGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGt                                 >  2:646150/1‑101 (MQ=255)                                                                  cGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGtt                                >  2:729971/1‑101 (MQ=255)                                                                  cGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGtt                                <  1:678525/101‑1 (MQ=255)                                                                   gtgtTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTTTCAGttt                               >  2:1119711/1‑101 (MQ=255)                                                                   gtgtTAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGttt                               >  1:1559294/1‑101 (MQ=255)                                                                      ttAATGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGc                            <  2:1554522/101‑1 (MQ=255)                                                                         aTGACAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTCTTGCCGc                         >  1:428099/1‑101 (MQ=255)                                                                             cAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACtt                     >  2:104658/1‑101 (MQ=255)                                                                             cAAAAGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACtt                     >  2:526970/1‑101 (MQ=255)                                                                              aaaaGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACttt                    >  1:374467/1‑101 (MQ=255)                                                                                 aGCAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTgcg                 >  2:979926/1‑101 (MQ=255)                                                                                   cAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTgcgcg               <  1:847448/101‑1 (MQ=255)                                                                                   cAGATAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTgcgcg               >  1:905210/1‑101 (MQ=255)                                                                                      aTAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCtt            <  1:574692/101‑1 (MQ=255)                                                                                      aTAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCtt            <  2:1233329/101‑1 (MQ=255)                                                                                      aTAAAAAAGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCtt            >  2:1018085/1‑101 (MQ=255)                                                                                        aaaaaaGTTGTTA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCtttt          >  1:481073/1‑101 (MQ=255)                                                                                              gttgttA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTa    >  2:1048919/1‑101 (MQ=255)                                                                                               ttgttATTTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTa    <  2:213511/101‑1 (MQ=255)                                                                                               ttgttA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTaa   >  1:998855/1‑101 (MQ=255)                                                                                               ttgttA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTaa   >  1:461103/1‑101 (MQ=255)                                                                                               ttgttA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTaa   >  1:21652/1‑101 (MQ=255)                                                                                                tgttA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTAAt  >  1:931729/1‑101 (MQ=255)                                                                                                tgttA‑TTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTAAt  >  2:125666/1‑101 (MQ=255)                                                                                                tgttA‑TTTTTTTCATAACATGATCAGGGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTAAt  >  1:149698/1‑101 (MQ=255)                                                                                                     |                                                                                                 GTAATGACGATTTCCGCACTGAAAGAAATCGAAATGCAGTTTTGTCAGATATTACGCCTGTGTGCCGTGTTAATGACAAAAGCAGATAAAAAAGTTGTTATTTTTTTTCATAACATGATCAGTGTCAGATTTTTACCGAATGGAAAACGATGATTTTTTTATCAGTTTTGCCGCACTTTGCGCGCTTTTCCCGTAAT  >  19/17893‑18089

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