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,4390AG94.7% 52.9 / ‑3.9 19Y518H (TAT→CAT) dosPoxygen‑sensing cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (7/11);  total (8/11)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

TGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAA  >  NZ_CP009273/1558353‑1558519
                                                                                      |                                                                                
tGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc                                                                               <  1:259280/90‑1 (MQ=255)
 gcgcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATcc                                                                              >  2:464403/1‑90 (MQ=255)
     caaaacCCAACCCCTTTTTTCCCTTTTCTTTCATCGCCGGGCTGAAGACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattg                                                                          <  2:405655/88‑1 (MQ=255)
     cTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCGGGGCTGAAGACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattg                                                                          <  1:659778/90‑1 (MQ=255)
          aCCAACCGCTTTTTTCCCTTTTCTTTCATCGCCGGGCTGAAAAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgcattg                                                                     <  2:566846/90‑1 (MQ=255)
               aCGCTTTTTTCCCTTTTCTTTCATCGCCGGGCTAAAAAACTCCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                <  1:525472/90‑1 (MQ=255)
               aCGCTTTTTTCCCTTTTCGTTCATCGCCGGGCTGAAAACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                <  2:708699/90‑1 (MQ=255)
                            aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAg                                                   >  1:293455/1‑90 (MQ=255)
                            aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAg                                                   >  1:667557/1‑90 (MQ=255)
                              ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTa                                                 >  2:74676/1‑90 (MQ=255)
                                    caccgcgGGGCTGAAGACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                                           <  1:741254/87‑1 (MQ=255)
                                    caccGCCGGGTTAAAAAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                                           <  1:776612/87‑1 (MQ=255)
                                      tCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGt                                         >  1:806116/1‑90 (MQ=255)
                                      tCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGt                                         >  1:657522/1‑90 (MQ=255)
                                              ttGAAGACCTCCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                 <  2:806116/89‑1 (MQ=255)
                                                gaagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAg                               >  2:141035/1‑90 (MQ=255)
                                                            gCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGAt                   <  2:360354/90‑1 (MQ=255)
                                                                          tCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGACGCCAATACTCAAgg     <  1:69758/90‑1 (MQ=255)
                                                                             tGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCTAATACTCAAGGTaa  >  2:644050/1‑90 (MQ=255)
                                                                                      |                                                                                
TGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAA  >  NZ_CP009273/1558353‑1558519

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