Predicted mutation
evidence seq id position mutation annotation gene description
RA NZ_CP009273 1,613,753 (A)6→5 coding (377/435 nt) 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,7480A.93.9% 121.6 / 2.0 33coding (372/435 nt)marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (20/11);  total (21/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 = 1.00e+00

GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGA  >  NZ_CP009273/1613659‑1613831
                                                                                         |                                                                                   
ggcgTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACa                                                                                     >  1:1158209/1‑90 (MQ=255)
ggcgTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACa                                                                                     >  1:309019/1‑90 (MQ=255)
           aaaaCTTACCACCGGCGGGGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGAc                                                                          >  1:1070764/1‑90 (MQ=255)
                ttACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGa                                                                    >  2:728983/1‑90 (MQ=255)
                ttACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGa                                                                    >  2:269418/1‑90 (MQ=255)
                  accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACg                                                                  >  1:541313/1‑90 (MQ=255)
                         gcggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGc                                                           >  2:8705/1‑90 (MQ=255)
                         gcggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGc                                                           <  2:309019/90‑1 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGCCCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                          <  1:955252/90‑1 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                          <  2:656493/90‑1 (MQ=255)
                               cggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACt                                                     <  1:483863/90‑1 (MQ=255)
                                              aaTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTaa                                      >  1:222195/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                     >  2:801722/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                     >  2:612967/1‑90 (MQ=255)
                                                tGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGa                                    >  2:113681/1‑90 (MQ=255)
                                                tGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGa                                    >  1:824468/1‑90 (MQ=255)
                                                 gCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGaa                                   >  1:206197/1‑90 (MQ=255)
                                                         ttAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGc                           >  1:957959/1‑90 (MQ=255)
                                                         ttAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGc                           >  1:11707/1‑90 (MQ=255)
                                                             ttGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGGGGACGAAGTGGCAACACTTGAGTTTTTGCTTAAGAAAGTCCTGCCGTa                       <  1:350661/90‑1 (MQ=255)
                                                                gCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAAc                    >  2:1020895/1‑90 (MQ=255)
                                                                gCCAGGACCTGCACCAAGAATTAAAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAAc                    <  1:65723/90‑1 (MQ=255)
                                                                      aCCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACaaaaaa              <  2:852544/90‑1 (MQ=255)
                                                                      aCCTGCACCAAGAATTAACAAAAACCGGACGGGGGACAAAGGGGCAACACTTGAGTATTTGCTTAAGAAAGCCCTGCCGTAAACaaaaaa              <  1:127824/90‑1 (MQ=255)
                                                                        cTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAaga            <  2:957959/90‑1 (MQ=255)
                                                                        cTGCACAAAGAATTAACAAAACCCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAaga            <  2:257926/90‑1 (MQ=255)
                                                                         tGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAagag           >  2:892395/1‑90 (MQ=255)
                                                                          gCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAgg          >  2:710078/1‑90 (MQ=255)
                                                                          gCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAgg          >  2:595955/1‑90 (MQ=255)
                                                                              cAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATg      >  1:217901/1‑90 (MQ=255)
                                                                               aaGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGa     >  2:670100/1‑90 (MQ=255)
                                                                               aaGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGa     >  2:862831/1‑90 (MQ=255)
                                                                               aaGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGa     >  1:986651/1‑90 (MQ=255)
                                                                                  aaTTAACAAAACCCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  <  2:541313/90‑1 (MQ=255)
                                                                                  aaTTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  <  2:810503/90‑1 (MQ=255)
                                                                                         |                                                                                   
GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGA  >  NZ_CP009273/1613659‑1613831

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