Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	157	3,893	T→G	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	157	3,893	0	T	G	98.5%	235.4 / ‑4.4	66	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (45/20);  total (45/21)
Fisher's exact test for biased strand distribution p-value = 3.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCA  >  157/3796‑3990                                                                                                  |                                                                                                  cgGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTa                                                                                                >  2:268034/1‑101 (MQ=255)  gggCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTAt                                                                                               >  2:1560615/1‑101 (MQ=255)   ggCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATc                                                                                              >  1:1656663/1‑101 (MQ=255)       tATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGt                                                                                          >  2:419044/1‑101 (MQ=255)       tATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGt                                                                                          <  2:435514/101‑1 (MQ=255)        aTTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTg                                                                                         >  2:1149515/1‑101 (MQ=255)            accaccGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTgccgc                                                                                     >  1:85193/1‑101 (MQ=255)            accaccGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTgccgc                                                                                     >  2:692038/1‑101 (MQ=255)                      gggCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCgg                                                                           >  2:1031775/1‑101 (MQ=255)                         cAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGcga                                                                        >  1:118806/1‑101 (MQ=255)                          aaTGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGcgac                                                                       >  1:1163514/1‑101 (MQ=255)                          aaTGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGcgac                                                                       <  1:1513581/101‑1 (MQ=255)                           aTGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGcgacg                                                                      >  1:1569294/1‑101 (MQ=255)                             gCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGCCGaa                                                                    >  2:727566/1‑101 (MQ=255)                              ccAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAat                                                                   <  2:561029/101‑1 (MQ=255)                               cAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAata                                                                  <  2:1133919/101‑1 (MQ=255)                                aaCACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAatat                                                                 >  1:1563112/1‑101 (MQ=255)                                aaCACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAatat                                                                 >  1:1206948/1‑101 (MQ=255)                                 acacGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATc                                                                <  1:738308/101‑1 (MQ=255)                                     gCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCa                                                            >  2:141163/1‑101 (MQ=255)                                      cTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCAc                                                           >  1:60425/1‑101 (MQ=255)                                        gttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCg                                                         <  2:852345/101‑1 (MQ=255)                                         ttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGCATATCGTCACCGc                                                        >  1:1700013/1‑101 (MQ=255)                                         ttgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGc                                                        >  1:541538/1‑101 (MQ=255)                                          tgttCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCt                                                       >  1:1654798/1‑101 (MQ=255)                                               aCCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCAc                                                  >  1:1694911/1‑101 (MQ=255)                                                  gCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAg                                               >  2:304577/1‑101 (MQ=255)                                                   cccGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGc                                              <  2:816891/101‑1 (MQ=255)                                                     cGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGa                                            <  2:986507/101‑1 (MQ=255)                                                      gAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGAt                                           <  1:1466851/101‑1 (MQ=255)                                                          gCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATc                                       >  1:605199/1‑101 (MQ=255)                                                           cTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCa                                      <  2:1049644/101‑1 (MQ=255)                                                              aaTAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCg                                   >  1:827632/1‑101 (MQ=255)                                                              aaTAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCg                                   <  1:692043/101‑1 (MQ=255)                                                               aTAGGGAAGTTCTTTGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCgg                                  <  2:168776/101‑1 (MQ=255)                                                               aTAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCgg                                  >  1:1099838/1‑101 (MQ=255)                                                                 aGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGa                                >  2:1653043/1‑101 (MQ=255)                                                                 aGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGa                                >  2:710519/1‑101 (MQ=255)                                                                  gggAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGAt                               >  2:1495713/1‑101 (MQ=255)                                                                  gggAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGAt                               >  2:361769/1‑101 (MQ=255)                                                                   ggAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATc                              <  1:397914/101‑1 (MQ=255)                                                                   ggAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATc                              <  2:85193/101‑1 (MQ=255)                                                                   ggAAGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATc                              >  1:933611/1‑101 (MQ=255)                                                                      aGTTCTTCGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGcc                           >  2:707640/1‑101 (MQ=255)                                                                         tcttcGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACg                        >  1:470730/1‑101 (MQ=255)                                                                           ttcGCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGtt                      >  1:1065284/1‑101 (MQ=255)                                                                              gcgttGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATc                   >  1:1427932/4‑101 (MQ=255)                                                                              gCTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATc                   <  2:22218/101‑1 (MQ=255)                                                                               cTTTGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCg                  >  1:1446158/1‑101 (MQ=255)                                                                                tttGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGt                 >  1:1641388/1‑101 (MQ=255)                                                                                 ttGTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTa                <  2:1198829/101‑1 (MQ=255)                                                                                   gTCGAGAAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGa              <  1:1024306/101‑1 (MQ=255)                                                                                   gTCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGa              >  2:1112186/1‑101 (MQ=255)                                                                                    tCCAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGAc             <  1:1226604/101‑1 (MQ=255)                                                                                     ccAGTAGCCCGACGTTATCCAGTTCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACa            >  1:542670/1‑101 (MQ=255)                                                                                      cAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAg           <  1:1031780/101‑1 (MQ=255)                                                                                      cAGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAg           >  2:629883/1‑101 (MQ=255)                                                                                       aGTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGt          >  2:387074/1‑101 (MQ=255)                                                                                        gTAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTa         >  1:446739/1‑101 (MQ=255)                                                                                         tAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTAc        >  2:1256814/1‑101 (MQ=255)                                                                                         tAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTAc        >  2:553043/1‑101 (MQ=255)                                                                                         tAGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACCAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTAc        <  2:113960/101‑1 (MQ=255)                                                                                          aGCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACg       <  2:630017/101‑1 (MQ=255)                                                                                           gCCCGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACgg      <  2:1569288/101‑1 (MQ=255)                                                                                            cccGACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGt     >  2:781616/1‑101 (MQ=255)                                                                                               gACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAGTCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCa  >  1:1121720/1‑101 (MQ=255)                                                                                               gACGTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCa  >  1:1246942/1‑101 (MQ=255)                                                                                                  |                                                                                                  CGGGCTTATTCACCACCGCGCGGGCAATGCCAACACGCTGTTGTTCACCGCCCGAAAGCTGAATAGGGAAGTTCTTCGCTTTGTCCAGTAGCCCGACTTTATCCAGTGCCGCCGACACCCGGCGACGAATATCGTCACCGCTGGCACCAGCGATAATCAGCGGGATCGCCACGTTATCGTAGACAGTACGGTCCA  >  157/3796‑3990

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