Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,486,264 C→T 84.4% P752P (CCG→CCA valS ← valyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,486,2640CT84.4% 74.0 / 8.1 32P752P (CCG→CCAvalSvalyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base T (10/17);  total (12/20)
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.01e-01

GTCGGCAGTGATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATCGGATGCGCGAGGCGCAGCAGACCTTCCAGTACAGTCACCAGCGTATGGCGAGTACCGCGCAGT  >  W3110S.gb/4486198‑4486327
                                                                  |                                                               
gtCGGCAGTGATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTgg                                                               >  1:633329/1‑69 (MQ=255)
 tCGGCAGTGATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGa                                                              <  1:542315/69‑1 (MQ=255)
 tCGGCAGTGATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGa                                                              <  1:328534/69‑1 (MQ=255)
 tCGGCAGTGATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGa                                                              <  1:2409287/69‑1 (MQ=255)
 tCGGCAGTGATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGa                                                              <  1:1871639/69‑1 (MQ=255)
 tCGGCAGTGATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGa                                                              <  1:1922148/69‑1 (MQ=255)
      aGTGATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGATgc                                                           >  1:2386341/1‑67 (MQ=255)
         gATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGATGCGCGa                                                       >  1:953049/1‑68 (MQ=255)
               gCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGATGCGCGAGgcgc                                                  >  1:3468788/1‑67 (MQ=255)
                 aaaGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCgcagca                                              <  1:2126962/69‑1 (MQ=255)
                 aaaGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCgcagca                                              <  1:724558/69‑1 (MQ=255)
                   aGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCGCAGCAGa                                            <  1:1992052/69‑1 (MQ=255)
                   aGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCGCAGCAGa                                            <  1:2157270/69‑1 (MQ=255)
                    gTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGa                                                              <  1:1206638/50‑1 (MQ=255)
                    gTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCGCAGCAGAc                                           <  1:754377/69‑1 (MQ=255)
                    gTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATCGGATGCGCGAGGCGCAGCAGAc                                           >  1:1018608/1‑69 (MQ=255)
                      aCTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCGCAGCAGACCt                                         >  1:3522728/1‑69 (MQ=255)
                      aCTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATCGGATGcgcg                                                        <  1:1055248/54‑1 (MQ=255)
                       cTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATCGGATGcgcg                                                        <  1:612972/53‑1 (MQ=255)
                         ttCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATCGGATGCGCGAGGCGCAGCAGACCTTcc                                      >  1:2362695/1‑69 (MQ=255)
                              acGCTGCCAGATGGTTTCGGTGATGAACGGAATGATTGGa                                                              <  1:2096454/40‑1 (MQ=255)
                                         tGGTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCGCAGCAGACCTTCCAGTACAGTCACCAGCg                      <  1:3268021/69‑1 (MQ=255)
                                          ggTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCGCAGCAGACCTTCCAg                                    <  1:2397534/54‑1 (MQ=255)
                                          ggTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCGCAGCAGACCTTCCAg                                    >  1:1959303/1‑54 (MQ=255)
                                          ggTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCGCAGCAGACCTTCCAg                                    <  1:387069/54‑1 (MQ=255)
                                          ggTTTCGGTGATGAACGGAATGATTGGATGCGCGAGGCGCAGCAGACCTTCCAg                                    >  1:723087/1‑54 (MQ=255)
                                                  tgatgaACGGAATGATTGGATGCGCGAGGCGCAGCAGACCTTCCAGTACAGTCACCAGCGTATGGCGAg             <  1:1031368/69‑1 (MQ=255)
                                                  tgatgaACGGAATGATTGGATGCGCGAGGCGCAGCAGACCTTCCAGTACAGTCACCAGCGTATGGCGAg             <  1:80831/69‑1 (MQ=255)
                                                   gatgaACGGAATGATCGGATGAGCGAGGCGCAGCAGAACTTCCAGTACAGTCACCAGCGTATGGCg               <  1:2165662/66‑1 (MQ=255)
                                                         cGGAATGATTGGATGCGCGAGGCGCAGCAGACCTTCCAGTACAGTCACCAGCGTATGGCGAGTAcc         >  1:3711121/1‑66 (MQ=255)
                                                             aTGATTGGATGCGCGAGGCGCAGCAGACCTTCCAGTACAGTCACCAGCGTATGGCGAGTACCGCGCAGt  >  1:1785912/1‑69 (MQ=255)
                                                             aTGATTGGATGCGCGAGGCGCAGCAGACCTTCCAGTACAGTCACCAGCGTATGGCGAGTACCGCGCAGt  >  1:2431911/1‑69 (MQ=255)
                                                                  |                                                               
GTCGGCAGTGATACCGCAAAGTACTTTCACACGCTGCCAGATGGTTTCGGTGATGAACGGAATGATCGGATGCGCGAGGCGCAGCAGACCTTCCAGTACAGTCACCAGCGTATGGCGAGTACCGCGCAGT  >  W3110S.gb/4486198‑4486327

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: