Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 3,852,078 T→G T467P (ACG→CCG)  yidK ← solute:sodium symporter family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092733,852,0780TG82.5% 107.9 / 15.0 40T467P (ACG→CCG) yidKsolute:sodium symporter family transporter
Reads supporting (aligned to +/- strand):  ref base T (4/3);  new base G (12/21);  total (16/24)
Fisher's exact test for biased strand distribution p-value = 4.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.73e-01

CCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  NZ_CP009273/3851991‑3852162
                                                                                       |                                                                                    
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  1:792566/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:637160/90‑1 (MQ=255)
ccGCAAACGCATCTTTGAAGGTAAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGAc                                                                                    <  2:327021/90‑1 (MQ=255)
  gCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAg                                                                                  >  2:43710/1‑90 (MQ=255)
    aaaCGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGc                                                                                <  1:501171/90‑1 (MQ=255)
    aaaCGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGc                                                                                <  2:642883/90‑1 (MQ=255)
     aaCGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGcc                                                                               <  1:290047/90‑1 (MQ=255)
     aaCGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGcc                                                                               <  2:794254/90‑1 (MQ=255)
        gCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGc                                                                            <  1:659754/90‑1 (MQ=255)
        gCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGc                                                                            <  2:591240/90‑1 (MQ=255)
          aTCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGcac                                                                          >  2:585088/1‑90 (MQ=255)
          aTCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGcac                                                                          >  2:520882/1‑90 (MQ=255)
              ttGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATag                                                                      <  2:490780/90‑1 (MQ=255)
              ttGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATag                                                                      <  1:335080/90‑1 (MQ=255)
                   ggTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGaa                                                                 >  1:579784/1‑90 (MQ=255)
                     tGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAAt                                                               >  2:255114/1‑90 (MQ=255)
                      gAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATg                                                              >  2:651333/1‑90 (MQ=255)
                             gTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCg                                                       >  2:829817/1‑90 (MQ=255)
                                gcgcgcGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAAc                                                    <  1:57718/90‑1 (MQ=255)
                                gcgcgcGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAAc                                                    <  1:668418/90‑1 (MQ=255)
                                   cgcgGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTc                                                 <  1:829817/90‑1 (MQ=255)
                                    gcgGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCa                                                >  1:290564/1‑90 (MQ=255)
                                     cgGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCAc                                               >  2:630848/1‑90 (MQ=255)
                                     cgGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCAc                                               <  2:365930/90‑1 (MQ=255)
                                     cgGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCAc                                               <  1:203602/90‑1 (MQ=255)
                                         ttGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGa                                           <  2:95446/90‑1 (MQ=255)
                                          tGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGAt                                          <  2:391962/90‑1 (MQ=255)
                                          tGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGAt                                          <  2:85872/90‑1 (MQ=255)
                                                 aCCGATCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGAt                                   >  1:386992/1‑90 (MQ=255)
                                                     aTCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTg                               >  2:385680/1‑90 (MQ=255)
                                                      tCACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGa                              >  2:596089/1‑90 (MQ=255)
                                                       cACCAGCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGAt                             <  1:331863/90‑1 (MQ=255)
                                                            gCATCACGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGc                        >  2:622875/1‑90 (MQ=255)
                                                                 aCGACCACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGAtt                   <  1:402489/90‑1 (MQ=255)
                                                                     ccACGTTGATGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGc               <  1:637237/90‑1 (MQ=255)
                                                                       aCGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGt                                                         <  2:177229/46‑1 (MQ=255)
                                                                       aCGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGt                                                         >  1:177229/1‑46 (MQ=255)
                                                                       aCGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCa             >  2:299158/1‑90 (MQ=255)
                                                                              tGCAGAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCa      <  2:110931/90‑1 (MQ=255)
                                                                                  gAACGGACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGc  >  1:58759/1‑90 (MQ=255)
                                                                                       |                                                                                    
CCGCAAACGCATCTTTGAAGGTGAACGGCGTGGCGCGCGGTTTGATAAAACCGATCACCAGCATCACGACCACGTTGATGCAGAACGTACAGGCCAGCACATAGAGGAAATGGAAGTCGAACTTCACCAGATAGTTGATGGTGATGTAGCTGATTATGCCAATCCCCATCGC  >  NZ_CP009273/3851991‑3852162

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