Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	19	30,831	(A)7→8	intergenic (–/–)	– / –	–/–

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	19	30,824	1	.	A	96.0%	168.6 / 0.8	50	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base . (1/1);  new base A (25/23);  total (26/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GGATGACTAAACTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑AAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTGAAGTGGGTAA  >  19/30727‑30919                                                                                                   |                                                                                                ggATGACTAAACTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑aaa                                                                                              <  2:901018/101‑1 (MQ=255) ggATGACTAAACTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑aaa                                                                                              >  2:171801/1‑101 (MQ=255) ggATGACTAAACTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑aaa                                                                                              <  2:1695816/101‑1 (MQ=255)   aTGACTAAACTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑aaaaa                                                                                            >  1:1598928/1‑101 (MQ=255)     gACTAAACTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑aaaaaaa                                                                                          <  1:332038/101‑1 (MQ=255)     gACTAAACTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑aaaaaaa                                                                                          <  1:1508113/101‑1 (MQ=255)       cTAAACTCGAGTGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑aaaaaaaag                                                                                        <  2:263062/101‑3 (MQ=255)       cTAAACTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑aaaaaaaag                                                                                        >  2:507145/1‑99 (MQ=255)        tAAACTCGCATGGGCGATAGGTGCGGTAGTACTGGTGACCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑aaaaaaaagt                                                                                       <  2:495998/101‑4 (MQ=255)           aCTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGAAGGATTACGAAAATGATGAAAAAAAAGTAtt                                                                                     <  1:67747/101‑1 (MQ=255)            cTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGTGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTc                                                                                    <  2:67495/101‑1 (MQ=255)            cTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTc                                                                                    <  1:960800/101‑1 (MQ=255)             tCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCt                                                                                   <  1:1672125/101‑1 (MQ=255)              cGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTg                                                                                  >  1:1203638/1‑101 (MQ=255)               gAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTgg                                                                                 <  2:292511/101‑1 (MQ=255)               gAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTgg                                                                                 <  2:1457578/101‑1 (MQ=255)                aaTGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGc                                                                                <  1:1285953/101‑1 (MQ=255)                 aTGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCg                                                                               >  1:1166005/1‑101 (MQ=255)                     gCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTc                                                                           >  2:536137/1‑101 (MQ=255)                      cGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCt                                                                          >  2:251888/1‑101 (MQ=255)                      cGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCt                                                                          <  1:255317/101‑1 (MQ=255)                       gATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTg                                                                         >  1:296487/1‑101 (MQ=255)                       gATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTg                                                                         >  2:123867/1‑101 (MQ=255)                          aGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTa                                                                      <  2:1586441/101‑1 (MQ=255)                           ggTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTAc                                                                     <  2:1136439/101‑1 (MQ=255)                            gTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACt                                                                    >  2:1188303/1‑101 (MQ=255)                              gCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCa                                                                  <  2:528126/101‑1 (MQ=255)                                ggTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCAcc                                                                <  2:658296/101‑1 (MQ=255)                                 gtagtaCTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑AAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAg                                                              >  2:1281154/1‑101 (MQ=255)                                 gtactaCTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCa                                                               <  1:578681/101‑1 (MQ=255)                                   agtaCTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGt                                                             >  1:987247/1‑101 (MQ=255)                                    gtaCTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGtt                                                            >  2:1098433/1‑101 (MQ=255)                                         ggTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGAATTCTGGCGTTTCTGTTACTCACCAGTTCTGCa                                                       >  1:1696077/1‑101 (MQ=255)                                          gTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAg                                                      <  1:1493691/101‑1 (MQ=255)                                           tGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAgc                                                     >  1:2522/1‑101 (MQ=255)                                               catcaGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAgcggcg                                                 <  1:86747/101‑1 (MQ=255)                                                atcaGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGgcgc                                                >  2:1103811/1‑101 (MQ=255)                                                atcaGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGgcgc                                                >  2:153829/1‑101 (MQ=255)                                                     ccTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTggcg                                           >  1:1189290/1‑101 (MQ=255)                                                        gCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCa                                        >  2:980009/1‑101 (MQ=255)                                                         cTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCAc                                       >  1:725259/1‑101 (MQ=255)                                                              gACTCTCGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAg                                  <  1:171801/101‑1 (MQ=255)                                                                   tcGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGAt                             >  2:821083/1‑101 (MQ=255)                                                                    cGAACCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGAtt                            <  2:178559/101‑1 (MQ=255)                                                                    cGAACCTTTTTGATGGATTACGAAAATGATG‑AAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTa                           <  2:305073/101‑1 (MQ=255)                                                                      aaCCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTAc                          <  1:426650/101‑1 (MQ=255)                                                                      aaCCTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTAc                          >  2:355725/1‑101 (MQ=255)                                                                         cTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGt                       >  1:829714/1‑101 (MQ=255)                                                                         cTTTTTGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGt                       <  1:980014/101‑1 (MQ=255)                                                                           ttttGATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCa                     <  2:1598922/101‑1 (MQ=255)                                                                               gATGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCg                 <  2:1332840/101‑1 (MQ=255)                                                                                aTGGATTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGt                >  1:1585221/1‑101 (MQ=255)                                                                                    aTTACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTg            >  2:1344734/1‑101 (MQ=255)                                                                                     ttACGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTGa           >  1:994988/1‑101 (MQ=255)                                                                                       aCGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTGAAg         >  1:414889/1‑101 (MQ=255)                                                                                        cGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTGAAGt        >  1:1337707/1‑101 (MQ=255)                                                                                        cGAAAATGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTGAAGt        <  1:227700/101‑1 (MQ=255)                                                                                          aaaaTGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTGAAGTgg      >  2:286068/1‑101 (MQ=255)                                                                                           aaaTGATGAAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTGAAGTggg     <  2:1166000/101‑1 (MQ=255)                                                                                              tgatgaAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTGAAGTGGGTaa  >  1:27412/1‑101 (MQ=255)                                                                                              tgatgaAAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGGGATTACCGTCAGCCGTTTTGAAGTGGGTaa  >  2:1542924/1‑101 (MQ=255)                                                                                                   |                                                                                                GGATGACTAAACTCGAATGGGCGATAGGTGCGGTAGTACTGGTGATCATCAGCCTGCTTGCGACTCTCGAACCTTTTTGATGGATTACGAAAATGATG‑AAAAAAAGTATTCTGGCGTTTCTGTTACTCACCAGTTCTGCAGCGGCGCTGGCGGCACCGCAGGTGATTACCGTCAGCCGTTTTGAAGTGGGTAA  >  19/30727‑30919

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