Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 2,621,216 A→G I361T (ATC→ACC)  pdeF ← cyclic‑guanylate‑specific phosphodiesterase PdeF

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092732,621,2160AG100.0% 58.0 / NA 17I361T (ATC→ACC) pdeFcyclic‑guanylate‑specific phosphodiesterase PdeF
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (7/10);  total (7/10)

AAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTC  >  NZ_CP009273/2621129‑2621301
                                                                                       |                                                                                     
aaGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTg                                                                                     <  1:356912/90‑1 (MQ=255)
            cAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACaa                                                                         >  1:224970/1‑90 (MQ=255)
                 agaAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACAAAgcgc                                                                    <  1:384863/90‑1 (MQ=255)
                 agaAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACAAAgcgc                                                                    <  1:195869/90‑1 (MQ=255)
                       tttttGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACAAAGCGCAGACCa                                                              <  1:367205/90‑1 (MQ=255)
                                          cGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAAcgcg                                           >  2:508226/1‑90 (MQ=255)
                                          cGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAAcgcg                                           <  2:133740/90‑1 (MQ=255)
                                           gcagcaTGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAAcgcgc                                          >  1:366484/1‑90 (MQ=255)
                                           gcagcaTGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAAcgcgc                                          >  2:258921/1‑90 (MQ=255)
                                            cagcaTGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAAcgcgcg                                         <  2:421865/90‑1 (MQ=255)
                                                   atgCCATAGTTCTTAACCAGCATTTCCATGCCAGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAgg                                  <  2:107409/90‑1 (MQ=255)
                                                                     aGCATTTCCATGCCAGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGa                >  1:458371/1‑90 (MQ=255)
                                                                      gCATTTCCATGCCAGGGGTGCGTAAATAACAAAgcgc                                                                    <  2:218096/37‑1 (MQ=38)
                                                                      gCATTTCCATGCCAGGGGTGCGTAAATAACAAAgcgc                                                                    >  1:218096/1‑37 (MQ=38)
                                                                              aTGCCAGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCTGGGCGCGAACATTTGGCAGATGAaccacc       <  1:508226/90‑1 (MQ=255)
                                                                                  cAGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGt   <  2:116249/90‑1 (MQ=255)
                                                                                   aGGGGTGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTc  >  1:430958/1‑90 (MQ=255)
                                                                                       |                                                                                     
AAGGGTGACAGCCAGTGAGAAAGTTTTTGCTTGTATTGAATCCGCAGCATGATGCCATAGTTCTTAACCAGCATTTCCATGCCAGGGATGCGTAAATAACAAAGCGCAGACCAGGGGGCATCACGCAACGCGCGATTCAGGGCGCGAACATTTGGCAGATGAACCACCGGGTC  >  NZ_CP009273/2621129‑2621301

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