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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	191	1,591	0	C	G	83.9%	124.2 / 22.1	56	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base C (6/3);  new base G (15/32);  total (21/35)
Fisher's exact test for biased strand distribution p-value = 6.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.06e-01

TGGTTTCACTGGTTTTTGCAGCTTTTTGACTGTTTGCTGATGCAGTGGCAGAAGCAGCCGCCGCACTTGCAGAACCAGCTGCAGCACTCTCGCTTTCGGCTGCTGCATCCTGACTGCTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTTGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTA  >  191/1495‑1685                                                                                                 |                                                                                               tGGTTTCACTGGTTTTTGCAGCTTTTTGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGgcc                                                                                            <  1:1468815/101‑2 (MQ=255) tGGTTTCACTGGTTTTTGCAGCTTTTTGACTGTTGGCTGATCCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTTCGGCACTCTCGCTTTGGgcc                                                                                            <  2:1366163/101‑2 (MQ=255)    tttCACTGGTTTTTGCAGCTTTTTGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGgccgcn                                                                                         <  1:710233/101‑2 (MQ=255)       cACTGGTTTTTGCAGCTTTTTGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCa                                                                                      <  2:1719712/101‑1 (MQ=255)         cTGGTTTTTGCAGCTTTTTGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATc                                                                                    >  1:66010/1‑101 (MQ=255)         cTGGTTTTTGCAGCTTTTTGACTGTTGGCTGATGCAATGGCATAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATc                                                                                    >  2:1299314/1‑101 (MQ=255)           ggTTTTTGCAGCTTTTTGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCt                                                                                  <  1:429185/101‑1 (MQ=255)            gTTTTTGCAGGTTATTGACTGTCGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGTCCGCTGCATCCTg                                                                                 <  2:876192/101‑1 (MQ=255)                ttGCAGCTTTTTGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTg                                                                             <  1:874089/101‑1 (MQ=255)                     gCTTTTTGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTtttt                                                                        <  1:1561670/101‑1 (MQ=255)                     gCTTTTTGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTtttt                                                                        <  1:1062566/101‑1 (MQ=255)                      cTTTTTGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTc                                                                       <  2:762976/101‑1 (MQ=255)                          ttGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTcgccg                                                                   <  1:1275731/101‑1 (MQ=255)                          ttGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTcgccg                                                                   >  1:1362203/1‑101 (MQ=255)                           tGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTcgccgc                                                                  <  2:944508/101‑1 (MQ=255)                           tGACTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTcgccgc                                                                  >  1:888214/1‑101 (MQ=255)                             aCTGTTGGCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAg                                                                <  2:1644338/101‑1 (MQ=255)                              cTGTTGGCTGATGCAGGGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCCCTCTCACTTTGGGCCGCTTCATCCTGACTGTTTTTCGCCGCAGt                                                               <  2:1583912/101‑1 (MQ=255)                                 ttggCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTc                                                            >  1:414609/4‑101 (MQ=25)                                 ttggCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTc                                                            <  1:1246806/98‑1 (MQ=25)                                 ttggCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTc                                                            <  2:1104751/98‑1 (MQ=25)                                  tggCTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCg                                                           >  1:897062/3‑101 (MQ=25)                                     cTGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTg                                                        >  2:402342/1‑101 (MQ=25)                                      tGATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTgg                                                       >  1:633231/1‑101 (MQ=25)                                       gATGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGc                                                      <  1:694332/101‑1 (MQ=25)                                         tGCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCtt                                                    >  1:163311/1‑101 (MQ=25)                                          gCAGTGGCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCttt                                                   <  1:415422/101‑1 (MQ=25)                                                gCAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCAtt                                             <  2:501253/101‑1 (MQ=21)                                                 cAGAAGAAGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCAttc                                            >  1:5962/1‑100 (MQ=21)                                                   gaagaaGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGt                                          >  1:673271/1‑101 (MQ=18)                                                   gaagaaGCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGt                                          <  2:810679/101‑1 (MQ=18)                                                    aagaagCCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGtt                                         <  2:983242/101‑1 (MQ=18)                                                     agaagcCGCCGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGttt                                        <  1:619794/98‑1 (MQ=18)                                                         gccgccGCACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCt                                    <  2:1525231/101‑1 (MQ=18)                                                               gcACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCtt                              <  2:295074/101‑1 (MQ=17)                                                                cACTTGCAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTc                             <  2:538872/101‑1 (MQ=17)                                                                      cAGAACCAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGcc                       <  2:505075/101‑1 (MQ=14)                                                                           ccAGCTGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTctg                  <  1:433886/101‑1 (MQ=14)                                                                            cggCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCctgg                 >  1:681485/3‑101 (MQ=11)                                                                              gCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCctggct               <  1:1323834/101‑1 (MQ=11)                                                                              gCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCctggct               <  1:311560/101‑1 (MQ=11)                                                                               cTGCGGCACTCTCGCTGTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTctggctg              <  1:957028/101‑1 (MQ=14)                                                                                tGCGGCACTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTctggctgg             >  2:745272/1‑101 (MQ=14)                                                                                 gccgcaCTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGa            >  2:524441/4‑101 (MQ=11)                                                                                   ggcaCTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACt          <  2:1400008/100‑1 (MQ=14)                                                                                    gcaCTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACtt         <  2:1593306/101‑1 (MQ=14)                                                                                    gcaCTCTCGCTTTCGGCTGCTGCAGCCTCACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACtt         <  2:438723/101‑1 (MQ=11)                                                                                     caCTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACttt        <  1:1649020/101‑1 (MQ=14)                                                                                     caCTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACttt        >  1:1268286/1‑101 (MQ=14)                                                                                      aCTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTg       >  1:500769/1‑101 (MQ=14)                                                                                      aCTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTg       <  1:402347/101‑1 (MQ=14)                                                                                      aCTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTg       >  1:978532/1‑101 (MQ=14)                                                                                      aCTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTg       <  1:1653789/101‑1 (MQ=11)                                                                                      aCTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTg       >  1:568867/1‑101 (MQ=11)                                                                                      aCTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTg       >  2:1090226/1‑101 (MQ=11)                                                                                       ctctCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGc      <  2:657565/101‑1 (MQ=14)                                                                                         ctcGCTTTGGGCCGCTGCATCCTGCCTGTTTTTCGCCGCAGTTTCGCTGGCTTGGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGt    <  1:50596/101‑1 (MQ=11)                                                                                         ctcGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGt    >  1:1313394/1‑101 (MQ=11)                                                                                           cGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGGTGGACTTTGCGTTa  >  2:107575/1‑101 (MQ=14)                                                                                           cGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTa  <  1:1317147/101‑1 (MQ=14)                                                                                                 |                                                                                               TGGTTTCACTGGTTTTTGCAGCTTTTTGACTGTTTGCTGATGCAGTGGCAGAAGCAGCCGCCGCACTTGCAGAACCAGCTGCAGCACTCTCGCTTTCGGCTGCTGCATCCTGACTGCTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTTGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTA  >  191/1495‑1685

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