Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 4,182,652 C→T 29.4% V464M (GTG→ATG)  gspD ← general secretory pathway component, cryptic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb4,182,6520CT29.4% 48.8 / 23.1 34V464M (GTG→ATG) gspDgeneral secretory pathway component, cryptic
Reads supporting (aligned to +/- strand):  ref base C (10/14);  new base T (7/3);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 2.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTTGCGTTCGACGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAGGAT  >  W3110S.gb/4182585‑4182718
                                                                   |                                                                  
gTTTTGCGTTCGACGGTATTAAAGACGTTATCCCCTGACGTGGTCTGTGACCCGGATAGCACCGGcaca                                                                   <  1:983173/69‑1 (MQ=255)
gTTTTGCGTTCGACGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGcaca                                                                   <  1:1190201/69‑1 (MQ=255)
gTTTTGCGTTCGACGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGcaca                                                                   <  1:2247013/69‑1 (MQ=255)
gTTTTGCGTTCGACGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGcaca                                                                   <  1:1503949/69‑1 (MQ=255)
 ttttGCGTTCGACGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACAt                                                                  <  1:2323878/69‑1 (MQ=255)
 ttttGCGTTCGACGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACAt                                                                  <  1:3087512/69‑1 (MQ=255)
   ttGCGTTCGACGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATcc                                                                >  1:865821/1‑69 (MQ=255)
             cGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCt                                                               >  1:3294526/1‑60 (MQ=255)
             cGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTg                                                              >  1:3298990/1‑61 (MQ=255)
                 aTTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATTGCACCGGCACAt                                                                  <  1:3166131/53‑1 (MQ=255)
                 aTTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACAt                                                                  <  1:1914647/53‑1 (MQ=255)
                 aTTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACAt                                                                  <  1:228880/53‑1 (MQ=255)
                 aTTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGa                                                   >  1:1752814/1‑68 (MQ=255)
                 aTTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGa                                                   >  1:643393/1‑68 (MQ=255)
                  ttAAAGACGTTATCCCATGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTGCCCCACGTTGAAg                                                 >  1:1636102/1‑69 (MQ=255)
                            tATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTGCCCCACGTTGAAGGACGCGAGtt                                       >  1:2797678/1‑69 (MQ=255)
                            tATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTGCCCCACGTTGAAGGACGCGAGtt                                       >  1:3627717/1‑69 (MQ=255)
                               ccccTGAAGTGGTCTGTGACCCGGATAGCACCGGCATATCCTGCCCCACGTTGAAGGACGCGAGTttat                                    >  1:1104272/1‑69 (MQ=255)
                                    gAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAg                                         >  1:1038329/1‑59 (MQ=255)
                                    gAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAg                                         >  1:2085998/1‑59 (MQ=255)
                                          gTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTttatta                                  >  1:2679745/1‑60 (MQ=255)
                                          gTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGtt                         >  1:956112/1‑69 (MQ=255)
                                          gTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGtt                         >  1:1957700/1‑69 (MQ=255)
                                          gTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGtt                         >  1:2922756/1‑69 (MQ=255)
                                          gTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGtt                         >  1:1376678/1‑69 (MQ=255)
                                            cTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTAc                       <  1:1276471/69‑1 (MQ=255)
                                              gtgACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGttt                                      <  1:3231262/52‑1 (MQ=255)
                                              gtgACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGttt                                      <  1:1268698/52‑1 (MQ=255)
                                              gtgACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGttt                                      <  1:3289815/52‑1 (MQ=255)
                                              gtgACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGttt                                      <  1:2954850/52‑1 (MQ=255)
                                                   ccGGATAGCACCGGCATATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCt                 <  1:1548381/68‑1 (MQ=255)
                                                   ccGGATAGCACCGGCATATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCt                 <  1:3369710/68‑1 (MQ=255)
                                                   ccGGATAGCACCGGCATATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCt                 <  1:1863583/68‑1 (MQ=255)
                                                     ggATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTgg              <  1:1282243/69‑1 (MQ=255)
                                                     ggATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTgg              <  1:3688606/69‑1 (MQ=255)
                                                      gATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTg               <  1:3183884/67‑1 (MQ=255)
                                                             ccGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGc        >  1:2644713/1‑67 (MQ=255)
                                                              cGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAg     <  1:1179563/69‑1 (MQ=255)
                                                                  acaTCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAGGAt  >  1:2791363/1‑68 (MQ=255)
                                                                   |                                                                  
GTTTTGCGTTCGACGGTATTAAAGACGTTATCCCCTGAAGTGGTCTGTGACCCGGATAGCACCGGCACATCCTGCCCCACGTTGAAGGACGCGAGTTTATTATCCAGCGTTACGATGCTTGGGGTGGCGAGGAT  >  W3110S.gb/4182585‑4182718

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