Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 854,906 A→G 20.0% H28H (CAT→CAC yciG ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported854,9060AG20.0% 105.1 / 22.2 55H28H (CAT→CACyciGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (21/23);  new base G (4/7);  total (25/30)
Fisher's exact test for biased strand distribution p-value = 7.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCAC  >  Exported/854836‑854973
                                                                      |                                                                   
cGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                                                                     <  1:2394448/71‑1 (MQ=255)
  cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATg                                                                   >  1:1814136/1‑71 (MQ=255)
     tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCt                                                                 >  1:599170/1‑70 (MQ=255)
         taccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                            >  1:5201179/2‑70 (MQ=255)
         gACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                                                                     <  1:262660/62‑1 (MQ=255)
          accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTgtc                                                              <  1:2902977/68‑3 (MQ=255)
          accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTccac                                                           >  1:2807254/1‑69 (MQ=255)
          accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                            <  1:1115050/70‑2 (MQ=255)
          accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                            <  1:2408153/70‑2 (MQ=255)
          accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                                            <  1:2414555/70‑2 (MQ=255)
          aCCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGAcctc                                                           >  1:574382/1‑69 (MQ=255)
             gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGTATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                        <  1:3742697/71‑1 (MQ=255)
             gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                        <  1:3456577/71‑1 (MQ=255)
             gCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                        <  1:49739/71‑1 (MQ=255)
              ccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttt                                                       >  1:504233/1‑71 (MQ=255)
               cTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCttt                                                       <  1:1081837/70‑1 (MQ=255)
                ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGc                                                           <  1:3452642/65‑1 (MQ=255)
                ttttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTAc                                                     >  1:226188/1‑71 (MQ=255)
                  ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                                     <  1:2860583/69‑3 (MQ=255)
                  ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                                     <  1:3709759/69‑3 (MQ=255)
                  ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                                     <  1:1200742/69‑3 (MQ=255)
                  ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACgg                                                   >  1:1301709/1‑71 (MQ=255)
                     tACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCtt                                                        <  1:1418166/63‑1 (MQ=255)
                      aCCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCt                                                <  1:2288744/70‑1 (MQ=255)
                         cGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                                                                     >  1:276554/1‑46 (MQ=255)
                            ttCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCt                                                                 >  1:4389781/1‑47 (MQ=255)
                              cAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGAt                                       >  1:3225039/1‑71 (MQ=255)
                              cAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGATATGCTGACCGCCTTTACGGCCTGCGTCGGAt                                       >  1:4576986/1‑71 (MQ=255)
                                 aTGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGcc                                    <  1:2214833/71‑1 (MQ=255)
                                       gTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTtctaa                              <  1:4430476/71‑3 (MQ=255)
                                       gTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCa                              <  1:1577361/71‑1 (MQ=255)
                                            ggATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTc                         <  1:1007233/71‑1 (MQ=255)
                                              aTCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCtt                       <  1:4524946/71‑1 (MQ=255)
                                              aTCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCtt                       <  1:4348405/71‑1 (MQ=255)
                                               tCATTCTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTc                      >  1:4203481/1‑71 (MQ=255)
                                                cATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCg                     >  1:4987754/1‑71 (MQ=255)
                                                cATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCg                     >  1:2239612/1‑71 (MQ=255)
                                                 aTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                                        >  1:122792/1‑51 (MQ=255)
                                                 aTTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCgg                    <  1:2795517/71‑1 (MQ=255)
                                                  tttttAAAATTACCGCCGCTATGCTGACCGCCttt                                                       <  1:1184598/35‑1 (MQ=255)
                                                  tttttAAAATTACCGCCGCTATGCTGACCGCCttt                                                       <  1:1575907/35‑1 (MQ=255)
                                                  tttttAAAATTACCGCCGCTATGCTGACCGCATTTACGGCCTGCGTCGGATGCCTTCTCAc                             <  1:5426690/61‑1 (MQ=255)
                                                   ttttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGc                   >  1:3910092/1‑70 (MQ=255)
                                                    tttAAAATTATCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGa                                        <  1:4904911/48‑1 (MQ=255)
                                                    tttAAAATTACCGCCGCTATGCTGACCTCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCg                  >  1:4801684/1‑70 (MQ=255)
                                                    tttAAAATTACCGCCGCTATGCTGACCGCCTTTAc                                                     >  1:5648954/1‑35 (MQ=255)
                                                    tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACTGTCTTCgg                    >  1:4192220/1‑68 (MQ=255)
                                                    tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:1497894/1‑62 (MQ=255)
                                                    tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:679450/1‑62 (MQ=255)
                                                    tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:1160322/1‑62 (MQ=255)
                                                    tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:136047/1‑62 (MQ=255)
                                                    tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCg                  >  1:70114/1‑70 (MQ=255)
                                                     ttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt                          >  1:1891397/1‑61 (MQ=255)
                                                     ttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGa                 <  1:3303030/70‑1 (MQ=255)
                                                     ttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGa                 <  1:1108846/70‑1 (MQ=255)
                                                                gccgcTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAac     <  1:5376772/71‑1 (MQ=255)
                                                                  cgcTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAacca   <  1:1263127/71‑1 (MQ=255)
                                                                   gcTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAaccac  <  1:2010800/71‑1 (MQ=255)
                                                                      aTGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTc           <  1:2353036/59‑1 (MQ=255)
                                                                      aTGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTg        >  1:1894013/1‑62 (MQ=255)
                                                                      aTGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTg        >  1:2560812/1‑62 (MQ=255)
                                                                      |                                                                   
CGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCGAAATTTCCTGAACCAC  >  Exported/854836‑854973

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: