Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,094,108 A→G E256G (GAA→GGA)  ghrA → glyoxylate/hydroxypyruvate reductase GhrA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,094,1080AG92.3% 33.9 / ‑3.2 13E256G (GAA→GGA) ghrAglyoxylate/hydroxypyruvate reductase GhrA
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (6/6);  total (7/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 = 2.77e-01

GTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAA  >  NZ_CP009273/1094026‑1094193
                                                                                  |                                                                                     
gtTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTa                                                                                <  1:110744/90‑1 (MQ=255)
        tgtgGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGa                                                                        >  2:240033/1‑90 (MQ=255)
        tgtgGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGa                                                                        >  2:51601/1‑90 (MQ=255)
          tgGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGaaa                                                                      >  2:244874/1‑90 (MQ=255)
             aaGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTc                                                                   <  1:332102/90‑1 (MQ=255)
                  gACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGctc                                                              <  1:55855/90‑1 (MQ=255)
                                  tGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCAcgcg                                              <  1:248368/90‑1 (MQ=255)
                                       aGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACg                                         <  2:292720/90‑1 (MQ=255)
                                             aaaGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAaca                                   >  2:387096/1‑90 (MQ=255)
                                             aaaGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAaca                                   >  2:84346/1‑90 (MQ=255)
                                                             tGTTGGATGTTTTTAATCGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGa                   <  2:77788/90‑1 (MQ=255)
                                                                tGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTa                >  1:144439/1‑90 (MQ=255)
                                                                              cGTGGACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGaa  >  2:342050/1‑90 (MQ=255)
                                                                                  |                                                                                     
GTTCATGTTGTGGAAGATGACCTGCTCGCGGCGCTGGATAGCGGCAAAGTTAAAGGCGCAATGTTGGATGTTTTTAATCGTGAACCCTTACCGCCTGAAAGTCCGCTCTGGCAACATCCACGCGTGACGATAACACCACATGTCGCCGCGATTACCCGTCCCGCTGAA  >  NZ_CP009273/1094026‑1094193

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