Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,814,267 A→G A368A (GCT→GCC alaS ← alanine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,814,2670AG87.5% 34.0 / ‑0.5 16A368A (GCT→GCCalaSalanine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (3/11);  total (4/12)
Fisher's exact test for biased strand distribution p-value = 4.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

GCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCAT  >  NZ_CP009273/2814187‑2814347
                                                                                |                                                                                
gCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGc                                                                         <  2:80278/90‑1 (MQ=255)
  aGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGTCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGctc                                                                       <  2:124421/90‑1 (MQ=255)
             ccAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTtcttcagtct                                                            >  1:239438/1‑90 (MQ=255)
                    aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                     <  1:165923/90‑1 (MQ=255)
                    aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                     <  2:33148/90‑1 (MQ=255)
                    aTCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCAc                                                     <  1:252165/90‑1 (MQ=255)
                        ccAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGc                                                 <  1:173922/90‑1 (MQ=255)
                          aGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTc                                               >  2:307933/1‑90 (MQ=255)
                                                         gCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCa                        <  1:284485/82‑1 (MQ=255)
                                                         gCCAGACCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCa                        >  2:284485/1‑82 (MQ=255)
                                                              aCCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGc           <  1:230529/90‑1 (MQ=255)
                                                              aCCGCGCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGc           <  1:55901/90‑1 (MQ=255)
                                                                cgcgCTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCag         <  1:311930/90‑1 (MQ=255)
                                                                  cgcTCCAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCagag       >  1:118936/1‑90 (MQ=255)
                                                                       cAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCAt  <  1:307933/90‑1 (MQ=255)
                                                                       cAGAGTACGGGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCAt  <  2:118936/90‑1 (MQ=255)
                                                                                |                                                                                
GCAGTTTCACCATCCAGCGTATCACCAGAAAGTTTTGCCAGCTCTTCATCCAGCAACGCCAGACCGCGCTCCAGAGTACGAGCAAACTGCTCTTCTTCAGTCTTCAGCACCTGCTCAACCTGCGCCTGCTGGCGTTTCAGGTCTTCACCCGCAGAGCCCAT  >  NZ_CP009273/2814187‑2814347

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