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,6180AG83.9% 62.9 / 9.5 31E118G (GAA→GGA) yadGABC transporter ATP‑binding protein
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (16/10);  total (18/13)
Fisher's exact test for biased strand distribution p-value = 6.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.19e-01

TTCAACCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGAAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTAATG  >  NZ_CP009273/139530‑139703
                                                                                        |                                                                                     
ttCAACCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGaa                                                                                      <  2:675096/90‑1 (MQ=255)
   aaCCCGTTTGAAACCGTGCAGCAAATTGTTGTGAATCCGGCAGGGTACCACGGCGCGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAg                                                                                   <  2:365990/90‑1 (MQ=255)
   aaCCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAg                                                                                   <  2:434375/90‑1 (MQ=255)
   aaCCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAg                                                                                   <  2:149957/90‑1 (MQ=255)
                  gTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTc                                                                    >  2:712306/1‑90 (MQ=255)
                  gTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTc                                                                    >  1:794523/1‑90 (MQ=255)
                  gTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGAAAAGTATCTTAAACAACTc                                                                    >  1:102785/1‑90 (MQ=255)
                    gcagcaAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGAAAAGTATCTTAAACAACTCGa                                                                  <  1:272739/90‑1 (MQ=255)
                       gcaAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCt                                                               <  2:279357/90‑1 (MQ=255)
                          aaTTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATg                                                            <  2:794523/90‑1 (MQ=255)
                                 gtgAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGaaaa                                                     >  2:70147/1‑90 (MQ=255)
                                      tCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGAAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCaa                                                <  2:361670/90‑1 (MQ=255)
                                        aggcaggGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCaacg                                              >  2:158965/1‑90 (MQ=255)
                                             gggTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGt                                         <  1:14666/90‑1 (MQ=255)
                                                  ctacGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGcgcg                                    <  1:139487/90‑1 (MQ=255)
                                                       gCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGt                               >  2:758731/1‑90 (MQ=255)
                                                       gCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGt                               >  2:497639/1‑90 (MQ=255)
                                                       gCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGt                               >  1:435092/1‑90 (MQ=255)
                                                       gCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGAAAAg                                                                                   <  2:517433/38‑1 (MQ=255)
                                                       gCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGAAAAg                                                                                   >  1:517433/1‑38 (MQ=255)
                                                         gTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTa                             >  2:38875/1‑90 (MQ=255)
                                                         gTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTa                             >  1:427588/1‑90 (MQ=255)
                                                         gTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTa                             >  2:603490/1‑90 (MQ=255)
                                                         gTGGAGCGCAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTa                             >  2:143852/1‑90 (MQ=255)
                                                                gcAAAGAAGCGTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCg                      >  2:231242/1‑90 (MQ=255)
                                                                          gTACATCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAgcg            <  1:65806/90‑1 (MQ=255)
                                                                               tCCGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGtt       >  2:256249/1‑90 (MQ=255)
                                                                                ccGCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGttt      <  1:266277/90‑1 (MQ=255)
                                                                                  gCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTaa    >  2:135680/1‑90 (MQ=255)
                                                                                  gCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTaa    >  1:632800/1‑90 (MQ=255)
                                                                                  gCAGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTaa    >  1:840104/1‑90 (MQ=255)
                                                                                    aGCGGAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTAATg  <  1:542432/90‑1 (MQ=255)
                                                                                        |                                                                                     
TTCAACCCGTTTGAAACCGTGCAGCAAATTGTGGTGAATCAGGCAGGGTACTACGGCGTGGAGCGCAAAGAAGCGTACATCCGCAGCGAAAAGTATCTTAAACAACTCGATCTATGGGGAAAACGCAACGAACGTGCGCGTATGTTATCTGGCGGGATGAAGCGCCGTTTAATG  >  NZ_CP009273/139530‑139703

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