Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	191	1,999	T→G	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	191	1,999	0	T	G	82.0%	179.9 / 21.8	61	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base T (6/5);  new base G (24/26);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.67e-01

CTATTTTTCGCGTTGTTTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGACGCTGCCGCCTGTGCGCTGTTAGCTGCATTCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGCGTTGCGTGC  >  191/1902‑2098                                                                                                  |                                                                                                    cTATTTTTCGCGTTGTTTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGt                                                                                                  >  1:1206536/1‑101 (MQ=255)     ttttCGCGTTGTTTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCt                                                                                              <  1:644135/101‑1 (MQ=255)         cgcgTTGTTTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAgg                                                                                          <  2:681480/101‑1 (MQ=255)           cgTTGTTTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGtt                                                                                        >  2:322734/1‑101 (MQ=255)           cgTTGTTTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGtt                                                                                        >  1:477888/1‑101 (MQ=255)           cgTTGTTTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGtt                                                                                        <  1:1624046/101‑1 (MQ=255)             ttgttTTCTGATGTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGtttt                                                                                      <  2:170456/101‑1 (MQ=255)              tgttTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTc                                                                                     <  2:1313388/101‑1 (MQ=255)              tgttTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTc                                                                                     >  2:743071/1‑101 (MQ=255)              tgttTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTc                                                                                     >  2:1134979/1‑101 (MQ=255)               gttTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTcg                                                                                    <  2:360556/101‑1 (MQ=255)               gttTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTcg                                                                                    >  1:46238/1‑101 (MQ=255)                 tttCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTcgcc                                                                                  >  1:60886/1‑101 (MQ=255)                 tttCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTcgcc                                                                                  <  1:961700/101‑1 (MQ=255)                      gATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGtt                                                                             <  2:982840/101‑1 (MQ=255)                       aTTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGttc                                                                            >  1:612245/1‑101 (MQ=255)                             tAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCgat                                                                      <  1:364344/101‑1 (MQ=255)                             tAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCgat                                                                      <  2:787245/101‑1 (MQ=255)                             tAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCgat                                                                      <  2:1050990/101‑1 (MQ=255)                              aGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCgatg                                                                     <  1:520903/101‑1 (MQ=255)                                ctgctgTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCgatgat                                                                   <  1:1090231/101‑1 (MQ=255)                                 tgctgTCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGatgatg                                                                  >  1:1278797/1‑101 (MQ=255)                                      tCGCGGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCTGCGTTCTTCGATGATGCCgct                                                             <  1:55742/101‑1 (MQ=255)                                         cgGAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCgctgct                                                          <  1:616430/101‑1 (MQ=255)                                          ggAGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGctgctg                                                         >  2:1366930/1‑101 (MQ=255)                                            aGTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGtt                                                       <  2:1278792/101‑1 (MQ=255)                                             gTTTGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGttt                                                      <  1:1430803/101‑1 (MQ=255)                                              tttGCCGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTc                                                     >  2:1488463/1‑101 (MQ=255)                                                   cGATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGCtt                                                >  2:1492651/1‑101 (MQ=255)                                                    gATGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGAttt                                               >  2:115120/1‑101 (MQ=255)                                                      tGCAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATctct                                             >  2:1474754/1‑101 (MQ=255)                                                        cAGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCttt                                           <  1:743076/101‑1 (MQ=255)                                                         aGTCTGCGAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTtg                                          >  2:556277/1‑101 (MQ=255)                                                                gAGGCCGCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTGTCTGATTTCTTTGCTGCCGc                                   >  1:580105/1‑101 (MQ=255)                                                                      gCTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGct                             >  2:845975/1‑101 (MQ=255)                                                                       cTGCCGCCTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGctc                            >  1:1674809/1‑101 (MQ=255)                                                                         gccgccTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTg                          <  2:900337/101‑1 (MQ=255)                                                                         gccgccTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTg                          >  2:456893/1‑101 (MQ=255)                                                                             ctttcgcgcTGTTACCTGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGagag                      <  2:95276/96‑1 (MQ=14)                                                                             ccTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGagag                      <  2:1206531/101‑1 (MQ=255)                                                                             ccTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGagag                      <  2:220604/101‑1 (MQ=255)                                                                             ccTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGagag                      <  1:845980/101‑1 (MQ=255)                                                                              cTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAgg                     >  2:477102/1‑101 (MQ=255)                                                                              cTGTGCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAgg                     <  2:1297079/101‑1 (MQ=255)                                                                               tgtgCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAggc                    <  2:859785/101‑1 (MQ=255)                                                                                 tgCGCTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAggcgg                  >  2:302246/1‑101 (MQ=255)                                                                                 cgcgcTGTTACCTGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAggcgg                  <  2:570725/100‑1 (MQ=17)                                                                                  gcgcTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAggcggc                 <  2:444398/101‑1 (MQ=255)                                                                                   cgcTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCt                <  1:1366936/101‑1 (MQ=255)                                                                                   cgcTGTTACCTGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCt                >  1:1182305/1‑101 (MQ=17)                                                                                    gcTGTTACCTGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCtt               <  2:1687890/101‑1 (MQ=18)                                                                                     cTGTTAGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTc              >  1:384712/1‑101 (MQ=255)                                                                                      tGTTACCTGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCa             <  2:1512521/101‑1 (MQ=18)                                                                                        ttACCTGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGc           >  2:1396220/1‑101 (MQ=17)                                                                                          accTGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGCGt         >  1:676536/3‑101 (MQ=18)                                                                                          aGCTGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGCGt         >  2:480106/1‑101 (MQ=255)                                                                                            cTGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGCGTTg       >  2:1077747/1‑101 (MQ=18)                                                                                             tGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGCGtttc      >  1:910530/1‑99 (MQ=21)                                                                                             tGCATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGCGTTGc      >  2:1451036/1‑101 (MQ=21)                                                                                             tGCATGCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGATTCAGCGTTGc      >  1:575234/1‑101 (MQ=255)                                                                                               cATTCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGCGTTGCGt    <  2:464749/101‑1 (MQ=21)                                                                                                 ttCGTCTCTGATGTTTTCGCCGCGTTCCTGGACGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGCGTTGCGtgc  >  1:1458626/1‑101 (MQ=25)                                                                                                  |                                                                                                    CTATTTTTCGCGTTGTTTTCTGATTTTTTAGCTGCTGTCGCGGAGTTTGCCGATGCAGTCTGCGAGGACGCTGCCGCCTGTGCGCTGTTAGCTGCATTCGTTTCTGAGGTTTTCGCCGCGTTCTTCGATGATGCCGCTGCTGTTTCTGATTTCTTTGCTGCCGCTGCGCTCTGAGAGGCGGCTTCAGCGTTGCGTGC  >  191/1902‑2098

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