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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1911,8220CT96.5% 59.2 / ‑3.7 29intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (10/18);  total (10/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.17e-01

AGCAGTTGCTGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCGCTCTCGCTTTCGGCTGCTGCATTCTGGCTATTTTTCGCCGCAGTTTCACTGGCTTTGGCATTCGTTTCGCTGGTTTTCGCTGCCGTCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCTATTTTTCGCGTTGTTT  >  191/1724‑1919
                                                                                                  |                                                                                                 
agcagTTGCTGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGttt                                                                                                 <  1:1357832/101‑1 (MQ=17)
   agTTGCTGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGGACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTgg                                                                                              >  2:569838/1‑101 (MQ=17)
      tGCTGATCCAGCTGCTGAAGTGGCAGACCCGGCTGCCGCACTCTCGCTTTCGGCTTCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCAt                                                                                           <  2:163311/101‑1 (MQ=14)
      tGCTGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCAt                                                                                           <  2:485577/101‑1 (MQ=14)
       gCTGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCAtt                                                                                          <  2:1362197/101‑1 (MQ=14)
        cTGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTc                                                                                         >  1:567270/1‑101 (MQ=17)
        cTGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTc                                                                                         >  2:716364/1‑101 (MQ=14)
        cTGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTc                                                                                         <  1:107575/101‑1 (MQ=14)
         tGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCg                                                                                        >  1:360556/1‑101 (MQ=17)
          gATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGt                                                                                       <  1:320482/101‑1 (MQ=14)
           aTCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGtt                                                                                      <  2:1031517/101‑1 (MQ=14)
            tCCAGCTGCTGAAGTCGCAGAACCGGCTGCTGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGttt                                                                                     <  1:638281/101‑1 (MQ=14)
              cAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCg                                                                                   <  2:1544540/101‑1 (MQ=14)
              cAGCTGCTGAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCg                                                                                   <  1:1299320/101‑1 (MQ=14)
                  tgctgAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCttg                                                                               >  1:1089836/1‑99 (MQ=14)
                    ctgAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTTCTGCAGCCTGACTGCTTTTCGCCGCCGTTTCACTGGTTTTGGCATTCGTTTCGCTTGtt                                                                             >  1:787250/1‑101 (MQ=11)
                    ctgAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGtt                                                                             >  1:447397/1‑101 (MQ=11)
                     tgAAGTCGCAGAACCGGCTGCCGCACTCTCGCTTTCGGCTGCTGCAGCCTGACTGCTTTTCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGttt                                                                            <  2:343351/101‑1 (MQ=11)
                                                                   gcctgactgcttttCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTtgc                              <  2:1392122/91‑1 (MQ=11)
                                                                    cctgactgcttttCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTtgct                             >  2:500675/10‑101 (MQ=11)
                                                                             ttttCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCg                    <  2:1613460/101‑1 (MQ=255)
                                                                              tttCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGc                   >  2:1430797/1‑101 (MQ=255)
                                                                               ttCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCt                  <  2:567265/101‑1 (MQ=255)
                                                                               ttCGCCGCAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGATGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCt                  >  2:314650/1‑101 (MQ=17)
                                                                                 cgccgcAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCTAt                <  1:913249/101‑1 (MQ=255)
                                                                                  gccgcAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCTAtt               <  2:792654/101‑1 (MQ=255)
                                                                                  gccgcAGTTTCACTGGCTTTGGCATTCGTTTCGCTGGTTTTCGCTGCCGTCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCGGCTGTCTCGCTGtt               <  2:699982/101‑1 (MQ=255)
                                                                                     gcAGTTTCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCTAttttt            <  2:373090/101‑1 (MQ=255)
                                                                                          ttCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCTATTTTTCGCgt       <  1:1181684/101‑1 (MQ=255)
                                                                                           tCACTGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCTATTTTTCGCgtt      >  1:444403/1‑101 (MQ=21)
                                                                                               tGGTTTTGGCATTCGTTTCGCTTGTTTTCGCTGCCGCCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCTATTTTTCGCGTTGttt  <  2:447392/101‑1 (MQ=255)
                                                                                                  |                                                                                                 
AGCAGTTGCTGATCCAGCTGCTGAAGTCGCAGAACCGGCTGCCGCGCTCTCGCTTTCGGCTGCTGCATTCTGGCTATTTTTCGCCGCAGTTTCACTGGCTTTGGCATTCGTTTCGCTGGTTTTCGCTGCCGTCTGGCTGGACTTTGCGTTGGTTTCGCTCGTCTTTGCTGCTGTCTCGCTATTTTTCGCGTTGTTT  >  191/1724‑1919

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