Predicted mutation
evidence seq id position mutation annotation gene description
RA 62 8,758 C→A intergenic (–/–)  /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*628,7580CA98.5% 226.5 / ‑3.9 67intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (38/28);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 4.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.35e-01

CCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGCGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAA  >  62/8662‑8855
                                                                                                |                                                                                                 
ccTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtgg                                                                                               <  1:1015476/101‑1 (MQ=255)
  tGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAgtggtg                                                                                             <  2:1671973/101‑1 (MQ=255)
    gTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTgcg                                                                                           >  2:523057/1‑101 (MQ=255)
      cACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCa                                                                                         >  1:569180/1‑101 (MQ=255)
           ggAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGc                                                                                    >  2:1113517/1‑101 (MQ=255)
              aaTTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACt                                                                                 <  1:1219635/101‑1 (MQ=255)
                ttGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTAt                                                                               <  1:456541/101‑1 (MQ=255)
                 tGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTAtt                                                                              <  1:1124267/101‑1 (MQ=255)
                   aaCTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTAc                                                                            <  2:646639/101‑1 (MQ=255)
                   aaCTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTAc                                                                            <  2:675421/101‑1 (MQ=255)
                   aaCTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTAc                                                                            <  1:512047/101‑1 (MQ=255)
                    aCTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACg                                                                           <  1:606439/101‑1 (MQ=255)
                     cTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGa                                                                          >  2:1188151/1‑101 (MQ=255)
                       ggCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGAcc                                                                        >  2:1350908/1‑101 (MQ=255)
                         cACTAATGCGCCTGTATGAGGTGACTGAAGAGCAGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCaa                                                                      >  2:1059474/1‑101 (MQ=255)
                          aCTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAg                                                                     <  2:1550728/101‑1 (MQ=255)
                            tAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGaa                                                                   >  1:1449481/1‑101 (MQ=255)
                               tGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAatat                                                                <  1:303974/101‑1 (MQ=255)
                                 cgcCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGa                                                              <  2:548210/101‑1 (MQ=255)
                                  gcCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGaa                                                             <  1:143679/101‑1 (MQ=255)
                                   gcTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGaaa                                                            <  1:1071028/100‑1 (MQ=255)
                                    cTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGaaaa                                                           <  2:569175/101‑1 (MQ=255)
                                    cTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGAGGTGCGCAACCGCACTATTACGACCAAGAATATGaaaa                                                           >  2:337437/1‑101 (MQ=255)
                                      tgtaggaGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAgc                                                         <  1:951086/95‑1 (MQ=255)
                                         tGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGcgcg                                                      >  1:14690/1‑101 (MQ=255)
                                         tGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGcgcg                                                      >  1:1586937/1‑101 (MQ=255)
                                           aaGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCggg                                                    >  2:217593/1‑101 (MQ=255)
                                            aGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGc                                                   >  2:1132533/1‑101 (MQ=255)
                                              tGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAg                                                 <  1:415063/101‑1 (MQ=255)
                                               gACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGa                                                <  2:1278224/101‑1 (MQ=255)
                                                aCTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGAc                                               >  2:1633439/1‑101 (MQ=255)
                                                 cTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACa                                              >  1:247718/1‑101 (MQ=255)
                                                  tGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACCGAGAGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACAt                                             >  2:577081/1‑101 (MQ=255)
                                                  tGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACAt                                             <  2:449416/101‑1 (MQ=255)
                                                  tGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGCGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACAt                                             >  1:1580347/1‑101 (MQ=255)
                                                     agagCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGc                                          >  1:1145284/1‑101 (MQ=255)
                                                     agagCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGc                                          <  1:350651/101‑1 (MQ=255)
                                                      gagCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCa                                         <  2:1015961/101‑1 (MQ=255)
                                                       agCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGCGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCAc                                        <  2:892244/101‑1 (MQ=255)
                                                       agCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCAc                                        >  1:1049036/1‑101 (MQ=255)
                                                       agCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCAc                                        <  2:602431/101‑1 (MQ=255)
                                                        gCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACt                                       >  1:1518895/1‑101 (MQ=255)
                                                          cgcgCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTgg                                     >  2:1471675/1‑101 (MQ=255)
                                                          cgcgCTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTgg                                     >  1:1475490/1‑101 (MQ=255)
                                                            cgcTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGca                                   >  1:1155784/1‑101 (MQ=255)
                                                            cgcTACCTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGca                                   >  1:939821/1‑101 (MQ=255)
                                                                 ccTGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCt                              <  2:440696/101‑1 (MQ=255)
                                                                  cTGGCGCTGACGAGCTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTa                             >  1:1352099/1‑101 (MQ=255)
                                                                   tGGCGCTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTAc                            <  1:1025074/101‑1 (MQ=255)
                                                                      cgcTGACGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGc                         <  1:1188156/101‑1 (MQ=255)
                                                                      cgcTGACGAACTATTTTGTCGAACAGAGGGGGGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGc                         >  2:366769/1‑101 (MQ=255)
                                                                           aCGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGc                    >  2:369523/1‑101 (MQ=255)
                                                                            cGAACTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCg                   >  2:1277432/1‑101 (MQ=255)
                                                                            cGAACTATTTTGTCGAACAGAGTGGGGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTCAGGCCCGGCg                   >  2:1323023/1‑101 (MQ=255)
                                                                              aaCTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCgcg                 >  2:1449511/1‑99 (MQ=255)
                                                                               aCTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTgg                >  1:443448/1‑101 (MQ=255)
                                                                                cTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGa               >  1:1618714/1‑101 (MQ=255)
                                                                                cTATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGa               >  2:1134217/1‑101 (MQ=255)
                                                                                 tATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGAt              <  1:217593/101‑1 (MQ=255)
                                                                                 tATTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGAt              >  2:832704/1‑101 (MQ=255)
                                                                                  aTTTTGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATg             >  2:275374/1‑101 (MQ=255)
                                                                                   ttttGTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATgg            <  2:627530/101‑1 (MQ=255)
                                                                                       gTCGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGaa        <  2:1449475/101‑1 (MQ=255)
                                                                                         cGAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCAGGCGCGGATGGGGaaat      >  2:265973/1‑100 (MQ=255)
                                                                                          gAACAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGa     >  1:818216/1‑101 (MQ=255)
                                                                                           aaCAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGAc    >  2:1677794/1‑101 (MQ=255)
                                                                                             cAGAGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACaa  >  2:736408/1‑101 (MQ=255)
                                                                                                |                                                                                                 
CCTGGTCACCCGGAAATTGAACTGGCACTAATGCGCCTGTATGAAGTGACTGAAGAGCCGCGCTACCTGGCGCTGACGAACTATTTTGTCGAACAGCGTGGTGCGCAACCGCACTATTACGACCAAGAATATGAAAAGCGCGGGCAGACATCGCACTGGCACACCTACGGCCCGGCGTGGATGGTGAAAGACAA  >  62/8662‑8855

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