Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,358,271 A→G 72.7% K1262R (AAA→AGA)  gltB → glutamate synthase, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb3,358,2710AG72.7% 53.4 / 20.3 33K1262R (AAA→AGA) gltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base A (3/6);  new base G (11/13);  total (14/19)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATCG  >  W3110S.gb/3358207‑3358336
                                                                |                                                                 
aaCGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAAcc                                                               <  1:2279693/69‑1 (MQ=255)
aaCGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAAcc                                                               <  1:914956/69‑1 (MQ=255)
           gAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGa                                                    >  1:1371091/1‑69 (MQ=255)
           gAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGa                                                    >  1:3696637/1‑69 (MQ=255)
            aaCGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAAc                                                                <  1:966823/56‑1 (MQ=255)
            aaCGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAAc                                                                <  1:2403041/56‑1 (MQ=255)
            aaCGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGat                                                   >  1:173988/1‑69 (MQ=255)
                    gTTGCTGCAACAGTCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCaa                                           <  1:871065/69‑1 (MQ=255)
                    gTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCaa                                           <  1:142080/69‑1 (MQ=255)
                    gTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCaa                                           >  1:219701/1‑69 (MQ=255)
                     ttGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAac                                          <  1:230478/69‑1 (MQ=255)
                     ttGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAac                                          <  1:2633523/69‑1 (MQ=255)
                      tgctgcAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGat                                                   >  1:1190603/1‑59 (MQ=255)
                      tgctgcAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAaca                                         >  1:540856/1‑69 (MQ=255)
                      tgctgcAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAac                                          >  1:1927361/1‑68 (MQ=255)
                       gctgcAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGa                                                    <  1:3114027/57‑1 (MQ=255)
                       gctgcAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAaca                                         >  1:975680/1‑68 (MQ=255)
                        ctgcAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCg                                                     <  1:1242331/55‑1 (MQ=255)
                             aCAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACaccgaccg                                  <  1:1832968/69‑1 (MQ=255)
                             aCAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACaccgaccg                                  <  1:1722204/69‑1 (MQ=255)
                                 gCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTc                                                            >  1:1995671/1‑39 (MQ=255)
                                          tttGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGgcgc                      <  1:1033695/68‑1 (MQ=255)
                                            tGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGa                                                    >  1:3515211/1‑36 (MQ=255)
                                            tGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGa                                                    >  1:1228001/1‑36 (MQ=255)
                                                gATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGc                                             >  1:2602520/1‑39 (MQ=255)
                                                      cgcCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGc         <  1:1984946/69‑1 (MQ=255)
                                                      cgcCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGc         <  1:3236234/69‑1 (MQ=255)
                                                      cgcCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGc         <  1:528527/69‑1 (MQ=255)
                                                          agagCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCta       <  1:1795221/67‑1 (MQ=255)
                                                          agagCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCta       <  1:3722427/67‑1 (MQ=255)
                                                          agagCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCta       <  1:751155/67‑1 (MQ=255)
                                                              cAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATCg  >  1:2297799/1‑68 (MQ=255)
                                                              cAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATCg  >  1:166706/1‑68 (MQ=255)
                                                                |                                                                 
AACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATCG  >  W3110S.gb/3358207‑3358336

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