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.4 / NA 21Y518H (TAT→CAT) dosPoxygen‑sensing cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (5/16);  total (5/16)

CTGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGG  >  NZ_CP009273/1558352‑1558523
                                                                                       |                                                                                    
cTGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGAt                                                                                    >  2:100499/1‑90 (MQ=255)
 tGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc                                                                                   <  1:79501/90‑1 (MQ=255)
      cTAAAACCAAACGCTTTTTTCCCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattg                                                                              <  1:627239/90‑1 (MQ=255)
            ccAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATtgt                                                                        >  2:446966/1‑90 (MQ=255)
                ccGCTTTTTTCCTTTTTCGTTCATCGCGGGGCTGAAGACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                    <  2:135714/89‑1 (MQ=255)
                ccGCTCTTTTCCCTTTTCTTTCATCGCCGGGCTGAAGACCTCCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                    <  1:401208/89‑1 (MQ=255)
                aCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                    <  1:176053/90‑1 (MQ=255)
                     tttttCCCTTTTCTTTCATCGCGGGGTTAAAAACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTgg                                                               <  2:636392/89‑1 (MQ=255)
                     tttttCCCTTTTCTTTCACCGCGGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTgg                                                               <  2:251618/89‑1 (MQ=255)
                       tttACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGgag                                                             <  1:446966/90‑1 (MQ=255)
                             aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAg                                                       >  1:636392/1‑90 (MQ=255)
                                               cTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                     <  1:32370/90‑1 (MQ=255)
                                                      aCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGt                              <  1:226307/90‑1 (MQ=255)
                                                           cAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTg                         <  2:23665/90‑1 (MQ=255)
                                                           cAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTg                         <  1:163503/90‑1 (MQ=255)
                                                                   aCCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAAt                 <  1:578219/90‑1 (MQ=255)
                                                                         aTTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCaa           <  1:245283/90‑1 (MQ=255)
                                                                            cTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGt        <  1:530226/90‑1 (MQ=255)
                                                                                cGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGc    <  2:621095/90‑1 (MQ=255)
                                                                                 gAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCg   >  1:512824/1‑90 (MQ=255)
                                                                                  aaTATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCgg  >  1:584829/1‑90 (MQ=255)
                                                                                       |                                                                                    
CTGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAAGCGG  >  NZ_CP009273/1558352‑1558523

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: