Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 139,618 A→G E118G (GAA→GGA)  yadG → ABC transporter ATP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP009273139,6180AG92.3% 67.3 / ‑0.6 26E118G (GAA→GGA) yadGABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/0);  new base G (14/10);  total (16/10)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

TCAACCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGAAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTAA  >  NZ_CP009273/139531‑139701
                                                                                       |                                                                                   
tCAACCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGgaa                                                                                   >  2:282986/1‑90 (MQ=255)
 cAACCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGgaaa                                                                                  >  1:139131/1‑90 (MQ=255)
  aaCCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGAACTACGGCGGGGAGCGCAAAGAAGCGGACATCCGCAGCGGAAAg                                                                                 >  2:505020/1‑90 (MQ=255)
       gTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCt                                                                            >  1:283992/1‑90 (MQ=255)
                cGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACt                                                                   <  2:199768/90‑1 (MQ=255)
                 gTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTc                                                                  >  1:559350/1‑90 (MQ=255)
                 gTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTc                                                                  >  1:389464/1‑90 (MQ=255)
                  tGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCg                                                                 <  1:53155/90‑1 (MQ=255)
                                       aggcaggGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAg                                                                                 <  1:135009/53‑1 (MQ=255)
                                       aggcaggGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAg                                                                                 >  2:135009/1‑53 (MQ=255)
                                                actacGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTgcgc                                   <  1:59115/90‑1 (MQ=255)
                                                 ctacGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGcgcg                                  <  1:264714/90‑1 (MQ=255)
                                                   acGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTa                                <  1:538745/90‑1 (MQ=255)
                                                      gCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGt                             >  1:35414/1‑90 (MQ=255)
                                                      gCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGt                             >  1:242210/1‑90 (MQ=255)
                                                      gCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGt                             >  1:56083/1‑90 (MQ=255)
                                                        gTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTa                           >  1:104330/1‑90 (MQ=255)
                                                        gTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTa                           >  1:496721/1‑90 (MQ=255)
                                                           gAGCGCAAAGAAGCGTACATCCGCAGCGAAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCt                        >  2:348449/1‑90 (MQ=255)
                                                             gcgcAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTgg                      >  1:362865/1‑90 (MQ=255)
                                                                       gCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAg            >  2:280166/1‑90 (MQ=255)
                                                                         gTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAgcg          >  1:104358/1‑90 (MQ=255)
                                                                          tACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAgcgc         <  2:303286/90‑1 (MQ=255)
                                                                          tACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAgcgc         <  2:486627/90‑1 (MQ=255)
                                                                               ccGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGttt    <  2:559350/90‑1 (MQ=255)
                                                                                cGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTa   <  2:104330/90‑1 (MQ=255)
                                                                                 gCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTaa  >  1:251751/1‑90 (MQ=255)
                                                                                 gCAGCGAAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGATATCTGGCGGGATGAAGCGCCGTTTaa  >  2:252932/1‑90 (MQ=255)
                                                                                       |                                                                                   
TCAACCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGAAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTAA  >  NZ_CP009273/139531‑139701

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