Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,613,576 A→C D67A (GAC→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,5760AC87.5% 44.5 / ‑0.1 16D67A (GAC→GCC) marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (8/6);  total (9/7)
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.31e-01

TTACCGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGG  >  NZ_CP009273/1613489‑1613660
                                                                                       |                                                                                    
ttACCGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGccc                                                                                    >  1:447702/1‑90 (MQ=255)
  aCCGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTg                                                                                  <  2:447702/90‑1 (MQ=255)
     gcgGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGa                                                                               >  2:593267/1‑90 (MQ=255)
                      gctctgctctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCt                                                              >  2:356704/1‑90 (MQ=255)
                       ctctgctctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTg                                                             <  1:480568/90‑1 (MQ=255)
                       ctctgctctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTg                                                             >  2:295339/1‑90 (MQ=255)
                              ctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGcc                                                      <  2:31345/90‑1 (MQ=255)
                                   cGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTc                                                 <  1:97799/90‑1 (MQ=255)
                                     cTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGGATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTg                                               >  2:311530/1‑90 (MQ=255)
                                            gcgTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGc                                        <  1:30303/90‑1 (MQ=255)
                                                      cTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGt                                                  >  1:88405/1‑70 (MQ=255)
                                                      cTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGt                                                  <  2:88405/70‑1 (MQ=255)
                                                      cTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGaaa                              >  1:559099/1‑90 (MQ=255)
                                                        ccGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAgg                            >  2:533902/1‑90 (MQ=255)
                                                          ggTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGtt                          >  1:125842/1‑90 (MQ=255)
                                                               aaCTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGa                     <  2:559099/90‑1 (MQ=255)
                                                                                  ggtcgCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCgcgg  >  1:413844/1‑90 (MQ=255)
                                                                                       |                                                                                    
TTACCGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGG  >  NZ_CP009273/1613489‑1613660

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: