Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 3,358,271 A→G 100% 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,2710AG92.9% 84.9 / ‑1.2 28K1262R (AAA→AGA) gltBglutamate synthase, large subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (16/10);  total (17/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATCGCC  >  W3110S.gb/3358203‑3358338
                                                                    |                                                                   
tGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCagaa                                                                   >  1:1880131/1‑71 (MQ=255)
  aTAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAAcc                                                                 <  1:2130303/71‑1 (MQ=255)
  aTAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAAcc                                                                 <  1:3040496/71‑1 (MQ=255)
          ctgctgAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCTATGAGCGCCAGAGCAGAACCTTCTGGt                                                          >  1:1959601/1‑70 (MQ=255)
          ctgctgAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGt                                                          >  1:1429661/1‑70 (MQ=255)
            gctgAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTc                                                        >  1:3503826/1‑70 (MQ=255)
             ctgAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGa                                                      <  1:642990/71‑1 (MQ=255)
                    cgcAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGc                                               <  1:985922/71‑1 (MQ=255)
                    cgcAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGc                                               <  1:3132966/71‑1 (MQ=255)
                           gctgcAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACacc                                         >  1:640242/1‑70 (MQ=255)
                                aaCAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGa                                                      <  1:3568108/52‑1 (MQ=255)
                                 aCAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGtt                                  >  1:1295861/1‑71 (MQ=255)
                                 aCAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGtt                                  >  1:190665/1‑71 (MQ=255)
                                   aGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCt                                >  1:2797763/1‑71 (MQ=255)
                                   aGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCt                                >  1:1836821/1‑71 (MQ=255)
                                          aCCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGgcg                         >  1:3451402/1‑71 (MQ=255)
                                          aCCGTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACAACGACCGTTCTGTCGgc                          >  1:2295933/1‑70 (MQ=255)
                                             gTTTGTCGATGAGCGCCAGAGCAGAACCTTCTGGTTc                                                        <  1:2604626/37‑1 (MQ=255)
                                                    gATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCttt                <  1:2532710/70‑1 (MQ=255)
                                                    gATGAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCttt                >  1:36368/1‑70 (MQ=255)
                                                       gAGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAgg            >  1:2279616/1‑71 (MQ=255)
                                                        aGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGc           >  1:1127968/1‑71 (MQ=255)
                                                        aGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGc           >  1:2222387/1‑71 (MQ=255)
                                                        aGCGCCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGc           >  1:3321371/1‑71 (MQ=255)
                                                         gcgcCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGc           <  1:285310/70‑1 (MQ=255)
                                                           gcCAGAGCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCtat        <  1:3277764/71‑1 (MQ=255)
                                                              agagCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATc     >  1:2759746/1‑71 (MQ=255)
                                                               gagCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCtat        <  1:2246924/67‑1 (MQ=255)
                                                                 gCAGAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATCGcc  >  1:3184897/1‑71 (MQ=255)
                                                                    |                                                                   
TGATAACGGCCTGCTGAACGCGCAGTTGCTGCAACAGGCGAAACCGTTTGTCGATGAGCGCCAGAGCAAAACCTTCTGGTTCGATATTCGCAACACCGACCGTTCTGTCGGCGCGTCGCTTTCAGGCTATATCGCC  >  W3110S.gb/3358203‑3358338

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: