Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,558,439 A→G Y518H (TAT→CAT)  dosP ← oxygen‑sensing cyclic‑di‑GMP phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,558,4390AG100.0% 67.2 / NA 19Y518H (TAT→CAT) dosPoxygen‑sensing cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/14);  total (5/14)

TGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGG  >  NZ_CP009273/1558353‑1558524
                                                                                      |                                                                                     
tGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc                                                                                    <  2:548055/90‑1 (MQ=255)
     cTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattg                                                                               <  2:210845/90‑1 (MQ=255)
            cAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATtgtg                                                                        >  1:530132/1‑90 (MQ=255)
             aaaCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcat                                                                                 >  2:78940/1‑80 (MQ=255)
             aaaCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcat                                                                                 <  1:78940/80‑1 (MQ=255)
               ccGCTCTTTTCCCTTTTCGTTCATCGCGGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                     <  1:501557/89‑1 (MQ=255)
               aCGCTCTTTTCCCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                     <  2:54344/90‑1 (MQ=255)
                    tttttCCTATTTCTTTCATCGCGGGCATAAAAAACTGCCGGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTgg                                                                <  1:543291/89‑1 (MQ=255)
                              ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTa                                                      >  2:435281/1‑90 (MQ=255)
                                   tCACCGCGGGGCTAAAGACCTCCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgc                                                 <  1:406351/90‑1 (MQ=255)
                                    cATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                                                <  2:83692/90‑1 (MQ=255)
                                              cTGAAGAACTGCCAGCCGTTACCGCCAGTCTTGCGAATATGATCCATTGCTTTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                      <  2:552892/90‑1 (MQ=255)
                                                      cTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                      <  2:362800/82‑1 (MQ=255)
                                                      cTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                      >  1:362800/1‑82 (MQ=255)
                                                          cAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTg                          <  2:530132/90‑1 (MQ=255)
                                                            gCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGAt                        <  2:216415/90‑1 (MQ=255)
                                                                ttACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCa                    <  2:61990/90‑1 (MQ=255)
                                                                         ttCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAg           >  2:239818/1‑90 (MQ=255)
                                                                                  atatGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCggg  <  1:545057/90‑1 (MQ=255)
                                                                                      |                                                                                     
TGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGGG  >  NZ_CP009273/1558353‑1558524

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