Predicted mutation
evidence seq id position mutation freq annotation gene description
RA REL606 1,778,501 T→C 12.4% E446G (GAA→GGA)  thrS ← threonyl‑tRNA synthetase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*REL6061,778,5010T→C12.4% 2.7 ‑3.9 ‑0.0391337thrSthreonyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  new base (0/5):  ref base (20/14):  total (20/19)
Fisher's exact test strand distribution p-value = 2.02e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 9.92e-01

ACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCC  >  REL606/1778467‑1778536
                                  |                                   
aCCCAGTTGATATTCAAACGGGATGTTGTTttcttc                                    <  1:1007300/36‑1 (MQ=255)
aCCCAGTTGATATTCAAACGGGATGTTGTTttcttc                                    >  1:1743690/1‑36 (MQ=255)
aCCCAGTTGATATTCAAACGGGATGTTGTTttcttc                                    >  1:4674281/1‑36 (MQ=255)
aCCCAGTTGATATTCAAACGGGATGTTGTTttcttc                                    >  1:4415704/1‑36 (MQ=255)
 cccAGTTGATATTCAAACGGGATGTTGTTTTCTTcc                                   <  1:5037207/36‑1 (MQ=255)
 cccAGTTGATATTCAAACGGGATGTTGTTTTCTTcc                                   <  1:4303215/36‑1 (MQ=255)
  ccAGTTGATATTCAAACGGGATGTTGTTTTCTTCCa                                  <  1:7617581/36‑1 (MQ=255)
   cAGTTGATATTCAAACGGGATGTTGTTTTCTTCTAg                                 >  1:8306675/1‑36 (MQ=255)
    aGTTGATATTCAAACGGGATGTTGTTTTCTTCCAgc                                >  1:7081995/1‑36 (MQ=255)
     gTTGATATTCAAACGGGATGTTGTTTTCTTCCAgcg                               >  1:7280635/1‑36 (MQ=255)
      ttGATATTCAAACGGGATGTTGTTTTCTTCCAgcgc                              <  1:2045162/36‑1 (MQ=255)
       tGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCa                             >  1:6859168/1‑36 (MQ=255)
       tGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCa                             <  1:6345788/36‑1 (MQ=255)
       tGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCa                             >  1:3828872/1‑36 (MQ=255)
       tGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCa                             >  1:3842977/1‑36 (MQ=255)
        gATATTCAAACGGGATGTTGTTTTCTTCCAGCGCaa                            <  1:6169261/36‑1 (MQ=255)
        gATATTCAAACGGGATGTTGTTTTCTTCCAGCGCaa                            <  1:3742721/36‑1 (MQ=255)
        gATATTCAAACGGGATGTTGTTTTCTTCCAGCGCaa                            >  1:2627625/1‑36 (MQ=255)
        gATATTCAAACGGGATGTTGTTTTCTTCCAGCGCaa                            >  1:2192301/1‑36 (MQ=255)
         atatTCAAACGGGATGTTGTTTTCTTCCAGCGCAAc                           <  1:4787966/36‑1 (MQ=255)
         atatTCAAACGGGATGTTGTTTTCTTCCAGCGCAAc                           <  1:8276545/36‑1 (MQ=255)
         atatTCAAACGGGATGTTGTTTTCTTCCAGCGCAAc                           >  1:5816017/1‑36 (MQ=255)
         atatTCAAACGGGATGTTGTTTTCTTCCAGCGCAAc                           >  1:5441296/1‑36 (MQ=255)
          tatTCAAGCGGGATGTTGTTTTCTTCCAGCGCAAcc                          <  1:5890334/36‑1 (MQ=255)
          tatTCAAACGGGATGTTGTTTTCTTCCAGCGCAAcc                          >  1:3011076/1‑36 (MQ=255)
          tatTCAAACGGGATGTTGTTTTCTTCCAGCGCAAcc                          <  1:6502620/36‑1 (MQ=255)
           atTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCg                         <  1:4799527/36‑1 (MQ=255)
            ttCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGc                        >  1:5224299/1‑36 (MQ=255)
                aaCGGGATGTTGTTTTCTTCCAGCGCAACCGCCAgg                    >  1:8322584/1‑36 (MQ=255)
                 aCGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGt                   <  1:2871850/36‑1 (MQ=255)
                 aCGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGt                   <  1:8176909/36‑1 (MQ=255)
                    ggATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCg                >  1:5625111/1‑36 (MQ=255)
                    ggATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCg                >  1:5881701/1‑36 (MQ=255)
                    ggATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCg                >  1:1688183/1‑36 (MQ=255)
                     gATTTTTTTTTCCCCCCGCGCAACCGCCAGGTCCGc               <  1:1790657/36‑1 (MQ=16)
                        gttgttTTCTTCCAGCGCAACCGCCAGGTCCGcccc            >  1:8103822/1‑34 (MQ=255)
                        gttgttTTCTTCCAGCGCAACCGCCAGGTCCGCCTc            >  1:5257432/1‑36 (MQ=255)
                         tttttttCCCCCACCACACCCGCCAGGTCCGCCTCa           <  1:4857217/33‑1 (MQ=12)
                         ttgttTTCACCCATCGCAACCGCCAGGTCCGCCTCa           <  1:6614529/36‑1 (MQ=37)
                          tgttTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAg          >  1:5480093/1‑36 (MQ=255)
                            ttttCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCa        >  1:7789835/1‑36 (MQ=255)
                             tttCTCCCACCGCAACCGCCAGGTCCGCCTCAGCAc       <  1:185065/36‑1 (MQ=37)
                             tttCTCCCACCCCAACCGCCAGGTCCGCCTCAGCAc       <  1:2539798/36‑1 (MQ=25)
                              ttctccCACCGCACCCGCCAGGTCCGCCTCAGCACg      <  1:7957277/36‑1 (MQ=25)
                               tcttcCAGCGCAACCGCCAGGTCCGCCTCAGCACGa     >  1:6967819/1‑36 (MQ=255)
                               tcttcCAGCGCAACCGCCAGGTACGCCCCAGCACGa     >  1:1995898/1‑36 (MQ=37)
                                cttcCAGCGCAACCGCCAGGTCCGCCTCAGCACGAt    >  1:5161201/1‑36 (MQ=255)
                                 ttcCAGCGCAACCGCCAGGTCCGCCTCAGCACGATc   >  1:5904663/1‑36 (MQ=255)
                                  tcCAGCGCAACCGCCAGGTCCGCCTCAGCACGATcc  >  1:2499022/1‑36 (MQ=255)
                                  |                                   
ACCCAGTTGATATTCAAACGGGATGTTGTTTTCTTCCAGCGCAACCGCCAGGTCCGCCTCAGCACGATCC  >  REL606/1778467‑1778536

Base quality scores:  ATCG  < 3 ≤  ATCG  < 11 ≤  ATCG  < 24 ≤  ATCG  < 32 ≤  ATCG  < 37 ≤  ATCG