Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,345,414 T→G H147P (CAC→CCC)  glpT ← glycerol‑3‑phosphate transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,345,4140TG97.6% 140.4 / ‑4.0 42H147P (CAC→CCC) glpTglycerol‑3‑phosphate transporter
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (23/18);  total (23/19)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

AGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC  >  NZ_CP009273/2345330‑2345500
                                                                                    |                                                                                      
agcagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCAcc                                                                                   <  1:642463/90‑1 (MQ=255)
agcagcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCAcc                                                                                   <  1:214135/90‑1 (MQ=255)
   agcGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATa                                                                                <  1:932175/90‑1 (MQ=255)
      ggcggAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATa                                                                                <  2:508568/87‑1 (MQ=255)
      ggcggAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATa                                                                                >  1:508568/1‑87 (MQ=255)
      ggcggAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTa                                                                             >  2:598350/1‑90 (MQ=255)
                 accaccGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACAcgg                                                                  >  1:868607/1‑90 (MQ=255)
                 accaccGACGCTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACAcgg                                                                  >  1:938271/1‑90 (MQ=255)
                      cGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGcc                                                             <  1:779830/90‑1 (MQ=255)
                            tgtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      >  2:412523/1‑59 (MQ=255)
                            tgtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  1:412523/59‑1 (MQ=255)
                            tgtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCa                                                       >  1:777093/1‑90 (MQ=255)
                                 gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATcccc                                                  <  1:733308/90‑1 (MQ=255)
                                    caGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCCCCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTgg                                               <  1:101968/90‑1 (MQ=255)
                                     aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATcccc                                                  >  1:424660/1‑86 (MQ=255)
                                     aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATcccc                                                  <  2:424660/86‑1 (MQ=255)
                                     aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGa                                              >  2:376346/1‑90 (MQ=255)
                                              cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCAt                                                      <  1:248760/73‑1 (MQ=255)
                                              cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCAt                                                      >  2:248760/1‑73 (MQ=255)
                                              cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGc                                     >  1:385352/1‑90 (MQ=255)
                                                gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGc                                                                                      <  2:600321/39‑1 (MQ=255)
                                                gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCAcc                                                                                   <  1:237547/42‑1 (MQ=255)
                                                gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCAcc                                                                                   >  2:237547/1‑42 (MQ=255)
                                                gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:149602/1‑90 (MQ=255)
                                                gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:737507/1‑90 (MQ=255)
                                                gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:403960/1‑90 (MQ=255)
                                                gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCag                                   >  2:29247/1‑90 (MQ=255)
                                                    caATGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGa                               <  1:551414/90‑1 (MQ=255)
                                                      aTGCCGCCACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaac                             >  2:597845/1‑90 (MQ=255)
                                                            ccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGt                       <  2:164706/90‑1 (MQ=255)
                                                            ccACGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGt                       >  1:527635/1‑90 (MQ=255)
                                                              aCGTTCTTTCTGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTAc                     >  2:568517/1‑90 (MQ=255)
                                                                gTTCTTTCTGCGCCCCCCAGGCCACCATAGTACGACCACACGGCGCCCCCCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                   <  1:598350/90‑1 (MQ=255)
                                                                  tctttctGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAAc                 <  1:543499/90‑1 (MQ=255)
                                                                     ttctGCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc              <  2:577220/90‑1 (MQ=255)
                                                                         gCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCg          >  2:992346/1‑90 (MQ=255)
                                                                         gCGACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGAAACAACAGTACAAACATCACCg          >  2:655746/1‑90 (MQ=255)
                                                                          cGACCACCAGGGCACCATAGTACGACCACACGGCGGcc                                                             <  2:357839/38‑1 (MQ=39)
                                                                          cGACCACCAGGGCACCATAGTACGACCACACGGCGGcc                                                             >  1:357839/1‑38 (MQ=39)
                                                                           gACCACCAGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCa        >  2:681089/1‑90 (MQ=255)
                                                                               accaGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGc    >  2:524236/1‑90 (MQ=255)
                                                                                 caGGGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc  <  1:403960/90‑1 (MQ=255)
                                                                                    |                                                                                      
AGCAGCGGCGGAATACCACCACCGACGTTGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTC  >  NZ_CP009273/2345330‑2345500

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: