Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,613,623 T→G W83G (TGG→GGG)  marR → multiple antibiotic resistance transcriptional regulator MarR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NZ_CP0092731,613,6230TG95.7% 82.3 / ‑3.2 23W83G (TGG→GGG) marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (14/8);  total (15/8)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.78e-01

TATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAAT  >  NZ_CP009273/1613538‑1613707
                                                                                     |                                                                                    
tATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCggggg                                                                                  >  1:323865/1‑89 (MQ=255)
       ccGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAgg                                                                           >  2:190220/1‑90 (MQ=255)
       ccGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGGGGAAAgg                                                                           >  1:187820/1‑90 (MQ=255)
       ccGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGGGGAAAgg                                                                           >  1:259973/1‑90 (MQ=255)
         ggTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAgggt                                                                         >  2:98396/1‑88 (MQ=255)
         ggTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGtt                                                                         >  2:545676/1‑90 (MQ=255)
         ggTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGtt                                                                         >  1:716376/1‑90 (MQ=255)
         ggTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGGGGGAAGGtt                                                                         >  1:250245/1‑90 (MQ=255)
               aCTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGGGGAAAGGTTGCCGaa                                                                   >  1:212382/1‑90 (MQ=255)
                               tcggtcgACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAgcgc                                                   <  2:716376/90‑1 (MQ=255)
                                 ggtcgACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCgcgg                                                 <  1:545676/90‑1 (MQ=255)
                                     gACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTa                                             <  2:237795/90‑1 (MQ=255)
                                     gACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTa                                             <  2:45229/90‑1 (MQ=255)
                                     gACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTa                                             <  1:302748/90‑1 (MQ=255)
                                     gACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGGGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTa                                             >  1:45229/1‑90 (MQ=255)
                                                aCTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGCGGGAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACt                                  >  1:243788/1‑90 (MQ=255)
                                                  tGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTa                                >  2:200792/1‑90 (MQ=255)
                                                    aCCCGTATGCTGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTacc                              <  2:620593/90‑1 (MQ=255)
                                                              tGGATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGgcggcg                    <  1:646389/90‑1 (MQ=255)
                                                               ggATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGgcggcgc                   >  1:68262/1‑90 (MQ=255)
                                                               ggATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGgcggcgc                   >  2:111457/1‑90 (MQ=255)
                                                               ggATCGCCTGGTCTGTAAAGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGGAAAACTTACCACCGgcggcgc                   >  1:738306/1‑90 (MQ=255)
                                                                         gTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATAtgt         >  1:151735/1‑90 (MQ=255)
                                                                                aaGGCGGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAAt  <  1:794829/90‑1 (MQ=255)
                                                                                     |                                                                                    
TATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAAT  >  NZ_CP009273/1613538‑1613707

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