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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1781,0140TG98.5% 235.0 / NA 67intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (22/44);  total (22/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGGGGATTTGAAGCCGATAAATCCTGATAAATATCCATGAACGTAAAAATCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTGCCGCCAGTGCTCCTGCCCCTGACATTCTCGATCAGTTTCTGGATGAAAAAGAGGGGAATCACACCACTGCATACCGCGATGGCTCCGGCATATGGACCATCTGTCGGG  >  178/916‑1111
                                                                                                  |                                                                                                 
cagggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGcc                                                                                                 >  1:1122916/10‑101 (MQ=255)
 agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt                                                                                                >  2:124465/9‑101 (MQ=255)
 agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt                                                                                                <  2:973764/93‑1 (MQ=255)
 agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt                                                                                                >  1:630300/9‑101 (MQ=255)
 agggttcagAAACCGATAAATCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCt                                                                                                >  1:754578/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                                                                             >  1:1344063/1‑99 (MQ=255)
                    aTCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGacc                                                                             <  2:591473/101‑3 (MQ=255)
                    aTCCTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGacc                                                                             >  1:1251641/1‑99 (MQ=255)
                      ccTGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAg                                                                           <  1:387059/101‑1 (MQ=255)
                        tGATAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGtt                                                                         <  1:154656/101‑1 (MQ=255)
                           tAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCt                                                                      <  2:1625755/101‑1 (MQ=255)
                           tAAATATCCATGAACGCAAAAATCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCt                                                                      <  2:1560666/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                                                      >  1:1088425/2‑100 (MQ=255)
                                           caaaaaTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAaggt                                                      <  2:981498/100‑2 (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:152024/98‑2 (MQ=255)
                                                       acggCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCAcacaacg                                          <  1:1656166/98‑2 (MQ=255)
                                                           ccTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCAt                                      >  1:1668162/1‑101 (MQ=255)
                                                           ccTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCAt                                      <  2:960522/101‑1 (MQ=255)
                                                            cTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATa                                     <  1:1595654/101‑1 (MQ=255)
                                                            cTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATa                                     >  1:1250666/1‑101 (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:1629744/6‑101 (MQ=255)
                                                                   ctgccgTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGAt                              >  2:527882/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)
                                                                                        ccgtcggTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATc         <  2:540766/95‑1 (MQ=255)
                                                                                         cgtcggTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGAGCATCt        >  1:773946/6‑101 (MQ=255)
                                                                                         cgtcggTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCt        <  2:1719668/96‑1 (MQ=255)
                                                                                         cgtcggTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCt        <  2:899452/96‑1 (MQ=255)
                                                                                            cggTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTc     >  1:32139/3‑101 (MQ=255)
                                                                                            cggTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTc     <  2:221257/99‑1 (MQ=255)
                                                                                               tGCGCCAGCGCCTGATATTCTCGAGCAGTTTCTGTATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCgtg  <  1:1631113/101‑3 (MQ=255)
                                                                                                  |                                                                                                 
CGGGGATTTGAAGCCGATAAATCCTGATAAATATCCATGAACGTAAAAATCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTGCCGCCAGTGCTCCTGCCCCTGACATTCTCGATCAGTTTCTGGATGAAAAAGAGGGGAATCACACCACTGCATACCGCGATGGCTCCGGCATATGGACCATCTGTCGGG  >  178/916‑1111

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