Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,903,549 T→C 10.9% 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,5490TC10.9% 140.7 / 10.0 55S104P (TCA→CCA) baxconserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/20);  new base C (4/2);  total (33/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 = 8.95e-01

GGACCGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTACA  >  W3110S.gb/3903484‑3903618
                                                                 |                                                                     
ggACCGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCaaaa                                                                  <  1:2044989/71‑1 (MQ=255)
   ccGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAg                                                               <  1:2811176/71‑1 (MQ=255)
   ccGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAg                                                               <  1:2213009/71‑1 (MQ=255)
      tAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTAt                                                            <  1:2653526/71‑1 (MQ=255)
         tGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAg                                                         <  1:1706830/71‑1 (MQ=255)
           ccTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAggg                                                       >  1:1799834/1‑71 (MQ=255)
           ccTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCAAAACAGTATCAggg                                                       >  1:1058506/1‑71 (MQ=255)
               aCATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAg                                                               <  1:1974441/59‑1 (MQ=255)
               aCATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCaa                                                   >  1:794773/1‑71 (MQ=255)
                cATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTa                                                             >  1:1242656/1‑60 (MQ=255)
                cATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCaa                                                   <  1:2715050/70‑1 (MQ=255)
                 aTTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTa                                                             >  1:1053927/1‑59 (MQ=255)
                 aTTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTa                                                             >  1:1081886/1‑59 (MQ=255)
                 aTTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCaa                                                                    <  1:671409/52‑1 (MQ=255)
                  ttACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATgg                                                <  1:240979/71‑1 (MQ=255)
                   tACCAGCCAAAATGGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGt                                               >  1:1716557/1‑70 (MQ=255)
                    aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAg                                                         >  1:1220815/1‑60 (MQ=255)
                    aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAg                                                         >  1:1706356/1‑60 (MQ=255)
                    aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAg                                                         >  1:1669841/1‑60 (MQ=255)
                    aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCaa                                                   >  1:979830/1‑66 (MQ=255)
                    aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGt                                               >  1:742955/1‑70 (MQ=255)
                    aCCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTCCAAAACAGTa                                                             >  1:1126488/1‑56 (MQ=255)
                     ccAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTc                                              <  1:745764/70‑1 (MQ=255)
                     ccAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCa                                             <  1:270191/71‑1 (MQ=255)
                        gCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGt                                                              <  1:2417679/51‑1 (MQ=255)
                         ccAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACctg                                         >  1:2272182/1‑71 (MQ=255)
                         ccAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACctg                                         >  1:1893790/1‑71 (MQ=255)
                          cAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTAt                                                            >  1:2428533/1‑51 (MQ=255)
                          cAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTAt                                                            >  1:1642966/1‑51 (MQ=255)
                            aaaTGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGt                                                              <  1:377291/47‑1 (MQ=255)
                                gCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCg                                  >  1:2699275/1‑71 (MQ=255)
                                gCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCg                                  >  1:2702459/1‑71 (MQ=255)
                                     cATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCaa                                                   >  1:2346183/1‑49 (MQ=255)
                                     cATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCaa                                                   >  1:24676/1‑49 (MQ=255)
                                      aTTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGcgcg                            <  1:152467/71‑1 (MQ=255)
                                       ttACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGcc                                                     <  1:2274634/45‑1 (MQ=255)
                                       ttACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTc                                              <  1:1290794/52‑1 (MQ=255)
                                       ttACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTc                                              <  1:1414854/52‑1 (MQ=255)
                                       ttACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTc                                              <  1:1475082/52‑1 (MQ=255)
                                        tACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATc                                                           <  1:63188/38‑1 (MQ=255)
                                        tACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTc                          >  1:190380/1‑71 (MQ=255)
                                         aCTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCAc                                            >  1:1150750/1‑52 (MQ=255)
                                         aCTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCt                         >  1:1053985/1‑71 (MQ=255)
                                         aCTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCt                         >  1:1863659/1‑71 (MQ=255)
                                         aCTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCt                         >  1:2696584/1‑71 (MQ=255)
                                         aCTGCGGAACGTAACTGGCTCATTTCAAAACAGTATAAGGGCCAATGGTCa                                             >  1:605597/1‑51 (MQ=255)
                                         aCTGCGGAACGTAACTGGCTCATTCCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCt                         >  1:2671618/1‑71 (MQ=255)
                                               gAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGa                   >  1:1768655/1‑71 (MQ=255)
                                                aaCGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGAc                  <  1:1861576/71‑1 (MQ=255)
                                                    tAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACctg                                         <  1:2565324/44‑1 (MQ=255)
                                                    tAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACctg                                         <  1:565954/44‑1 (MQ=255)
                                                     aaCTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGc             >  1:1179878/1‑71 (MQ=255)
                                                     aaCTGGCTCATTCCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGc             >  1:1149509/1‑71 (MQ=255)
                                                      aCTGGTTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGc             >  1:1537715/1‑70 (MQ=255)
                                                         ggCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCaaa         >  1:2346783/1‑71 (MQ=255)
                                                               tttCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAg                                    <  1:631872/38‑1 (MQ=255)
                                                               tttCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAg                                    <  1:97789/38‑1 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGaa                      >  1:882350/1‑50 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACAt                >  1:806492/1‑56 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACAt                >  1:2720731/1‑56 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCACACGCTAc   >  1:1604416/1‑69 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAAc        >  1:809514/1‑64 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   >  1:1236939/1‑69 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   >  1:1870693/1‑69 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTAc   >  1:1589430/1‑69 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTACa  <  1:1223818/70‑1 (MQ=255)
                                                                 tCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTACa  <  1:945269/70‑1 (MQ=255)
                                                                 |                                                                     
GGACCGTAATGCCTTACATTACCAGCCAAAATGCGGCCATTACTGCGGAACGTAACTGGCTCATTTCAAAACAGTATCAGGGCCAATGGTCACCTGCTGAGCGTGCGCGTCTGAAAGACATCGCCAAACGCTACA  >  W3110S.gb/3903484‑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: