Predicted mutation
evidence seq id position mutation freq annotation gene description
RA W3110S.gb 2,198,825 G→A 100% A398T (GCG→ACG)  metG → methionyl‑tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*W3110S.gb2,198,8250GA82.9% 81.1 / 9.8 35A398T (GCG→ACG) metGmethionyl‑tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (6/0);  new base A (15/14);  total (21/14)
Fisher's exact test for biased strand distribution p-value = 6.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

GGAAGATTTCGTTCAGCGTGTGAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATGCGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAGCGAACTGGCTG  >  W3110S.gb/2198758‑2198876
                                                                   |                                                   
ggAAGATTTCGTTCAGCGTGTGAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATGCgg                                                  >  1:2620731/1‑71 (MQ=255)
        tCGTTCAGCGTGTGAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATGCGGGCTTTa                                            >  1:942769/1‑69 (MQ=255)
        tCGTTCAGCGTGTGAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATGCGGGCTTTAt                                           >  1:3517308/1‑70 (MQ=255)
         cGTTCAGCGTGTGAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATc                                          <  1:3312783/70‑1 (MQ=255)
         cGTTCAGCGTGTGAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATc                                          <  1:3179602/70‑1 (MQ=255)
         cGTTCAGCGTGTGAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATc                                          >  1:1837866/1‑70 (MQ=255)
         cGTTCAGCGTGTGAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATc                                          >  1:760264/1‑70 (MQ=255)
                    tgAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATGCGGGCTTTATCAACAAGCGtttt                              >  1:3385634/1‑71 (MQ=255)
                    tgAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATGCGGGCTTTATCAACAAGCGtttt                              >  1:2420327/1‑71 (MQ=255)
                    tgAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGtttt                              <  1:1626370/71‑1 (MQ=255)
                    tgAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGtttt                              <  1:276092/71‑1 (MQ=255)
                      aaTGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTg                             <  1:2350871/70‑1 (MQ=255)
                          ccGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGc                        >  1:62848/1‑71 (MQ=255)
                               atCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCt                   >  1:76200/1‑71 (MQ=255)
                               atCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCt                   >  1:353126/1‑71 (MQ=255)
                               atCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCt                   >  1:2359648/1‑71 (MQ=255)
                                tCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCt                   <  1:125252/70‑1 (MQ=255)
                                tCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCt                   <  1:1386421/70‑1 (MQ=255)
                                tCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTg                  >  1:820248/1‑71 (MQ=255)
                                  gTTAACAAAGTGGTTAACCTGGCCTCCCGTAATGCGGGCTTTATcaa                                        >  1:2830077/1‑47 (MQ=255)
                                  gTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGAc                           >  1:455684/1‑60 (MQ=255)
                                  gTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGAc                           >  1:3420794/1‑60 (MQ=255)
                                  gTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGAc                           >  1:299214/1‑60 (MQ=255)
                                  gTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGAc                           >  1:3275959/1‑60 (MQ=255)
                                  gTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTgg                 >  1:3155565/1‑70 (MQ=255)
                                  gTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTgg                 >  1:3563991/1‑70 (MQ=255)
                                  gTTAACAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTgg                 >  1:408977/1‑70 (MQ=255)
                                     aaCAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAg             >  1:1781693/1‑71 (MQ=255)
                                       cAAAGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAGCg           <  1:2762463/71‑1 (MQ=255)
                                       cAAAGTGGTTAACATGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAGCg           <  1:1586329/71‑1 (MQ=255)
                                        aaaGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAGCg           <  1:1259072/70‑1 (MQ=255)
                                        aaaGTGGTTAACCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAGCGa          <  1:3251857/71‑1 (MQ=255)
                                                 aaCCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAGCGAACTGGCTg  <  1:2703005/70‑1 (MQ=255)
                                                 aaCCTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAGCGAACTGGCTg  <  1:245935/70‑1 (MQ=255)
                                                    cTGGCCTCCCGTAATACGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAGCGa          <  1:1669871/59‑1 (MQ=255)
                                                                   |                                                   
GGAAGATTTCGTTCAGCGTGTGAATGCCGATATCGTTAACAAAGTGGTTAACCTGGCCTCCCGTAATGCGGGCTTTATCAACAAGCGTTTTGACGGCGTGCTGGCAAGCGAACTGGCTG  >  W3110S.gb/2198758‑2198876

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: