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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1781,0560GA98.6% 234.5 / ‑4.5 69intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (28/40);  total (28/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

TAAAAATCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTGCCGCCAGTGCTCCTGCCCCTGACATTCTCGATCAGTTTCTGGATGAAAAAGAGGGGAATCACACCACTGCATACCGCGATGGCTCCGGCATATGGACCATCTGTCGGGGAGCAACGATGGTGGACGGTAAACCTGTTATTCCGGGAATGA  >  178/959‑1153
                                                                                                 |                                                                                                 
caaaaaTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAaggt                                                                                                >  2:1190206/2‑100 (MQ=255)
caaaaaTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAaggt                                                                                                <  2:981498/100‑2 (MQ=255)
caaaaaTCAGATACGGCCTGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAaggt                                                                                                >  1:1088425/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                                                                              <  1:124465/101‑1 (MQ=255)
                  tGTCGGCTGCCGTTCTGGCACTGATTGCCGTCGGTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATAc                                                                              <  2:1488230/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                                               <  2:221257/99‑1 (MQ=255)
                                                 cggTGCGCCTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTc                                               >  1:32139/3‑101 (MQ=255)
                                                    tGCGCCAGCGCCTGATATTCTCGAGCAGTTTCTGTATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCgtg                                            <  1:1631113/101‑3 (MQ=255)
                                                      cgccTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCgtggt                                          <  1:1446082/99‑2 (MQ=255)
                                                      cgccTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCgtggt                                          <  1:558763/99‑2 (MQ=255)
                                                       gccTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCgtggtg                                         <  1:1308888/100‑3 (MQ=255)
                                                        ccTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGc                                        >  2:1425292/1‑101 (MQ=255)
                                                        ccTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGc                                        >  1:940400/1‑101 (MQ=255)
                                                        ccTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGGAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGc                                        >  2:893566/1‑101 (MQ=255)
                                                         cTGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGcc                                       >  1:1576104/1‑100 (MQ=255)
                                                          tGCGCCTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGcca                                      >  2:239314/1‑99 (MQ=255)
                                                            cgccTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAATGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCgtggttccaca                                    >  1:924926/3‑94 (MQ=255)
                                                              ccTGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAAt                                  >  2:1557330/1‑101 (MQ=255)
                                                                tGATATTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAAtgg                                >  1:1057150/1‑101 (MQ=255)
                                                                  ataTTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAAtggtg                              >  1:375192/3‑101 (MQ=255)
                                                                    aTTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGa                            >  1:306868/1‑101 (MQ=255)
                                                                    aTTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGa                            >  2:260702/1‑101 (MQ=255)
                                                                    aTTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGa                            >  1:1003558/1‑101 (MQ=255)
                                                                    aTTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGa                            >  2:1508765/1‑101 (MQ=255)
                                                                    aTTCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGCATCTGGACCATCTGCCGTGGTGCCATCATGGTGGa                            >  2:553664/1‑101 (MQ=255)
                                                                     ttCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGat                           <  2:1344057/101‑2 (MQ=255)
                                                                     ttCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGat                           <  2:863412/101‑2 (MQ=255)
                                                                     ttCTCGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGat                           >  1:238827/1‑100 (MQ=255)
                                                                      tctcGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGatg                          <  1:427833/101‑3 (MQ=255)
                                                                         cGACCAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTa                       <  2:1385940/101‑1 (MQ=255)
                                                                          gaccAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTaa                      >  2:896720/4‑101 (MQ=255)
                                                                           accAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTaag                     <  2:779067/99‑2 (MQ=255)
                                                                           accAGTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTaag                     <  2:32139/99‑2 (MQ=255)
                                                                               gTTTCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGcccg                 <  2:1430999/101‑3 (MQ=255)
                                                                                tttCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGt                >  2:345800/1‑101 (MQ=255)
                                                                                  tCTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGtca              <  2:874989/101‑3 (MQ=255)
                                                                                   cTGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGTCAt             >  2:1491552/1‑101 (MQ=255)
                                                                                    tGGATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGTCAta            >  2:115537/1‑100 (MQ=255)
                                                                                     ggATGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGTCAtac           <  1:1373759/101‑3 (MQ=255)
                                                                                       aTGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGTCATACCg         <  2:855774/101‑1 (MQ=255)
                                                                                       aTGAAAAAGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGTCATACCg         <  1:241240/101‑1 (MQ=255)
                                                                                            aaaGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGTCATACCGGGAAt    <  1:345800/101‑1 (MQ=255)
                                                                                            aaaGAAGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGTCATACCGGGAAt    <  1:1196651/101‑1 (MQ=255)
                                                                                              agaaGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGTCATACCGGGAATGa  <  2:1242356/101‑1 (MQ=255)
                                                                                              agaaGGTAACCACACAACGGCATACCGCGATGGTTCCGGTATATGGACCATCTGTCGTGGTGCCACAATGGTGGATGGTAAGCCCGTCATACCGGGAATGa  <  1:577048/101‑1 (MQ=255)
                                                                                                 |                                                                                                 
TAAAAATCAGATACAGCCTGTCGGCCGCTGTTCTGGCACTGATTGCCGCCAGTGCTCCTGCCCCTGACATTCTCGATCAGTTTCTGGATGAAAAAGAGGGGAATCACACCACTGCATACCGCGATGGCTCCGGCATATGGACCATCTGTCGGGGAGCAACGATGGTGGACGGTAAACCTGTTATTCCGGGAATGA  >  178/959‑1153

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