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,4390AG97.1% 105.5 / ‑3.5 34Y518H (TAT→CAT) dosPoxygen‑sensing cyclic‑di‑GMP phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (14/19);  total (14/20)
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.93e-01

TGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTAA  >  NZ_CP009273/1558353‑1558519
                                                                                      |                                                                                
tGCGCCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc                                                                               >  2:1141416/1‑90 (MQ=255)
tGCCCCAAAACCCAACCGCTTTTTTCCCTTTTCTTTCACCCCCGGGCTAAAAAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATc                                                                               <  1:1141416/90‑1 (MQ=255)
 gcgcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATcc                                                                              >  2:250349/1‑90 (MQ=255)
 gcgcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATcc                                                                              >  1:469429/1‑90 (MQ=255)
 gcgcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATcc                                                                              >  1:52389/1‑90 (MQ=255)
 gcgcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATcc                                                                              >  1:544590/1‑90 (MQ=255)
 gcgcCTAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATcc                                                                              >  1:86064/1‑90 (MQ=255)
      tAAAACCAAACGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgc                                                                         >  2:89590/1‑90 (MQ=255)
          aCCAACCGCTTTTTTCCCTTTTCTTTCATCGCCGGGCTAAAAACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCcattgcattg                                                                     <  1:562590/90‑1 (MQ=255)
               aCGCTTTTTTCCCTTTTCTTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                <  2:892472/90‑1 (MQ=255)
               aCGCTCTTTTACCATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGc                                                                <  2:6444/90‑1 (MQ=255)
                  cttttttCCCTTTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGt                                                             <  2:1043451/87‑1 (MQ=255)
                        tcccATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGagag                                                       <  2:1080033/88‑1 (MQ=255)
                           cATTTCGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCaa                                                    <  2:52389/90‑1 (MQ=255)
                                cGTTCATCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAAt                                               >  2:464920/1‑90 (MQ=255)
                                    caccGCCGGGTTGAAAACCTCCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                                           <  1:886132/87‑1 (MQ=255)
                                    caccGCCGGGCTAAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                                           <  2:747402/87‑1 (MQ=255)
                                    cATCGCGGGGCTGAAAACCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATcgcg                                           <  1:386000/90‑1 (MQ=255)
                                      tCGCCGGGCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGt                                         >  1:753387/1‑90 (MQ=255)
                                           gggCTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTa                                    >  2:937919/1‑90 (MQ=255)
                                              ctaaagacCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                 <  1:237305/87‑1 (MQ=255)
                                              ctaaagaaCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                 <  2:1069988/87‑1 (MQ=255)
                                              cggaagacCTCCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                 <  1:723645/88‑1 (MQ=255)
                                              cTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                 <  2:233964/90‑1 (MQ=255)
                                              cTGAAGAACTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTAccc                                 <  1:451231/90‑1 (MQ=255)
                                                     aCTGCCAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGt                          >  1:922227/1‑90 (MQ=255)
                                                          cAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTg                     >  2:802069/1‑90 (MQ=255)
                                                          cAGCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTg                     >  1:304705/1‑90 (MQ=255)
                                                            gCCGTTACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGAt                   <  1:652943/90‑1 (MQ=255)
                                                            gCCGTTACCGCCATTCTTGCGAATATAATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGAt                   <  2:209564/90‑1 (MQ=255)
                                                                ttACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCa               >  2:211774/1‑90 (MQ=255)
                                                                ttACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCa               <  1:937919/90‑1 (MQ=255)
                                                                 tACCGCCATTCTTGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCaa              <  2:69479/90‑1 (MQ=255)
                                                                             tGCGAATATGATCCATTGCATTGTGAGCAGTGGAGAGCAAGTAATCGCGGTTTTTACCCAGGTCGTAGCTGATGCCAATACTCAAGGTaa  <  1:210332/90‑1 (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 < 19 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: