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,4390AG80.0% 14.0 / 0.7 10Y518H (TAT→CAT) dosPoxygen‑sensing cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (3/5);  total (4/6)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

GCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATAC  >  NZ_CP009273/1558354‑1558510
                                                                                     |                                                                       
gcgcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATcc                                                                     >  1:5892/1‑90 (MQ=255)
    cTAAACCCAACCGCTTTTTTCCCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattg                                                                 <  1:304080/90‑1 (MQ=255)
         aCCAACCGCTTTTTTCCCTTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgcattg                                                            <  2:77842/90‑1 (MQ=255)
         aCCAAACGCTTTTTTCCCTTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgcattg                                                            <  1:377849/90‑1 (MQ=255)
                   ctTTTCCCTTTTCGTTCACCGCGGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTgg                                                  <  2:80312/90‑1 (MQ=255)
                           aTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAg                                          >  1:93633/1‑90 (MQ=255)
                             ttCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTa                                        >  1:366764/1‑90 (MQ=255)
                                                                tcccGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCaa     <  1:208016/88‑1 (MQ=255)
                                                                tACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCaa     <  2:384753/90‑1 (MQ=255)
                                                                   cGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATAc  >  1:381203/1‑90 (MQ=255)
                                                                                     |                                                                       
GCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATAC  >  NZ_CP009273/1558354‑1558510

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