Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,717,969 T→C E237G (GAA→GGA)  glyQ ← glycine‑‑tRNA ligase subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,717,9690TC86.8% 86.3 / 8.2 38E237G (GAA→GGA) glyQglycine‑‑tRNA ligase subunit alpha
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base C (17/16);  total (20/18)
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.22e-01

GATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT  >  NZ_CP009273/3717887‑3718051
                                                                                  |                                                                                  
gaTGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAgcgc                                                                             >  2:301853/1‑90 (MQ=255)
  tGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca                                                                           >  1:199058/1‑90 (MQ=255)
  tGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCca                                                                           >  1:653015/1‑90 (MQ=255)
       ttACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCagcagc                                                                      >  1:78719/1‑90 (MQ=255)
             gcATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCtgctgc                                                                >  2:781733/1‑90 (MQ=255)
                  cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATGCGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc                                                           <  2:718390/90‑1 (MQ=255)
                  cagcagGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc                                                           <  2:810707/90‑1 (MQ=255)
                        gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc                                                           >  1:551131/1‑84 (MQ=255)
                        gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc                                                           <  2:551131/84‑1 (MQ=255)
                        gTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTtctc                                                     >  1:801587/1‑90 (MQ=255)
                            aaGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTAc                                                 <  2:653015/90‑1 (MQ=255)
                                 gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc                                            >  2:631718/1‑90 (MQ=255)
                                 gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc                                            >  2:67112/1‑90 (MQ=255)
                                 gtgGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTc                                            >  1:689769/1‑90 (MQ=255)
                                   gggCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGa                                          >  1:387456/1‑90 (MQ=255)
                                      cggcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc                                       >  1:523339/1‑90 (MQ=255)
                                      cggcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc                                       <  2:517303/90‑1 (MQ=255)
                                        gcTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAg                                     >  1:163107/1‑90 (MQ=255)
                                           ttCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTg                                  >  2:327349/1‑90 (MQ=255)
                                                aaTACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAg                             <  2:501346/90‑1 (MQ=255)
                                                  tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa                           <  2:288569/90‑1 (MQ=255)
                                                  tACGCTCGTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGa                           <  2:10334/90‑1 (MQ=255)
                                                         gTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc                                       >  1:631754/1‑71 (MQ=255)
                                                         gTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGc                                       <  2:631754/71‑1 (MQ=255)
                                                         gTAGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCcac                    >  2:712828/1‑90 (MQ=255)
                                                           aGGCTGGCAGCGGCAGCGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACAt                  <  2:689769/90‑1 (MQ=255)
                                                                ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc                                                           >  2:340245/1‑44 (MQ=255)
                                                                ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc                                                           <  2:445254/44‑1 (MQ=255)
                                                                ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc                                                           <  1:340245/44‑1 (MQ=255)
                                                                ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGcttc                                                           >  1:445254/1‑44 (MQ=255)
                                                                ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACa                              >  2:825222/1‑73 (MQ=255)
                                                                ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACa                              <  1:825222/73‑1 (MQ=255)
                                                                ggcagcggcagcGGATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCcgcg             >  1:121735/1‑90 (MQ=255)
                                                                       gcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGa      <  1:298582/90‑1 (MQ=255)
                                                                       gcagcggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGa      <  2:265751/90‑1 (MQ=255)
                                                                           cggATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt  <  2:523339/90‑1 (MQ=255)
                                                                           cggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt  <  1:508479/90‑1 (MQ=255)
                                                                           cggATTTCCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGtt  <  2:45414/90‑1 (MQ=255)
                                                                                  |                                                                                  
GATGGCTTTACGCGCATCCAGCAGGTTGAAGCTGTGGGCGGCTTTCAGAATACGCTCGTAGGCTGGCAGCGGCAGCGGATTTTCCAGCGCCAGCAGCTGCTGCGCTTCTTTCTCGTACTGCTCGAAGCAGGTGAACAGGAAGTCCACATCCGCGTATTCGAAGTT  >  NZ_CP009273/3717887‑3718051

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