Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,785,687 G→C G75G (GGC→GGG kbl ← glycine C‑acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,785,6870GC90.0% 26.4 / ‑6.7 10G75G (GGC→GGGkblglycine C‑acetyltransferase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base C (5/4);  total (5/5)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAA  >  NZ_CP009273/3785612‑3785771
                                                                           |                                                                                    
cAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTTCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTgaga                                                                        >  1:193639/1‑90 (MQ=255)
                                    gCTGTCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAgg                                    <  2:237570/90‑1 (MQ=255)
                                        tCCTGAGTGCCGCAAATAAACCGCCCCGAAGCCATGCCGAAACCGTGAGAACCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGgatga                                <  2:263673/90‑1 (MQ=255)
                                              gTGCCGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGcc                          <  2:159157/90‑1 (MQ=255)
                                                 ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            >  1:234896/1‑69 (MQ=255)
                                                 ccGCAAATAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAatc                                            <  2:234896/69‑1 (MQ=255)
                                                        tAAAACGCACCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGAt                >  1:284506/1‑90 (MQ=255)
                                                                aCCGAAGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGttgttg        >  2:120076/1‑90 (MQ=255)
                                                                    aaGCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGcac    >  1:305681/1‑90 (MQ=255)
                                                                      gCCATCCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACaa  <  1:291300/90‑1 (MQ=255)
                                                                           |                                                                                    
CAGGAAGGCCGCCAGTTTTTGTTCAAGCTCTTTATGGCTGTCCTGAGTGCCGCAAATAAAACGCACCGAAGCCATGCCGAAACCGTGAGAATCCATTCCCGCCTTTGCCGCCGCAATCAGATCAGGATGATTCGCCAGCCCGAGATAGTTGTTGGCACAA  >  NZ_CP009273/3785612‑3785771

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