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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	191	1,643	0	T	C	98.2%	142.2 / ‑4.2	56	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (25/30);  total (26/30)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

GAAGCAGCCGCCGCACTTGCAGAACCAGCTGCAGCACTCTCGCTTTCGGCTGCTGCATCCTGACTGCTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTTGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCCGCAGCAGTTGCTGATCCAG  >  191/1545‑1740                                                                                                   |                                                                                                  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:311560/101‑1 (MQ=11)                            gCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCctggct                                                                      <  1:1323834/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:978532/1‑101 (MQ=14)                                    aCTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTg                                                              >  1:500769/1‑101 (MQ=14)                                    aCTCTCGCTTTGGGCCGCTGCATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTg                                                              <  1:402347/101‑1 (MQ=14)                                    aCTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTg                                                              >  1:568867/1‑101 (MQ=11)                                    aCTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTg                                                              <  1:1653789/101‑1 (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)                                            tttggtgcgctggATCCTGACTGTTTTCCGCCGCCGTTTCGCTTGCTTTGGCATTCGTTTCGTTGGCCTTCGATGCCGTCTGGCTGGACTTTGCGTTAGtt                                                      <  2:83297/93‑1 (MQ=14)                                              tgggCCGCTGCATCCTGGCTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTc                                                    <  2:280209/99‑1 (MQ=11)                                                   cgctgcATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCg                                               <  1:1000159/100‑1 (MQ=255)                                                   cgctgcATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCg                                               <  2:1092013/100‑1 (MQ=255)                                                    gctgcATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGt                                              >  2:1278739/1‑101 (MQ=255)                                                        cATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTc                                          >  2:1357826/1‑101 (MQ=255)                                                        cATCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTc                                          <  1:469345/101‑1 (MQ=255)                                                         aTCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTNCGTTAGTTTCGCTCGTCTTCg                                         >  2:789041/1‑101 (MQ=255)                                                         aTCCTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCg                                         >  1:1248837/1‑101 (MQ=255)                                                            cTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCgctg                                      >  1:859790/1‑101 (MQ=255)                                                            cTGACTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCgctg                                      >  2:817168/1‑101 (MQ=255)                                                               aCTGTTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCtt                                   >  2:439112/1‑101 (MQ=255)                                                                  gtttttCGCCGCAGTTTCGCTGCCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCt                                <  2:1563879/99‑1 (MQ=255)                                                                   cTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTTGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTg                               >  1:699987/1‑101 (MQ=255)                                                                        cgccgcAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGt                          <  1:741424/101‑1 (MQ=255)                                                                          ccgcAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTa                        >  2:1653783/1‑101 (MQ=255)                                                                           cgcAGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAg                       >  1:1031522/1‑101 (MQ=255)                                                                              aGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCTCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGctg                    >  1:792659/1‑99 (MQ=255)                                                                              aGTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGctg                    <  1:936759/101‑3 (MQ=255)                                                                               gTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGTTTTCTGGCTGTTAGCTgc                   >  2:186682/1‑101 (MQ=255)                                                                               gTTTCGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCTgc                   <  2:5962/101‑1 (MQ=255)                                                                                   cGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCTgcagca               <  1:851897/101‑1 (MQ=255)                                                                                   cGCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCCgcagca               >  2:1031462/1‑101 (MQ=255)                                                                                    gCTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCTGcagcag              <  1:745277/101‑1 (MQ=255)                                                                                     cTGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCTGCAGCAGt             >  2:1636284/1‑101 (MQ=255)                                                                                      tGGCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCTGCAGCAGtt            >  1:79316/1‑101 (MQ=255)                                                                                       ggCTTTGGCATTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCTGCAGCAGTTg           <  2:1467809/101‑1 (MQ=255)                                                                                                aTTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCTGCAGCAGTTGCTGATCCAg  >  1:1544546/1‑101 (MQ=255)                                                                                                aTTCGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCTGCAGCAGTTGCTGATCCAg  >  1:1358388/1‑101 (MQ=255)                                                                                                   |                                                                                                  GAAGCAGCCGCCGCACTTGCAGAACCAGCTGCAGCACTCTCGCTTTCGGCTGCTGCATCCTGACTGCTTTTCGCCGCAGTTTCGCTGGCTTTGGCATTTGTTTCGCTGGTCTTCGCTGCCGTCTGGCTGGACTTTGCGTTAGTTTCGCTCGTCTTCGCTGCTTTCTGGCTGTTAGCCGCAGCAGTTGCTGATCCAG  >  191/1545‑1740

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