Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 454,926 C→T 9.7% G190G (GGC→GGT tig → peptidyl‑prolyl cis/trans isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb454,9260CT 9.7% 114.5 / 2.1 41G190G (GGC→GGTtigpeptidyl‑prolyl cis/trans isomerase
Reads supporting (aligned to +/- strand):  ref base C (14/23);  new base T (2/2);  total (16/25)
Fisher's exact test for biased strand distribution p-value = 6.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

CTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATC  >  W3110S.gb/454857‑454995
                                                                     |                                                                     
cTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGcc                                                                      <  1:1455300/71‑1 (MQ=255)
   cACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCa                                                                     <  1:1302839/69‑1 (MQ=255)
       ggTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCg                                                               <  1:3027453/71‑1 (MQ=255)
       ggTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCg                                                               <  1:1540825/71‑1 (MQ=255)
         ttCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCATCTGATTTCGTACTGGCGATGGGCCAGGGTCGTa                                                             >  1:3269932/1‑71 (MQ=255)
          tCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAg                                                                    >  1:2709148/1‑63 (MQ=255)
           cTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTAt                                                            <  1:3149465/70‑1 (MQ=255)
           cTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTAt                                                            >  1:3332298/1‑70 (MQ=255)
           cTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTAt                                                            <  1:1503854/70‑1 (MQ=255)
              tAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGc                                                                       <  1:2023710/56‑1 (MQ=255)
                gACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGTCAGGGTCGTATGCTcc                                                       >  1:3703470/1‑70 (MQ=255)
                gACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGTCAGGGTCGTATGATccc                                                      >  1:2835870/1‑71 (MQ=255)
                gACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTAt                                                            >  1:2153593/1‑65 (MQ=255)
                gACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATccc                                                      <  1:1469700/71‑1 (MQ=255)
                 aCGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCTATGGGTCAGGGTCGTATGATCCCg                                                     <  1:561653/71‑1 (MQ=255)
                 aCGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGTCAGGGTCGTATGATCCCg                                                     <  1:3469991/71‑1 (MQ=255)
                       aaGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCttt                                               >  1:3037574/1‑71 (MQ=255)
                        agagTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTg                                              <  1:3565122/71‑1 (MQ=255)
                         gagTTCGAAGGCGGTAAAGCGTCTGATTTGGTACTGGCGATGGGCCAGGGTCg                                                               <  1:139026/53‑1 (MQ=255)
                          agTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCg                                                               <  1:884931/52‑1 (MQ=255)
                           gTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGaa                                            <  1:403566/70‑1 (MQ=255)
                            ttCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGcc                                                                      <  1:3358033/43‑1 (MQ=255)
                            ttCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGcc                                                                      <  1:1202051/43‑1 (MQ=255)
                              cGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACg                                        >  1:843522/1‑71 (MQ=255)
                                      gTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCaa                                 <  1:1037026/70‑1 (MQ=255)
                                         aaGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGc                             >  1:2654491/1‑71 (MQ=255)
                                            cGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCcac                          <  1:1724349/71‑1 (MQ=255)
                                             gTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCgg                                                    >  1:1878043/1‑44 (MQ=255)
                                             gTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCcaca                         <  1:1080338/71‑1 (MQ=255)
                                               cTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCaaa                                >  1:1626884/1‑62 (MQ=255)
                                                  aTTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTg                                              <  1:2585269/45‑1 (MQ=255)
                                                  aTTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAg                               >  1:1766813/1‑60 (MQ=255)
                                                  aTTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCt                    <  1:2258332/71‑1 (MQ=255)
                                                   tttCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACgg                                       <  1:3572725/51‑1 (MQ=255)
                                                     tCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGc                 >  1:2852663/1‑71 (MQ=255)
                                                       gTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCg                <  1:2585238/70‑1 (MQ=255)
                                                           tGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTg                                              <  1:1111366/36‑1 (MQ=255)
                                                           tGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATc                                   <  1:3543929/47‑1 (MQ=255)
                                                           tGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACaaa                       <  1:1998321/59‑1 (MQ=255)
                                                             gCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAag             >  1:1735077/1‑67 (MQ=255)
                                                             gCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGt          >  1:2004207/1‑70 (MQ=255)
                                                              cGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTc        <  1:3483577/71‑1 (MQ=255)
                                                                   ggCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCaaa                                >  1:1644084/1‑42 (MQ=255)
                                                                    gCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACaaa                       <  1:36900/50‑1 (MQ=255)
                                                                    gCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATc  <  1:1437612/71‑1 (MQ=255)
                                                                     |                                                                     
CTTCACCGGTTCTGTAGACGGCGAAGAGTTCGAAGGCGGTAAAGCGTCTGATTTCGTACTGGCGATGGGCCAGGGTCGTATGATCCCGGGCTTTGAAGACGGTATCAAAGGCCACAAAGCTGGCGAAGAGTTCACCATC  >  W3110S.gb/454857‑454995

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: