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,5760AC89.5% 52.6 / ‑0.6 19D67A (GAC→GCC) marRmultiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base C (9/8);  total (10/9)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.20e-01

CGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGC  >  NZ_CP009273/1613493‑1613658
                                                                                   |                                                                                  
cgcgGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTggg                                                                              <  1:224321/90‑1 (MQ=255)
 gcgGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGa                                                                             >  1:333366/1‑90 (MQ=255)
                 tgctctgctcTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGc                                                             <  1:158705/90‑1 (MQ=255)
                   ctctgctctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTg                                                           >  1:241355/1‑90 (MQ=255)
                   ctctgctctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTg                                                           >  2:130005/1‑90 (MQ=255)
                   ctctgctctATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTg                                                           >  2:375602/1‑90 (MQ=255)
                   ctctgctctATCCGCTGAGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTg                                                           >  2:190069/1‑90 (MQ=255)
                            aTCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTg                                                  >  1:385462/1‑90 (MQ=255)
                              ccGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGt                                                <  1:219452/90‑1 (MQ=255)
                               cGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTc                                               <  2:385462/90‑1 (MQ=255)
                                gCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCt                                              <  1:293622/90‑1 (MQ=255)
                                    agcggcgTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTaa                                          <  1:190069/89‑1 (MQ=255)
                                      cggcgTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTagag                                        >  2:459375/1‑90 (MQ=255)
                                        gcgTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGc                                      >  2:502/1‑90 (MQ=255)
                                         cgTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCt                                     <  2:180477/90‑1 (MQ=255)
                                                    ccGGTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAgg                          >  2:361290/1‑90 (MQ=255)
                                                      ggTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGa                                                         <  2:115075/57‑1 (MQ=255)
                                                      ggTTGAACTGAAAAAGGTATTGTCGGTCGCCCTGGGAGCACTGACCCGTATGCTGGa                                                         >  1:115075/1‑57 (MQ=255)
                                                                            tcggtcgCCCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAgcgc  <  2:403546/90‑1 (MQ=255)
                                                                                   |                                                                                  
CGCGGCACAGTTTAAGGTGCTCTGCTCTATCCGCTGCGCGGCGTGTATTACTCCGGTTGAACTGAAAAAGGTATTGTCGGTCGACCTGGGAGCACTGACCCGTATGCTGGATCGCCTGGTCTGTAAAGGCTGGGTGGAAAGGTTGCCGAACCCGAATGACAAGCGC  >  NZ_CP009273/1613493‑1613658

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