Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,971,401 A→T 56.2% V373E (GTG→GAG)  ygeD ← predicted inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,971,4010AT56.2% 19.9 / 45.6 41V373E (GTG→GAG) ygeDpredicted inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (13/5);  new base T (14/9);  total (27/14)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGATCCACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTG  >  W3110S.gb/2971341‑2971463
                                                            |                                                              
cAGATCCACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGAt                                                      <  1:457446/71‑1 (MQ=255)
cAGATCCACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGAt                                                      <  1:454323/71‑1 (MQ=255)
   aTCCACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCgg                                                         <  1:643596/65‑1 (MQ=255)
    tCCACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGa                                                       >  1:1997118/1‑66 (MQ=255)
     ccACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTa                                                 >  1:47836/1‑71 (MQ=255)
     ccACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGa                                                       >  1:760530/1‑65 (MQ=255)
     ccACAGCGCCGTTATTGCCAGCGCAAACAGCGAACCGAAGCCAATGCCAATGGGCACGACCGGGa                                                       >  1:1465853/1‑65 (MQ=255)
          gcgcCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCAtt                                            <  1:737915/71‑1 (MQ=255)
                     gCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCATTACCGCCa                                     >  1:791500/1‑67 (MQ=255)
                     gCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCATTACCGCCa                                     >  1:705005/1‑67 (MQ=255)
                      ccAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAg                                >  1:948435/1‑71 (MQ=255)
                      ccAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAg                                >  1:1538532/1‑71 (MQ=255)
                        aGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTa                              <  1:1091755/71‑1 (MQ=255)
                        aGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTa                              <  1:2267498/71‑1 (MQ=255)
                                 aGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCAt                                             >  1:805552/1‑47 (MQ=255)
                                 aGCGCACCGAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCAt                                             >  1:804316/1‑47 (MQ=255)
                                    gcACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCAt                  >  1:1054332/1‑71 (MQ=255)
                                         gAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCa                      >  1:39134/1‑62 (MQ=255)
                                         gAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCa                      >  1:335640/1‑62 (MQ=255)
                                         gAAGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATCCCCa                      >  1:1670706/1‑62 (MQ=255)
                                          aaGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATcaacaa            <  1:2570126/71‑1 (MQ=255)
                                          aaGCCAATGCCAATGGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATcaacaa            <  1:2593413/71‑1 (MQ=255)
                                                atgccaatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGt                               <  1:217414/46‑1 (MQ=255)
                                                  gccaatgGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGgcgc    >  1:562885/1‑71 (MQ=255)
                                                  gccaatgGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGgcgc    >  1:559748/1‑71 (MQ=255)
                                                  gccaatgGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGgcgc    >  1:2208355/1‑71 (MQ=255)
                                                  gccaatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTa                              <  1:1355415/45‑1 (MQ=255)
                                                  gccaatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGgcgc    >  1:1727839/1‑71 (MQ=255)
                                                   ccaatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAacaaca          >  1:2484700/1‑64 (MQ=255)
                                                   ccaatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCt   >  1:1387058/1‑71 (MQ=255)
                                                   ccaatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCt   <  1:1224618/71‑1 (MQ=255)
                                                   ccaatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCt   >  1:1957077/1‑71 (MQ=255)
                                                   ccaatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCt   >  1:2551383/1‑71 (MQ=255)
                                                   ccaatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCt   >  1:2175473/1‑71 (MQ=255)
                                                   ccaatgGGCACGACCGGGATACCTATCATTACCGCCAGGGAGTAAATGCCCAGCATCAACAACATGGCGCt   >  1:2636925/1‑71 (MQ=255)
                                                    caatgGGCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTg  <  1:239194/71‑1 (MQ=255)
                                                    caatgGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTa                              <  1:1792140/43‑1 (MQ=255)
                                                        gggCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTaaa                            >  1:2109283/1‑41 (MQ=255)
                                                        gggCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTaaa                            >  1:296906/1‑41 (MQ=255)
                                                        gggCTCGACCGGGATGCCTATCATTACCGCCAGCGAGTaaa                            >  1:2112490/1‑41 (MQ=255)
                                                        gggCACGACCGGGATGCCTATCATTACCGCCAGCGAGTa                              <  1:926827/39‑1 (MQ=255)
                                                            aCGACCGGGATGCCTATCATTACCGCCAGCGAGTa                              <  1:74084/35‑1 (MQ=255)
                                                            |                                                              
CAGATCCACAGCGCCGTTATTGCCAGCGCAAACAGCGCACCGAAGCCAATGCCAATGGGCACGACCGGGATGCCTATCATTACCGCCAGCGAGTAAATGCCCAGCATCAACAACATGGCGCTG  >  W3110S.gb/2971341‑2971463

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: