Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,888,711 A→G V143A (GTC→GCC)  ygcQ ← electron transfer flavoprotein subunit alpha/FixB family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,888,7110AG93.3% 38.6 / ‑3.2 15V143A (GTC→GCC) ygcQelectron transfer flavoprotein subunit alpha/FixB family protein
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (10/4);  total (11/4)
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.63e-01

CTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGAAAGGCATAGC  >  NZ_CP009273/2888638‑2888791
                                                                         |                                                                                
cTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGttgttg                                                                  <  1:156927/90‑1 (MQ=255)
       aaTGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCa                                                           >  2:367729/1‑90 (MQ=255)
            gTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCATTCCAGAggg                                                      >  1:65432/1‑90 (MQ=255)
               gcgGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCAt                                                               >  1:177040/1‑78 (MQ=255)
               gcgGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCAt                                                               <  2:177040/78‑1 (MQ=255)
                                             aaCCAGCCAGTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGAcg                     >  1:339257/1‑90 (MQ=255)
                                                cagccagTCCGGAAGCGCCCCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACgcgc                  <  2:277401/90‑1 (MQ=255)
                                                 agccagTCCGGAAGCGCCCCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGcc                 <  2:195420/90‑1 (MQ=255)
                                                  gccagTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCa                >  1:109210/1‑90 (MQ=255)
                                                  gccagTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCa                >  1:127403/1‑90 (MQ=255)
                                                  gccagTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCa                >  1:219485/1‑90 (MQ=255)
                                                  gccagTCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCa                >  2:10578/1‑90 (MQ=255)
                                                       tCCGGAAGCGCACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGaa           >  2:278780/1‑90 (MQ=255)
                                                                gcACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGAAAGGCATAGc  >  2:186709/1‑90 (MQ=255)
                                                                gcACCCGGGGCAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGAAAGGCATAGc  >  2:214338/1‑90 (MQ=255)
                                                                         |                                                                                
CTCAGCTAATGGGTCGCGGGTAACATTTTTCAGATCCTCTGTGCTAACCAGCCAGTCCGGAAGCGCACCCGGGACAATGTTCAGTTGTTGCATTCCAGAGGGCAACGTGGCGTTTTTCGCTGCTCCCGCCTGACGCGCCAGGGAAAGGCATAGC  >  NZ_CP009273/2888638‑2888791

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