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.96.1% 194.4 / 3.4 51coding (372/435 nt)marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base . (32/17);  total (33/18)
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.61e-01

GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGA  >  NZ_CP009273/1613659‑1613831
                                                                                         |                                                                                   
ggcgTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACa                                                                                     >  1:135172/1‑90 (MQ=255)
             aaCTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGAcggc                                                                       >  2:1188709/1‑90 (MQ=255)
                 tACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGAc                                                                   <  1:1405539/90‑1 (MQ=255)
                  accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGGCCTGCACCAAGAATTAACAAAAACCTGACGGCGGACg                                                                  >  2:1439029/1‑90 (MQ=255)
                  accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACgacggacg                                                                  >  1:369765/1‑90 (MQ=255)
                  accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACg                                                                  >  2:94791/1‑90 (MQ=255)
                     accGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAg                                                               <  1:1150810/90‑1 (MQ=255)
                     accGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAg                                                               <  2:544293/90‑1 (MQ=255)
                      ccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGt                                                              >  2:1112196/1‑90 (MQ=255)
                        ggcggcGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTgg                                                            >  1:1339301/1‑90 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                          <  2:542480/90‑1 (MQ=255)
                          cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCa                                                          <  2:1438923/90‑1 (MQ=255)
                            gcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAac                                                        <  1:757241/90‑1 (MQ=255)
                             cgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAaca                                                       <  1:527197/90‑1 (MQ=255)
                                 gcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTg                                                   <  2:941432/90‑1 (MQ=255)
                                       tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                             <  1:1329396/90‑1 (MQ=255)
                                           aaCAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCt                                         >  2:800259/1‑90 (MQ=255)
                                              aaTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTaa                                      >  2:1477138/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGAGGAAGTGGCAACACTTGAGGATTTGCTTaaa                                     >  1:1021218/1‑89 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                     >  2:387360/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                     >  1:1343130/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAg                                     >  1:585395/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGGGGCAACACTTGAGTATTTGCTTAAg                                     >  2:1077995/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGGGGCAACACTTGAGGATTTGCTTAAg                                     >  1:919854/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACAAAGTGGCAACACTTGAGGATTTGCTTAAg                                     >  1:1425251/1‑90 (MQ=255)
                                               aTGCCATCAATTAGTTGGACAGGACCTGCACCAAGAATTAACAAAAACCTGACGGGGGGCTAAGTGGCAACACTTGAGTATTTGCttttt                                     >  1:1549383/1‑87 (MQ=255)
                                                tGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGa                                    >  1:1246650/1‑90 (MQ=255)
                                                    aTCAATTAGTTGGCCAGGACCTGAACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGt                                >  1:538465/1‑90 (MQ=255)
                                                              tGGCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTaa                      >  1:58280/1‑90 (MQ=255)
                                                               ggCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAac                                                        <  2:713079/56‑1 (MQ=255)
                                                               ggCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAac                                                        >  1:713079/1‑56 (MQ=255)
                                                               ggCCAGGACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTaaa                     >  1:818617/1‑90 (MQ=255)
                                                                     gACCTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACaaaaa               >  1:45767/1‑90 (MQ=255)
                                                                      gccTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGCCCTGCCGTAAACaaaaaa              <  1:1439029/89‑1 (MQ=255)
                                                                        cTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAaga            <  1:252166/90‑1 (MQ=255)
                                                                        cTGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAAAAAGTCCTGCCGTAAACAAAAAaga            <  1:1112196/90‑1 (MQ=255)
                                                                         tGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAagag           >  2:105841/1‑90 (MQ=255)
                                                                         tGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAagag           >  1:1230956/1‑90 (MQ=255)
                                                                         tGCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAagag           >  1:1196671/1‑90 (MQ=255)
                                                                          gCACCAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAgg          >  1:1112083/1‑90 (MQ=255)
                                                                              cAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATg      >  2:1429267/1‑90 (MQ=255)
                                                                              cAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATg      >  2:379062/1‑90 (MQ=255)
                                                                              cAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATg      >  1:726611/1‑90 (MQ=255)
                                                                              cAAGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATg      >  1:229789/1‑90 (MQ=255)
                                                                               aaGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGa     >  2:1441391/1‑90 (MQ=255)
                                                                               aaGAATTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGa     >  1:501413/1‑90 (MQ=255)
                                                                                  aaTTAACAAAACCCTGCCGGCGGACGAAGTGGAAACACTTGAGTATTTGCTTAAAAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  <  2:1196671/90‑1 (MQ=255)
                                                                                  aaTTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  <  2:1493732/90‑1 (MQ=255)
                                                                                  aaTTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  <  2:1246650/90‑1 (MQ=255)
                                                                                  aaTTAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTAAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  <  1:1441391/90‑1 (MQ=255)
                                                                                    ttAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  <  1:926693/88‑1 (MQ=255)
                                                                                    ttAACAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  >  2:926693/1‑88 (MQ=255)
                                                                                         |                                                                                   
GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGA  >  NZ_CP009273/1613659‑1613831

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