Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,360,083 T→C 40.0% 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,0830TC40.0% 14.6 / 32.4 30F344S (TTC→TCC) gltDglutamate synthase, 4Fe‑4S protein, small subunit
Reads supporting (aligned to +/- strand):  ref base T (9/9);  new base C (9/3);  total (18/12)
Fisher's exact test for biased strand distribution p-value = 2.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGAAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAATGGTGC  >  W3110S.gb/3360014‑3360148
                                                                     |                                                                 
atgaAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTc                                                                  <  1:3545157/71‑1 (MQ=255)
  gaAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCaa                                                                >  1:3965760/1‑71 (MQ=255)
  gaAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCaa                                                                >  1:1706603/1‑71 (MQ=255)
  gaAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCaa                                                                >  1:2918227/1‑71 (MQ=255)
            tGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAAc                                                               >  1:1309806/1‑62 (MQ=255)
             gCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCaa                                                                >  1:2967732/1‑60 (MQ=255)
                       cgccgcGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGa                                           <  1:822984/71‑1 (MQ=255)
                           gcgAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTg                                       <  1:2077514/71‑1 (MQ=255)
                               aGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTg                                            <  1:2352106/62‑1 (MQ=255)
                                  gAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGaa                                          >  1:1889779/1‑61 (MQ=255)
                                  gAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTATTGaa                                          >  1:288282/1‑61 (MQ=255)
                                   aaaaaCGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTaa                               >  1:3482376/1‑71 (MQ=255)
                                        cgcgcgGGAAGAAGGCGTAGAGTTCAAATCCAACGTCCAGCCGCTGGGTAt                                              >  1:1736190/1‑51 (MQ=255)
                                                aagaagGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgc                  >  1:1390028/1‑71 (MQ=255)
                                                aagaagGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgc                  >  1:1584438/1‑71 (MQ=255)
                                                         tAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGa                                           <  1:278734/37‑1 (MQ=255)
                                                         tAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGa                                           <  1:1044773/37‑1 (MQ=255)
                                                         tAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCaaa                        <  1:2734526/56‑1 (MQ=255)
                                                         tAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCaaa                        <  1:2967049/56‑1 (MQ=255)
                                                         tAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCaaa                        <  1:381569/56‑1 (MQ=255)
                                                         tAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCa                          <  1:2123141/54‑1 (MQ=255)
                                                           gagTTCAAATTCAACGTCCAGCCGCTGGGTATTGaa                                          >  1:1436404/1‑36 (MQ=255)
                                                           gagTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAgcgg                >  1:2272121/1‑62 (MQ=255)
                                                           gagTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAAt       <  1:3718964/71‑1 (MQ=255)
                                                           gagTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAAt       >  1:3447861/1‑71 (MQ=255)
                                                           gagTTCAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAAt       <  1:1704859/71‑1 (MQ=255)
                                                              ttCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAATGGt    >  1:1583358/1‑71 (MQ=255)
                                                                cAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTaa                               >  1:759358/1‑42 (MQ=255)
                                                                cAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAATGGTGc  >  1:848089/1‑71 (MQ=255)
                                                                cAAATCCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAATGGTGc  >  1:3989245/1‑71 (MQ=255)
                                                                     |                                                                 
ATGAAGAGAACATGCCGGGTTCCCGCCGCGAAGTGAAAAACGCGCGGGAAGAAGGCGTAGAGTTCAAATTCAACGTCCAGCCGCTGGGTATTGAAGTGAACGGTAACGGCAAAGTCAGCGGCGTAAAAATGGTGC  >  W3110S.gb/3360014‑3360148

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: