Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 120,483 A→C K327T (AAA→ACA)  aceE → pyruvate dehydrogenase (acetyl‑transferring), homodimeric type

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273120,4830AC87.5% 40.1 / ‑0.1 16K327T (AAA→ACA) aceEpyruvate dehydrogenase (acetyl‑transferring), homodimeric type
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (9/5);  total (10/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

GGATGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGC  >  NZ_CP009273/120400‑120569
                                                                                   |                                                                                      
ggATGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGaa                                                                                  >  1:278339/1‑90 (MQ=255)
      aCTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGAAAGATgg                                                                            <  2:199092/90‑1 (MQ=255)
         gctgcGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGAAAGATGGTGc                                                                         >  1:230454/1‑90 (MQ=255)
         gctgcGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTAACATCGAAAGAGGGTGc                                                                         >  1:139082/1‑90 (MQ=255)
           tgcGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGAAAGATGGTGCGt                                                                       <  1:7127/90‑1 (MQ=255)
             cGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGAAAGATGGTGCGTAc                                                                     <  1:554091/90‑1 (MQ=255)
                      aCCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGAGACTACCAGACCTTCACATCGAAAGATGGTGCGTACGTTCGTGaa                                                            >  2:410097/1‑90 (MQ=255)
                            ggTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGAAAGATGGTGCGTACGTTCGTGAACActtc                                                      >  1:246481/1‑90 (MQ=255)
                            ggTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGAAAGATGGTGCGTACGTTCGTGAACActtc                                                      >  1:330107/1‑90 (MQ=255)
                              tAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGAAAGATGGTGCGTACGTTCGTGAACActtctt                                                    <  2:246481/90‑1 (MQ=255)
                                    gATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGAAAGATGGGGCGTACGTTCGTGAACACTTCTTCGGTaa                                              >  2:589311/1‑90 (MQ=255)
                                      tCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATgg                                                                            >  1:205265/1‑58 (MQ=255)
                                      tCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATgg                                                                            <  2:205265/58‑1 (MQ=255)
                                         aGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCACATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATc                                         >  1:276588/1‑90 (MQ=255)
                                                                 gCGACTACCAGACCTTCACATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAg                 >  1:121210/1‑90 (MQ=255)
                                                                                tCACATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGc  <  1:410097/90‑1 (MQ=255)
                                                                                   |                                                                                      
GGATGAACTGCTGCGTAAGGATACCAGCGGTAAACTGATCCAGCTGATGAACGAAACCGTTGACGGCGACTACCAGACCTTCAAATCGAAAGATGGTGCGTACGTTCGTGAACACTTCTTCGGTAAATATCCTGAAACCGCAGCACTGGTTGCAGACTGGACTGACGAGC  >  NZ_CP009273/120400‑120569

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