Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,360,083 T→C 27.7% F344S (TTC→TCC)  gltD → glutamate synthase, 4Fe‑4S protein, small subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,360,0830TC27.7% 69.9 / 32.9 47F344S (TTC→TCC) gltDglutamate synthase, 4Fe‑4S protein, small subunit
Reads supporting (aligned to +/- strand):  ref base T (20/14);  new base C (6/7);  total (26/21)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

ATGAAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAAT  >  W3110S.gb/3360014‑3360143
                                                                     |                                                            
atgaAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTc                                                             <  1:686703/71‑1 (MQ=255)
 tgaAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTc                                                             <  1:453384/70‑1 (MQ=255)
     gagaACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGt                                                        <  1:835677/71‑1 (MQ=255)
       gaACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTc                                                       <  1:3852132/70‑1 (MQ=255)
                gggTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTg                                              >  1:3098244/1‑70 (MQ=255)
                 ggTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTggg                                            <  1:1252437/71‑1 (MQ=255)
                   ttCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCAACGTc                                                       <  1:3002026/58‑1 (MQ=255)
                      ccgccgCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTc                                                             <  1:882197/49‑1 (MQ=255)
                      ccgccgCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTg                                       <  1:1669437/71‑1 (MQ=255)
                         ccgcGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCATCCg                                                 >  1:339044/1‑58 (MQ=255)
                          cgcgAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGt                                   >  1:2610306/1‑71 (MQ=255)
                          cgcgAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGt                                   <  1:2415270/71‑1 (MQ=255)
                            cgAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGa                                 <  1:2836503/71‑1 (MQ=255)
                                gTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCa                                                            <  1:913156/40‑1 (MQ=255)
                                 tGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGt                            <  1:1424354/71‑1 (MQ=255)
                                 tGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGt                            <  1:1753856/71‑1 (MQ=255)
                                   aaaaaCTCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGCTATTGAAGTGAACGGTaa                          >  1:6294/1‑71 (MQ=255)
                                   aaaaaCGCGCGGGAAGAAGGCGTATAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTaa                          >  1:3821765/1‑71 (MQ=255)
                                   aaaaaCGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCg                                                 >  1:4506880/1‑48 (MQ=255)
                                   aaaaaCGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTaa                          >  1:4402941/1‑71 (MQ=255)
                                   aaaaaCGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTaa                          >  1:2323809/1‑71 (MQ=255)
                                     aaaCGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTg                                  <  1:959737/61‑1 (MQ=255)
                                      aaCGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTg                                  <  1:2475609/60‑1 (MQ=255)
                                        cgcgcgGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGc                      >  1:1079185/1‑70 (MQ=255)
                                        cgcgcgGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCa                     <  1:4444843/71‑1 (MQ=255)
                                             gggAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTc                <  1:619876/71‑1 (MQ=255)
                                                 agaagGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCa                     <  1:3134597/62‑1 (MQ=255)
                                                 agaagGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgcg            >  1:3178033/1‑71 (MQ=255)
                                                  gaagGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgcgg           >  1:4476900/1‑71 (MQ=255)
                                                   aagGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgcggc          >  1:2325415/1‑71 (MQ=255)
                                                     ggCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGa                                      <  1:2660904/41‑1 (MQ=255)
                                                     ggCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgcggcg         >  1:3448540/1‑70 (MQ=255)
                                                     ggCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgcggcg         >  1:192274/1‑70 (MQ=255)
                                                     ggCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgcggcg         >  1:3591955/1‑70 (MQ=255)
                                                       cGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGaa                                >  1:2882987/1‑45 (MQ=255)
                                                       cGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTa       <  1:554496/70‑1 (MQ=255)
                                                         tAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGa                                      <  1:565514/37‑1 (MQ=255)
                                                          agagTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTaaaaa   >  1:802219/1‑71 (MQ=255)
                                                           gagTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgc             >  1:2449265/1‑60 (MQ=255)
                                                           gagTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTaaaaa   >  1:1954456/1‑70 (MQ=255)
                                                           gagTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTaaaa    >  1:2765519/1‑69 (MQ=255)
                                                           gagTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAAt  >  1:4477996/1‑71 (MQ=255)
                                                           gagTTCAAATCCAACGTCCAGCCGCTGGGTATTGaa                                     >  1:2993273/1‑36 (MQ=255)
                                                                cAAATTCCACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgcgg           >  1:3622880/1‑57 (MQ=255)
                                                                cAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCa               >  1:646268/1‑53 (MQ=255)
                                                                 aaaTTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTaa                          >  1:2337992/1‑41 (MQ=255)
                                                                 aaaTTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTaa                          >  1:1031198/1‑41 (MQ=255)
                                                                     |                                                            
ATGAAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAAT  >  W3110S.gb/3360014‑3360143

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: