Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,478,630 T→G 19.5% R261S (AGA→AGC)  pgm ← phosphoglycerate mutase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,478,6300TG19.5% 68.6 / 3.1 43R261S (AGA→AGC) pgmphosphoglycerate mutase
Reads supporting (aligned to +/- strand):  ref base T (10/24);  new base G (7/2);  total (17/26)
Fisher's exact test for biased strand distribution p-value = 1.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTACATAGCCGTCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATCAACAACATCCAGGGC  >  BS168_reference_genome/3478493‑3478767
                                                                                                                                         |                                                                                                                                         
gCTACATAGCCGTCAACTGGTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGcggtccggtc                                                                                                                                          <  8:164695/139‑1 (MQ=255)
   aCATAGCCGTCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCGTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGTTGTTTTTGAACGTGGTGGGAAAATGGGAGGCGcggtccggtcgg                                                                                                                                        >  4:175039/1‑136 (MQ=255)
          cGTCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATtata                                                                                                                                <  3:41908/139‑1 (MQ=255)
            tCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTGTTTGTTTGTGAACGTGTTTGAAATCTGGATGGCGCGGGCCGGGCTGAAATTATAg                                                                                                                               >  6:290556/1‑138 (MQ=255)
            tCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTGGTTGGTTGTTGACGTGGTTGAAAAGTGGCTTGCGCGGGCCGGTCGTGAATTATAg                                                                                                                               >  8:170975/1‑138 (MQ=255)
            tCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGTTGTTTGTGAACGTGTTGGAAATATGGTTGGCGCGGTGCGGTCGGAAATTATAg                                                                                                                               >  5:319465/1‑138 (MQ=255)
              aaCTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTGATTGTTTGTGTACCTGTTGGAAATCAGGCTGGCGCGGTCCGGGGTGAAATTATAGaa                                                                                                                             >  1:336708/1‑138 (MQ=255)
              aaCAGTTTCACTGAAGGGAGTCAGGCAGACGAAATACCGGTTCTTAGGATACTTCTCGCCGCTGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGCTGGCGCGGTCCGGTCTACAATTATAGaa                                                                                                                             <  8:121108/138‑1 (MQ=255)
                cTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGaa                                                                                                                             <  4:142988/136‑1 (MQ=255)
                   tttCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTGGTGGTTTGTGATCGGGTTGGTAATCTGGGTGGGGGGGGCCGGGGGGCAATTATTGTAAATgac                                                                                                                       >  2:302396/1‑139 (MQ=255)
                           aaGTGAGTCAGGCAGACGAAATAGAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCg                                                                                                                <  3:175039/138‑1 (MQ=255)
                             gTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGGTGGTTGTTTTCCTGTTGGCAATCGGGTTTGCGCGGTGCGGGCGGATTTTTTTTAAAATCACAGaaaacccg                                                                                                             >  4:51963/1‑133 (MQ=255)
                                aGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCt                                                                                                          >  1:11854/1‑139 (MQ=255)
                                aGTCAGGCAGACGAAATACAGGTTCGTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTTTTTGTGAACCTGTTGGCAAAATGGATGGCGCGGTCCGGGCTGAAATTATAGAAAATAACAGAATAGCCctct                                                                                                          >  6:330292/1‑139 (MQ=255)
                                     ggCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGAt                                                                                                     >  7:426985/1‑139 (MQ=255)
                                      gCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCCTCGTggtat                                                                                                    >  7:139220/1‑136 (MQ=255)
                                              aaTACTGGTTCTTTGGATAATTCTCGCCGCTGTCGAAGTCACGGAAGTCTTTGTTTGTGAACTTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCaa                                                                                             <  5:51353/138‑1 (MQ=255)
                                                tACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTTGATTTTCGCAACAg                                                                                          >  1:2248/1‑139 (MQ=255)
                                                   aGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCt                                                                                       <  1:175229/139‑1 (MQ=255)
                                                    ggTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGGCTGAAATTATTGAAAATCACAGAATCGCCGTCTTGGTTTTTCGAACCAGGCTc                                                                                      >  7:64730/1‑139 (MQ=255)
                                                     gTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCGTTGTTTTTGTACGTTTTGGAAATCTGGATGGCGCGGGCCGGGCTGAAATTATATAAAAAAACACAATAGCCGTGGTGGTTTTTTGCAACAGGCGCa                                                                                     >  1:405211/1‑139 (MQ=255)
                                                        cTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGTAATCTGGATGGCGCGGTCCGGTCTGAAATTATATAAAATAACAGAATCGCCGTGTTGGTTTTTCTCAACAGGCGCACCg                                                                                  >  3:403069/1‑139 (MQ=255)
                                                           tGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGtt                                                                                <  8:72827/138‑1 (MQ=255)
                                                               tAATTCTCGCCGCGGGCGAAGGCTCGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGGCTTGGATTTTCGCAACAGGCTCACCGTTTTCtt                                                                           <  4:337415/139‑1 (MQ=255)
                                                                aaTTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGGCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCtt                                                                           <  2:336708/138‑1 (MQ=255)
                                                                  ttCTCGCCGCGGTCGAAGTCACTGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTtg                                                                        <  8:53599/139‑1 (MQ=255)
                                                                  ttCTCGCCGCGGTCGAAGTCACGGAAGTGTTTGTTTGTGTACGTGTTTGATATGTGGATGGGGCGGGCCGGTGTGAAATTATAGAAAATCACAGAATCGCCGGGGTGGGTTTTGGCAACAGGGGCACCGGGTTGTTTtg                                                                        >  1:121725/1‑139 (MQ=255)
                                                                                      aCGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCAc                                                     <  2:11854/138‑1 (MQ=255)
                                                                                         gAAGTGTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAAt                                                 <  5:330292/139‑1 (MQ=255)
                                                                                         gAAGTCTTAGTTTGTGAACGTGTTGGAAATCTGGCTGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAAt                                                 <  6:103060/139‑1 (MQ=255)
                                                                                             tCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCAt                                             <  5:430380/139‑1 (MQ=255)
                                                                                             tCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCAt                                             <  6:186754/139‑1 (MQ=255)
                                                                                                         atcGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCATGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAt                                 <  2:133326/137‑1 (MQ=255)
                                                                                                         aaCGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAt                                 <  6:192207/139‑1 (MQ=255)
                                                                                                         aaCGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAt                                 >  5:192207/1‑139 (MQ=255)
                                                                                                         aaCGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAt                                 <  7:170975/139‑1 (MQ=255)
                                                                                                         aaCGTGTTGGAAATCTGGATGGCGCGGGCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCGCACCGTTTTCTTTGGTGATGACAGATGGAATCACGAATTCATCGTAGATACCAt                                 <  1:296562/139‑1 (MQ=255)
                                                                                                                    aaTCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGaa                      <  6:224295/139‑1 (MQ=255)
                                                                                                                          gATGGCGCTGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAatc                    <  4:225752/135‑1 (MQ=255)
                                                                                                                          gATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGTCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAatc                    <  2:142765/135‑1 (MQ=255)
                                                                                                                              tcgcgGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaaca             <  3:51963/137‑1 (MQ=255)
                                                                                                                              gcgcGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaaca             <  8:139220/138‑1 (MQ=255)
                                                                                                                              gcgcGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGAGTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaaca             <  3:290681/138‑1 (MQ=255)
                                                                                                                              gcgcGGTACGGGATTAAATTATAGAAAATCACAGAATCGCCGTCTTTGATGTTCGAAACATGATCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTATATACCATTTGCATATGAATCATcaaca             <  5:207478/138‑1 (MQ=255)
                                                                                                                               cgcgGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaacaa            >  1:213173/1‑138 (MQ=255)
                                                                                                                                    tccggtcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATAGTAGATACCATTTGCATATGAATAATCAACAACATc        >  7:300576/1‑137 (MQ=255)
                                                                                                                                    tccggtcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATATTCGCAACAGGCTCACCGTTTTATTTTTTGATGACAGATGGAATCACGAATTCAAAGTAGATACAAATTTCATATGAATCATCAAACACATACa      >  4:128333/1‑139 (MQ=255)
                                                                                                                                        gtcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATCAACAACATCCAGGgc  >  5:343105/1‑139 (MQ=255)
                                                                                                                                         tcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaac              >  7:337391/1‑126 (MQ=255)
                                                                                                                                         tcTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATcaac              <  8:337391/126‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                         
GCTACATAGCCGTCAACTGTTTCACTGAAGTGAGTCAGGCAGACGAAATACAGGTTCTTTGGATAATTCTCGCCGCGGTCGAAGTCACGGAAGTCTTTGTTTGTGAACGTGTTGGAAATCTGGATGGCGCGGTCCGGTCTGAAATTATAGAAAATCACAGAATCGCCGTCTTGGATTTTCGCAACAGGCTCACCGTTTTCTTTTGTGATGACAGATGGAATCACGAATTCATCGTAGATACCATTTGCATATGAATCATCAACAACATCCAGGGC  >  BS168_reference_genome/3478493‑3478767

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