Predicted mutation
evidence seq id position mutation annotation gene description
RA 10 6,234 T→C intergenic (–/–)  /  –/–

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*106,2340TC98.2% 191.7 / ‑4.3 58intergenic (–/–)–/––/–
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (36/20);  total (38/20)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.44e-01

TATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGA  >  10/6136‑6329
                                                                                                  |                                                                                               
tataCGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCa                                                                                               >  1:213429/1‑101 (MQ=255)
 ataCGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCaa                                                                                              >  2:752146/1‑101 (MQ=255)
   aCGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaac                                                                                            <  1:1671397/101‑1 (MQ=255)
      cTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaacaac                                                                                         >  1:765891/1‑101 (MQ=255)
      cTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAaacaac                                                                                         <  1:1315015/101‑1 (MQ=255)
        ttGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTc                                                                                       >  2:783539/1‑101 (MQ=255)
        ttGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTc                                                                                       >  1:1410459/1‑101 (MQ=255)
          gATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAc                                                                                     <  2:1193095/101‑1 (MQ=255)
            tCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAccc                                                                                   <  2:1160668/101‑1 (MQ=255)
            tCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCAccc                                                                                   <  2:721891/101‑1 (MQ=255)
                cacaGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACgg                                                                               >  2:113211/1‑101 (MQ=255)
                    gCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAAt                                                                           <  2:993406/101‑1 (MQ=255)
                     cATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATa                                                                          <  1:1026985/101‑1 (MQ=255)
                     cATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATa                                                                          >  2:941254/1‑101 (MQ=255)
                      aTGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATaa                                                                         >  2:1634393/1‑101 (MQ=255)
                      aTGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGACAAACAACTCACCCACGGCAATaa                                                                         >  2:1545367/1‑101 (MQ=255)
                        gAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAAcc                                                                       <  1:942770/101‑1 (MQ=255)
                         aaTGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCt                                                                      >  2:910851/1‑101 (MQ=255)
                          aTGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTg                                                                     <  2:528698/101‑1 (MQ=255)
                              ggTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAg                                                                 >  1:285534/1‑101 (MQ=255)
                               gTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGt                                                                >  1:467185/1‑101 (MQ=255)
                                    ccGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTg                                                           >  2:233539/1‑101 (MQ=255)
                                     cGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAGCTCGCCCACGGCAATAACCTGAAAGTGGCTGc                                                          <  1:752151/101‑1 (MQ=255)
                                         atatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATc                                                      >  2:316814/1‑101 (MQ=255)
                                         atatatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATc                                                      <  1:674201/101‑1 (MQ=255)
                                            tatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGcc                                                   >  1:454832/1‑101 (MQ=255)
                                            tatTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGcc                                                   <  1:225958/101‑1 (MQ=255)
                                              ttCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCgg                                                 >  1:752470/1‑101 (MQ=255)
                                                ccATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGGCATCGCCGGaa                                               >  1:1641122/1‑101 (MQ=255)
                                                ccATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGaa                                               >  2:1137613/1‑101 (MQ=255)
                                                 cATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGTCTGCCATCGCCGGaaa                                              >  2:1192941/1‑101 (MQ=255)
                                                 cATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAAc                                             <  1:596050/101‑1 (MQ=255)
                                                 cATCAGCGAACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGaaa                                              >  2:1535345/1‑101 (MQ=255)
                                                  aTCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAAc                                             >  1:397581/1‑101 (MQ=255)
                                                     aGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTg                                          <  1:484037/101‑1 (MQ=255)
                                                     aGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTg                                          >  1:193507/1‑101 (MQ=255)
                                                       cGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGa                                        >  1:700383/1‑101 (MQ=255)
                                                        ggACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGAc                                       <  1:1634399/101‑1 (MQ=255)
                                                         gACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACt                                      <  2:1531194/101‑1 (MQ=255)
                                                         gACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACt                                      >  2:1492648/1‑101 (MQ=255)
                                                              tAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTg                                 <  2:1410453/101‑1 (MQ=255)
                                                                 aCCGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggag                              <  2:710731/101‑1 (MQ=255)
                                                                   cGTCTGCTGTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggaggg                            >  2:26465/1‑101 (MQ=255)
                                                                     tCTGCTGTTTGTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTggagggag                          >  2:1135775/1‑101 (MQ=255)
                                                                      ctgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgaga                         >  2:1628132/1‑101 (MQ=255)
                                                                       tgctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgagag                        >  2:1096519/1‑101 (MQ=255)
                                                                        gctgTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGgagaga                       >  1:920079/1‑101 (MQ=255)
                                                                           gTTTTTTAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCg                    >  1:257318/1‑101 (MQ=255)
                                                                              ttttAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtc                 <  1:1635849/101‑1 (MQ=255)
                                                                                 tAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtcatc              >  1:1468263/1‑101 (MQ=255)
                                                                                 tAACCCGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTtaatc              >  1:1667872/1‑101 (MQ=255)
                                                                                      cGACTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAacac         >  1:1421197/1‑101 (MQ=255)
                                                                                      cGACGCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAacac         >  2:1718365/1‑101 (MQ=255)
                                                                                         cTCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGct      <  1:650253/101‑1 (MQ=255)
                                                                                          tCATGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGctc     >  1:1151934/1‑101 (MQ=255)
                                                                                            aTGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTg   >  2:455911/1‑101 (MQ=255)
                                                                                            aTGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTg   <  1:316814/101‑1 (MQ=255)
                                                                                             tGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGa  >  1:47578/1‑101 (MQ=255)
                                                                                             tGCCGCCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGa  >  2:517604/1‑101 (MQ=255)
                                                                                                  |                                                                                               
TATACGCTTTGATCGACACAGCATGAATGCGGTTATCCGCAATATATTCCATCAGCGGACCATAGACCGTCTGCTGTTTTTTAACCCGACTCATGCCGTCAAACAACTCACCCACGGCAATAACCTGAAAGTGGCTGCCATCGCCGGAAACGTGGACTTCCTGGAGGGAGAGAGCGTTCATCAACACGCTCTGA  >  10/6136‑6329

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