Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,613,723 G→C G116A (GGC→GCC)  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,7230GC90.0% 66.5 / 1.0 20G116A (GGC→GCC) marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base C (7/11);  total (8/12)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.07e-01

TGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTA  >  NZ_CP009273/1613636‑1613810
                                                                                       |                                                                                       
tGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGccc                                                                                       >  1:637231/1‑90 (MQ=255)
 gCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGCCCa                                                                                      <  2:616683/90‑1 (MQ=255)
      aCCCGAATGACAAGCGCGGCCTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGAcc                                                                                 <  2:890646/90‑1 (MQ=255)
       cccGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCt                                                                                <  2:822552/90‑1 (MQ=255)
             tGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCa                                                                          <  2:524095/90‑1 (MQ=255)
              gACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCaa                                                                         <  1:321071/90‑1 (MQ=255)
                                       ttACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGAcggcgg                                                >  1:784510/1‑90 (MQ=255)
                                       ttACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGAcggcgg                                                >  1:813607/1‑90 (MQ=255)
                                         accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCa                                                                          <  1:335433/62‑1 (MQ=255)
                                         accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCa                                                                          >  2:335433/1‑62 (MQ=255)
                                                  ggcgcggcAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCa                                     <  2:165079/90‑1 (MQ=255)
                                                  ggcgcggcAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCa                                     <  1:510056/90‑1 (MQ=255)
                                                    cgcggcAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAac                                   <  2:300237/90‑1 (MQ=255)
                                                    cgcggcAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAac                                   <  2:196163/90‑1 (MQ=255)
                                                    cgcggcAATATGTGAACAATGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAac                                   <  2:893414/90‑1 (MQ=255)
                                                                  aaCAATGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGc                     >  2:564459/1‑90 (MQ=255)
                                                                    cAATGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCtt                   <  2:224797/90‑1 (MQ=255)
                                                                     aaTGCCATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTa                  >  2:783646/1‑90 (MQ=255)
                                                                         ccATCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGaa              >  2:755598/1‑90 (MQ=255)
                                                                           aTCAATTAGTTGCCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAg            >  1:76011/1‑90 (MQ=255)
                                                                                     tGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTa  >  1:33467/1‑90 (MQ=255)
                                                                                       |                                                                                       
TGCCGAACCCGAATGACAAGCGCGGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTA  >  NZ_CP009273/1613636‑1613810

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