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,5760AC96.3% 88.9 / NA 28D67A (GAC→GCC) marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base A (1/0);  major base C (15/11);  minor base G (1/0);  total (17/11)
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

TTACCGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGC  >  NZ_CP009273/1613489‑1613661
                                                                                       |                                                                                     
ttACCGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGccc                                                                                     >  1:60601/1‑90 (MQ=255)
ttACCGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACGGAAAAAGGTATTGTCGGACGccc                                                                                     >  2:386566/1‑90 (MQ=255)
     gcgGCACGGTTTAAGGAGGTATGCTCTATCCGCCGCGCGGCGTGCATTACTCCGGTTGAACTGAAAAAGGTATTCTCGGTCGGCCTgggg                                                                                >  1:352443/1‑89 (MQ=255)
      cgGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAg                                                                               <  2:157641/90‑1 (MQ=255)
             gTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGAc                                                                        >  2:76803/1‑90 (MQ=255)
              tttAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGAcc                                                                       <  1:86599/90‑1 (MQ=255)
                tAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCg                                                                     <  2:216636/90‑1 (MQ=255)
                  aGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTa                                                                   >  1:70344/1‑90 (MQ=255)
                       ctctgctctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTg                                                              >  2:238522/1‑90 (MQ=255)
                       ctctgctctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTg                                                              >  2:119957/1‑90 (MQ=255)
                              ctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGcc                                                       <  1:119957/90‑1 (MQ=255)
                                  ccGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGt                                                   <  2:70344/90‑1 (MQ=255)
                                    gCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCt                                                 <  1:125178/90‑1 (MQ=255)
                                     cTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTg                                                >  1:382799/1‑90 (MQ=255)
                                              gtgtATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTg                                       <  2:352443/90‑1 (MQ=255)
                                                      cTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCAc                                                                            >  1:273548/1‑45 (MQ=255)
                                                      cTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCAc                                                                            <  2:273548/45‑1 (MQ=255)
                                                      cTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGaaa                               >  2:125671/1‑90 (MQ=255)
                                                      cTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGaaa                               >  2:531175/1‑90 (MQ=255)
                                                          ggTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGtt                           >  2:169706/1‑90 (MQ=255)
                                                          ggTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGtt                           >  2:169730/1‑90 (MQ=255)
                                                          ggTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGtt                           >  1:523792/1‑90 (MQ=255)
                                                             tGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGcc                        >  2:177626/1‑90 (MQ=255)
                                                               aaCTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTCGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGa                      <  2:198236/90‑1 (MQ=255)
                                                                aCTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGaa                     >  1:136504/1‑90 (MQ=255)
                                                                   gAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAAccc                  >  1:355426/1‑90 (MQ=255)
                                                                             ttGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAg        >  1:13122/1‑90 (MQ=255)
                                                                              tGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAgc       <  1:125671/90‑1 (MQ=255)
                                                                                  ggtcgCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCgcgg   <  2:382799/90‑1 (MQ=255)
                                                                                   gtcgCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCgcggc  >  2:10344/1‑90 (MQ=255)
                                                                                       |                                                                                     
TTACCGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGCGGC  >  NZ_CP009273/1613489‑1613661

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