Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,903,549 T→C 19.1% S104P (TCA→CCA)  bax → conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,903,5490TC19.1% 147.8 / 30.9 68S104P (TCA→CCA) baxconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/16);  new base C (6/7);  total (45/23)
Fisher's exact test for biased strand distribution p-value = 1.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.39e-01

TCTCCGGACCGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTACA  >  W3110S.gb/3903479‑3903618
                                                                      |                                                                     
tctcCGGACCGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCAttt                                                                       >  1:3098023/1‑71 (MQ=255)
tctcCGGACCGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCAttt                                                                       >  1:2219501/1‑71 (MQ=255)
    cGGACCGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCaa                                                                    <  1:3058719/70‑1 (MQ=255)
    cGGACCGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCaa                                                                    <  1:2396864/70‑1 (MQ=255)
           tAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTAt                                                            <  1:3115036/71‑1 (MQ=255)
           tAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCAAAACAGt                                                              <  1:1319605/69‑1 (MQ=255)
           tAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCAAAACAGt                                                              <  1:1408987/69‑1 (MQ=255)
             aTGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTAt                                                            >  1:2276262/1‑69 (MQ=255)
             aTGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTAt                                                            >  1:3767785/1‑69 (MQ=255)
             aTGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTAt                                                            >  1:1868882/1‑69 (MQ=255)
             aTGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTAt                                                            >  1:1332092/1‑69 (MQ=255)
             aTGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGATCATTTCAAAACAGTAt                                                            >  1:2035323/1‑69 (MQ=255)
                    aCATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGcc                                                     <  1:2031925/69‑1 (MQ=255)
                    aCATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGcc                                                     <  1:210124/69‑1 (MQ=255)
                    aCATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCAAAACa                                                                >  1:2672573/1‑58 (MQ=255)
                    aCATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCAAAACAGTATCAGGGCCaa                                                   >  1:343497/1‑71 (MQ=255)
                      aTTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATg                                                 <  1:2443806/71‑1 (MQ=255)
                        tACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCaaaa                                                                  >  1:2027490/1‑52 (MQ=255)
                        tACCAGCCAAAATGCGGCAATTACTGCGGAACGTAACTGGCTCATTTCaaaa                                                                  >  1:2001283/1‑52 (MQ=255)
                        tACCAGACAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCaaaa                                                                  >  1:7449/1‑52 (MQ=255)
                         aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATgg                                                >  1:871225/1‑69 (MQ=255)
                         aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGt                                               >  1:3927302/1‑70 (MQ=255)
                         aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGt                                               >  1:3494579/1‑70 (MQ=255)
                         aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGt                                               >  1:434957/1‑70 (MQ=255)
                             gCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTAt                                                            <  1:91507/53‑1 (MQ=255)
                               cAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACctgc                                        >  1:214949/1‑71 (MQ=255)
                               cAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACctgc                                        >  1:2035842/1‑71 (MQ=255)
                                 aaaTGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGt                                                              <  1:3104113/47‑1 (MQ=255)
                                 aaaTGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTa                                                             <  1:3322301/48‑1 (MQ=255)
                                 aaaTGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTa                                                             <  1:4074723/48‑1 (MQ=255)
                                  aaTGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATc                                                           >  1:268754/1‑49 (MQ=255)
                                  aaTGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGcc                                                     >  1:448824/1‑55 (MQ=255)
                                    tGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCaaa                                                                   <  1:3662414/39‑1 (MQ=255)
                                      cGGCCATTACTGCGGAACGTAACTGGCTCATTTCaaa                                                                   >  1:3387082/1‑37 (MQ=255)
                                      cGGCCATTACTGCGGAACGTAACTGGCTCATTTCaaa                                                                   >  1:1924406/1‑37 (MQ=255)
                                      cGGCCATTACTGCGGAACGTAACTGGCTCATTTCaaa                                                                   >  1:2119038/1‑37 (MQ=255)
                                      cGGCCATTACTGCGGAACGTAACTGGCTCATTTCaaa                                                                   >  1:846903/1‑37 (MQ=255)
                                      cGGCCATTACTGCGGAACGTAACTGGCTCATTTCaaa                                                                   >  1:1832955/1‑37 (MQ=255)
                                      cGGCCATTACTGCGGAACGTAACTGGCTCATTTCaaa                                                                   >  1:3936981/1‑37 (MQ=255)
                                           aTTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTc                                              <  1:2159874/53‑1 (MQ=255)
                                            ttACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGt                                 <  1:1366301/65‑1 (MQ=255)
                                            ttACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTg                                <  1:3333971/66‑1 (MQ=255)
                                            ttACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTg                                <  1:2583260/66‑1 (MQ=255)
                                            ttACTGCGGAACGTAACTGGCTCATTCCAAAACAGt                                                              <  1:117099/36‑1 (MQ=255)
                                            ttACTGCGGAACGTAACTGGCTCATTCCAAAACAGt                                                              <  1:2421122/36‑1 (MQ=255)
                                             tACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCAc                                            >  1:2455898/1‑53 (MQ=255)
                                             tACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCAc                                            >  1:1378113/1‑53 (MQ=255)
                                             tACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCAc                                            >  1:3695777/1‑53 (MQ=255)
                                             tACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCAc                                            >  1:61641/1‑53 (MQ=255)
                                              aCTGCGGAACGTAACTGGCTCGTTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCt                         >  1:76663/1‑71 (MQ=255)
                                              aCTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCt                         >  1:438803/1‑71 (MQ=255)
                                              aCTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCt                         >  1:699831/1‑71 (MQ=255)
                                              aCTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCt                         >  1:3635048/1‑71 (MQ=255)
                                              aCTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCt                         >  1:1571921/1‑71 (MQ=255)
                                              aCTGCGGAACGTAACTGGCTCATTCCAAAACAGTAt                                                            >  1:662085/1‑36 (MQ=255)
                                              aCTGCGGAACGTAACTGGCTCATTCCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTc                          >  1:3348728/1‑70 (MQ=255)
                                                tGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTg                        >  1:1034523/1‑70 (MQ=255)
                                                  cGGAACGTAACTGGCTCATTCCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGa                       <  1:692559/69‑1 (MQ=255)
                                                     aaCGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGAc                  <  1:2545154/71‑1 (MQ=255)
                                                     aaCGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGAc                  <  1:2519412/71‑1 (MQ=255)
                                                       cGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGaaa                     >  1:1229064/1‑66 (MQ=255)
                                                       cGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGaa                      >  1:735337/1‑65 (MQ=255)
                                                           aCTGGCTCATTCCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGc             >  1:775411/1‑70 (MQ=255)
                                                            cTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCa           <  1:1949717/71‑1 (MQ=255)
                                                              ggCTCATTCCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGt                           >  1:654642/1‑53 (MQ=255)
                                                                  cATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCTTGCGCGTCTGAAAGa                   >  1:766494/1‑57 (MQ=255)
                                                                  cATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGa                   >  1:3438211/1‑57 (MQ=255)
                                                                  cATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGa                   >  1:1457784/1‑57 (MQ=255)
                                                                   aTTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGGGTCTGAAAGACATCGCCaa          >  1:2591958/1‑65 (MQ=255)
                                                                   aTTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCaa          >  1:1449504/1‑65 (MQ=255)
                                                                    tttCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGAc                  <  1:972815/56‑1 (MQ=255)
                                                                    tttCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGAc                  <  1:1888064/56‑1 (MQ=255)
                                                                    tttCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGAc                  <  1:1518110/56‑1 (MQ=255)
                                                                    tttCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   <  1:2867058/71‑1 (MQ=255)
                                                                    tttCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   <  1:1710998/71‑1 (MQ=255)
                                                                    tttCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   <  1:899150/71‑1 (MQ=255)
                                                                      tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   >  1:378282/1‑69 (MQ=255)
                                                                      tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   >  1:306279/1‑69 (MQ=255)
                                                                      tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   >  1:773278/1‑69 (MQ=255)
                                                                      tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   >  1:284094/1‑69 (MQ=255)
                                                                      tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   >  1:1301992/1‑69 (MQ=255)
                                                                      tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTACa  <  1:626429/70‑1 (MQ=255)
                                                                      tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTACa  <  1:554136/70‑1 (MQ=255)
                                                                      |                                                                     
TCTCCGGACCGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTACA  >  W3110S.gb/3903479‑3903618

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