Predicted mutation
evidence seq id position mutation freq annotation gene description
RA Exported 854,888 T→C 26.3% K34K (AAA→AAG yciG ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*Exported854,8880TC26.3% 84.3 / 33.6 57K34K (AAA→AAGyciGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (22/20);  new base C (7/8);  total (29/28)
Fisher's exact test for biased strand distribution p-value = 7.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

TGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCG  >  Exported/854818‑854957
                                                                      |                                                                     
tGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCAttt                                                                       >  1:1673781/1‑71 (MQ=255)
    tGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                                                    >  1:1199558/1‑70 (MQ=255)
      aTTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaa                                                                   <  1:4841389/69‑1 (MQ=255)
         tattaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                              >  1:4708965/1‑71 (MQ=255)
           ttaCCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                             <  1:4865645/70‑1 (MQ=255)
              ccaccGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAccgc                                                          <  1:1450435/70‑1 (MQ=255)
                accGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACcgccgc                                                       <  1:4981776/71‑1 (MQ=255)
                  cGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTa                                                     <  1:2394448/71‑1 (MQ=255)
                    cTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATg                                                   >  1:1814136/1‑71 (MQ=255)
                       tgttgACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCt                                                 >  1:599170/1‑70 (MQ=255)
                          ttaccaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTAccccc                                                         >  1:2537712/3‑59 (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:2414555/70‑2 (MQ=255)
                            accaccTTTTTTACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTcca                                            <  1:2408153/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)
                                   tttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAAt                                                                <  1:3374548/43‑1 (MQ=255)
                                    ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTa                                                              >  1:223896/1‑44 (MQ=255)
                                    ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                     <  1:3709759/69‑3 (MQ=255)
                                    ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                     <  1:1200742/69‑3 (MQ=255)
                                    ttttACCTGCTTCAGATGCGCGCTGCGGGTCATTCTTGAAATTACCCCCGCTGTGCTGTCCACCttttt                                     <  1:2860583/69‑3 (MQ=255)
                                    ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTaaaa                                                                 >  1:2976150/1‑41 (MQ=255)
                                    ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                                                    >  1:3626140/1‑38 (MQ=255)
                                    ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTa                                                                    >  1:2882276/1‑38 (MQ=255)
                                    ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                              >  1:4340307/1‑44 (MQ=255)
                                    ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                              >  1:3777310/1‑44 (MQ=255)
                                    ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                              >  1:4627231/1‑44 (MQ=255)
                                    ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTa                                                              >  1:4983566/1‑44 (MQ=255)
                                    ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                             >  1:1949034/1‑45 (MQ=255)
                                    ttttACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTAc                                                             >  1:4608047/1‑45 (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:4348405/71‑1 (MQ=255)
                                                                aTCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCtt       <  1:4524946/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:136047/1‑62 (MQ=255)
                                                                      tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt          >  1:1497894/1‑62 (MQ=255)
                                                                      tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt          >  1:1160322/1‑62 (MQ=255)
                                                                      tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGt          >  1:679450/1‑62 (MQ=255)
                                                                      tttAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCg  >  1:70114/1‑70 (MQ=255)
                                                                      |                                                                     
TGCCTGATTTATTACCACCGCTTTGTTGACCGCCTTTTTTACCCGCTTCAGATGCGCGTTGCGGATCATTTTTAAAATTACCGCCGCTATGCTGACCGCCTTTACGGCCTGCGTCGGATGCCTTCTCACGGTCTTCGGCG  >  Exported/854818‑854957

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: