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.95.8% 90.6 / ‑1.4 24coding (372/435 nt)marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base . (10/13);  total (11/13)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATG  >  NZ_CP009273/1613659‑1613833
                                                                                         |                                                                                     
ggcgTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACa                                                                                       >  1:155369/1‑90 (MQ=255)
      cTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaaacc                                                                                 >  1:233781/1‑88 (MQ=255)
      cTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaaaaacc                                                                                 >  1:129976/1‑88 (MQ=255)
                ttACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGa                                                                      >  2:654505/1‑90 (MQ=255)
                 tACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGa                                                                      >  2:426734/1‑90 (MQ=255)
                  accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACg                                                                    >  2:981672/1‑90 (MQ=255)
                    caccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGaa                                                                  <  2:119301/90‑1 (MQ=255)
                     accGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAg                                                                 <  1:594229/90‑1 (MQ=255)
                         gcggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGc                                                             <  2:122362/90‑1 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                            <  1:313312/90‑1 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                            <  1:654505/90‑1 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                            <  1:663062/90‑1 (MQ=255)
                            gcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAac                                                          <  2:696186/90‑1 (MQ=255)
                             cgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAaca                                                         <  1:146245/90‑1 (MQ=255)
                                       tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                               <  2:984001/90‑1 (MQ=255)
                                       tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                               <  1:58108/90‑1 (MQ=255)
                                       tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                               <  1:548506/90‑1 (MQ=255)
                                               atgtcATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                       >  2:209180/5‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGTCGAAAAGGCAAGCCTTGATTTTTTACTTAAg                                       >  1:633414/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                       >  1:988827/1‑90 (MQ=255)
                                                    aTCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGt                                  >  1:696587/1‑90 (MQ=255)
                                                                gCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAAc                      <  2:96206/90‑1 (MQ=255)
                                                                    ggACCTGCTCCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGAATTTGCTTAAGAAAGTCCTGCCGTAAACaaaa                  >  1:815229/1‑90 (MQ=255)
                                                                         tGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAagag             >  2:236565/1‑90 (MQ=255)
                                                                              cAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATg        >  1:622010/1‑90 (MQ=255)
                                                                                  aaTTAACAAAACCCTGACGGCGGACGAAGGGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa    <  1:809878/90‑1 (MQ=255)
                                                                                    ttAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATg  >  2:109915/1‑90 (MQ=255)
                                                                                         |                                                                                     
GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGATG  >  NZ_CP009273/1613659‑1613833

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