Predicted mutation
evidence seq id position mutation annotation gene description
RA REL606 30,364 G→C A6P (GCC→CCC)  ispH → 4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*REL60630,3640GC97.6% 160.1 / NA 42A6P (GCC→CCC) ispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base C (20/21);  minor base T (0/1);  total (20/22)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

GACCGTTCATTTTGATATTGAAGTGCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGGCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTAC  >  REL606/30293‑30435
                                                                       |                                                                       
gACCGTTCATTTTGATATTGAAGTGCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGTTCCTGTTGCCCaa                                                                     <  1:77444/76‑1 (MQ=255)
gACCGTTCATTTTGATATTGAAGTGCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCaa                                                                     <  1:190366/76‑1 (MQ=255)
     ttCATTTTGATATTGAAGTGCTGGAAACCGATCCGGCACTGGAGNCGNAACATGCAGATCCTNTTGCCCAACCcgc                                                                >  2:1085138/1‑76 (MQ=255)
       cATTTTGATATTGAAGTGCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGt                                                              >  1:291140/1‑76 (MQ=255)
             gATATTGAAGTGCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGtttt                                                        <  1:1031688/76‑1 (MQ=255)
             gATATTGAAGTGCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGtttt                                                        <  1:776860/76‑1 (MQ=255)
             gATATTGAAGTGCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATACTGTTGCCCANCCCGCGTGGtgtg                                                        >  1:152083/1‑73 (MQ=255)
               tatTGAAGTGCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTtg                                                      >  2:519681/1‑76 (MQ=255)
                        gCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGt                                             <  1:596304/76‑1 (MQ=255)
                        gCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTNTTGTGCCGGGGt                                             <  1:11297/76‑1 (MQ=255)
                          tGGAAATCNAGCCGGCACTNGAGGCGTAACATGCAGATCCTGNTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAg                                           <  2:1008266/76‑1 (MQ=255)
                          tGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAg                                           <  1:708191/76‑1 (MQ=255)
                           ggAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGa                                          >  1:307737/1‑76 (MQ=255)
                             aaaTCGATCCGGCACTGGAGGCGTAACATGCAGNTCCTGTTGCCCAACCCGCGTGGNTTTTGTGCAGNGGTAGAcc                                        >  2:1018041/1‑76 (MQ=255)
                                     ccGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATc                                >  1:234629/1‑76 (MQ=255)
                                     ccGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATc                                <  1:928654/76‑1 (MQ=255)
                                      cGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGAGGTAGACCGCGCTATCa                               >  1:136920/1‑76 (MQ=255)
                                      cGGCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGCGGCGGGGTAGACCGCGCTATCa                               <  1:6069/76‑1 (MQ=255)
                                        gCACTGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGc                             <  1:187969/76‑1 (MQ=255)
                                          aCTNGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGNTTTTGTGCCGGGGTAGACCGCGCTATCAGCAt                           <  2:1133926/76‑1 (MQ=255)
                                          aCTGGAGGCGTAACATGCAGATCCTGTTGCCCCACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTCTCAGCAt                           >  2:546719/1‑76 (MQ=255)
                                            tGGAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCAttg                         >  2:850345/1‑76 (MQ=255)
                                             ggAGGCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCAttgt                        >  1:855936/1‑76 (MQ=255)
                                                ggCGTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGa                     <  2:833256/76‑1 (MQ=255)
                                                   gTAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGAACGCGCTATCCGCATTGTTGaaaa                  >  1:1039274/1‑76 (MQ=255)
                                                    tAACATGCAGATCCTGTTGCCCGACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAAc                 >  1:310666/1‑76 (MQ=255)
                                                    tAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAAc                 >  2:465195/1‑76 (MQ=255)
                                                    tAACATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAAc                 >  2:839646/1‑76 (MQ=255)
                                                     aaCATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAAcg                <  1:844759/76‑1 (MQ=255)
                                                     aaCATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAAcg                <  1:786242/76‑1 (MQ=255)
                                                       cATGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAAcgcg              >  1:267206/1‑76 (MQ=255)
                                                         tGCAGATCCTGTTGTCCAACCCGCCTGGGTTNTGTGCCGGTGTAGANCGCGCNATCAGCATTGTTGTAAACGCGCt            <  2:151578/76‑1 (MQ=255)
                                                         tGCAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCt            >  1:682932/1‑76 (MQ=255)
                                                           cAGATCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTgg          >  2:398140/1‑76 (MQ=255)
                                                               tCCTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCAt      <  2:11238/76‑1 (MQ=255)
                                                                ccTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCAtt     >  2:122529/1‑76 (MQ=255)
                                                                ccTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCAtt     <  2:359530/76‑1 (MQ=255)
                                                                 cTGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCAttt    <  2:421977/76‑1 (MQ=255)
                                                                  tGTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTa   <  1:805944/76‑1 (MQ=255)
                                                                  gtttGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTa   <  1:969223/74‑1 (MQ=255)
                                                                   gTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTAc  >  1:163888/1‑76 (MQ=255)
                                                                   gTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTAc  <  2:749039/76‑1 (MQ=255)
                                                                   gTTGCCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTAc  <  1:943678/76‑1 (MQ=255)
                                                                       |                                                                       
GACCGTTCATTTTGATATTGAAGTGCTGGAAATCGATCCGGCACTGGAGGCGTAACATGCAGATCCTGTTGGCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTAC  >  REL606/30293‑30435

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