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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	178	1,007	0	C	T	97.0%	221.2 / NA	67	intergenic (–/–)	–/–	–/–
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (23/42);  total (23/44)
Fisher's exact test for biased strand distribution p-value = 5.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TATCCGGGGATTTGAAGCCGATAAATCCTGATAAATATCCATGAACGTAAAAATCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTGCCGCCAGTGCTCCTGCCCCTGACATTCTCGATCAGTTTCTGGATGAAAAAGAGGGGAATCACACCACTGCATACCGCGATGGCTCCGGCATATGGACCAT  >  178/912‑1103                                                                                                |                                                                                                 tatCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTg                                                                                             <  2:1313376/101‑1 (MQ=255) tatCCAGGGTTCAGAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTg                                                                                             <  2:1074918/101‑1 (MQ=255)     cagggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGcc                                                                                         >  1:1122916/10‑101 (MQ=255)      agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt                                                                                        <  2:973764/93‑1 (MQ=255)      agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt                                                                                        >  1:754578/9‑101 (MQ=255)      agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt                                                                                        >  2:124465/9‑101 (MQ=255)      agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt                                                                                        >  1:630300/9‑101 (MQ=255)        ggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGc                                                                                      >  2:241240/7‑101 (MQ=255)            cagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTgcgcctgcgcct                                                                                  >  2:152024/3‑101 (MQ=255)            aagAACCCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTgcgcctgcgcct                                                                                  <  2:322058/99‑1 (MQ=255)             agAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGcgcctgcgcctg                                                                                 >  2:773086/2‑101 (MQ=255)             agAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGcgcctgcgcctg                                                                                 >  2:1449259/2‑101 (MQ=255)              gAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGa                                                                                <  2:402908/101‑1 (MQ=255)               aaaccGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGat                                                                               <  2:280837/98‑2 (MQ=255)                 accGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGatat                                                                             <  1:224897/100‑1 (MQ=255)                   cGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATtc                                                                           <  2:1409888/101‑1 (MQ=255)                      tAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCg                                                                        <  1:588304/101‑1 (MQ=255)                      tAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCg                                                                        <  2:1012972/101‑1 (MQ=255)                      tAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCg                                                                        >  2:1463263/1‑101 (MQ=255)                         aTCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGacc                                                                     >  2:1595648/1‑99 (MQ=255)                         aTCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGacc                                                                     <  2:591473/101‑3 (MQ=255)                         aTCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGacc                                                                     >  1:1251641/1‑99 (MQ=255)                         aTCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGacc                                                                     >  1:1344063/1‑99 (MQ=255)                           ccTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAg                                                                   <  1:387059/101‑1 (MQ=255)                             tGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGtt                                                                 <  1:154656/101‑1 (MQ=255)                                tAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCt                                                              <  2:1560666/101‑1 (MQ=255)                                tAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCt                                                              <  2:1625755/101‑1 (MQ=255)                                 aaaTATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTg                                                             <  1:998209/101‑1 (MQ=255)                                  aaTATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTgg                                                            <  2:1059682/101‑1 (MQ=255)                                        cATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGaaa                                                      <  2:192131/101‑1 (MQ=255)                                        cATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGaaa                                                      <  2:1297011/101‑1 (MQ=255)                                         aTGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGaaaa                                                     <  2:896212/101‑1 (MQ=255)                                          tGAACGCAAAAATCAGATACGGCCTGTCGTCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGaaaaa                                                    >  1:874994/1‑101 (MQ=255)                                          tGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGaaaaa                                                    <  1:434989/101‑1 (MQ=255)                                           gAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAag                                                   <  1:371834/101‑1 (MQ=255)                                           gAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAag                                                   >  2:640000/1‑101 (MQ=255)                                             aCGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGGTTGCCGCCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAagaa                                                 <  2:820048/101‑2 (MQ=255)                                             aCGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAagaa                                                 <  2:672587/101‑2 (MQ=255)                                                caaaaaTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAaggt                                              <  2:981498/100‑2 (MQ=255)                                                caaaaaTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAaggt                                              >  1:1088425/2‑100 (MQ=255)                                                caaaaaTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAaggt                                              >  2:1190206/2‑100 (MQ=255)                                                 aaaaaTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAaggta                                             <  1:79404/101‑3 (MQ=255)                                                    aaTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTaacc                                          <  1:1398889/101‑3 (MQ=255)                                                    aaTCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTGCCGCCAGTGCTCCTGCCCCTGACATTCTCGATCAGTTTCTGGATGAAAAAGAGGGGAATc                                          <  2:1586210/101‑1 (MQ=255)                                                     aTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACca                                         <  2:1197176/101‑1 (MQ=255)                                                     aTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACca                                         <  2:204391/101‑1 (MQ=255)                                                       cAGATACGTCCTCTCGGCTGACGTCCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACcaca                                       <  2:245876/101‑1 (MQ=255)                                                           tACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCAcacaac                                   >  1:855779/1‑101 (MQ=255)                                                            acggCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCAcacaacg                                  <  1:1656166/98‑2 (MQ=255)                                                            acggCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCAcacaacg                                  <  1:152024/98‑2 (MQ=255)                                                                ccTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCAt                              >  1:1668162/1‑101 (MQ=255)                                                                ccTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCAt                              <  2:960522/101‑1 (MQ=255)                                                                 cTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATa                             >  1:1250666/1‑101 (MQ=255)                                                                 cTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATa                             <  1:1595654/101‑1 (MQ=255)                                                                  tGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATAc                            <  2:1488230/101‑1 (MQ=255)                                                                  tGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATAc                            <  1:124465/101‑1 (MQ=255)                                                                     cGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACcgc                         >  1:891989/1‑101 (MQ=255)                                                                        ctgccgTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGAt                      >  2:527882/6‑101 (MQ=255)                                                                        ctgccgTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGAt                      >  2:1629744/6‑101 (MQ=255)                                                                         tgccgTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATg                     <  2:932766/97‑1 (MQ=255)                                                                           ccgTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATggt                   <  2:533685/99‑2 (MQ=255)                                                                            cgTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATggtt                  <  2:1567722/100‑3 (MQ=255)                                                                            cgTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATggtt                  <  1:640005/100‑3 (MQ=255)                                                                                 tGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCggt             <  2:1088420/101‑2 (MQ=255)                                                                                     aCTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTatat         <  1:1411191/101‑1 (MQ=255)                                                                                      cTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATg        <  2:1325529/101‑1 (MQ=255)                                                                                       tGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATgg       <  2:1181018/101‑1 (MQ=255)                                                                                         aTTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGAc     >  2:1241583/1‑101 (MQ=255)                                                                                            gccgtcGGTGCGCCTGCGACTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCAt  <  1:113560/101‑1 (MQ=255)                                                                                                |                                                                                                 TATCCGGGGATTTGAAGCCGATAAATCCTGATAAATATCCATGAACGTAAAAATCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTGCCGCCAGTGCTCCTGCCCCTGACATTCTCGATCAGTTTCTGGATGAAAAAGAGGGGAATCACACCACTGCATACCGCGATGGCTCCGGCATATGGACCAT  >  178/912‑1103

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