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_CP009273	1,613,748	0	A	.	88.5%	76.2 / 10.0	26	coding (372/435 nt)	marR	multiple antibiotic resistance transcriptional regulator MarR
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base . (7/16);  total (9/17)
Fisher's exact test for biased strand distribution p-value = 2.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.36e-01

GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGA  >  NZ_CP009273/1613659‑1613831                                                                                          |                                                                                    ggcgTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGGAAACAATGCCATGAATTAGGTGGCCAGGACCTGCACCAAGAATTAACa                                                                                     >  2:67036/1‑90 (MQ=255)  gcgTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACaa                                                                                    >  2:618286/1‑90 (MQ=255)              aaCTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGAcggc                                                                       >  2:791360/1‑90 (MQ=255)                   accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAAGCTGACGGCGGACg                                                                  >  2:731086/1‑90 (MQ=255)                   accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACg                                                                  >  2:675578/1‑90 (MQ=255)                   accaccGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACg                                                                  >  2:294056/1‑90 (MQ=255)                           cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCCGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCa                                                          <  1:633431/90‑1 (MQ=255)                           cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCa                                                          <  1:195044/90‑1 (MQ=255)                           cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCa                                                          <  1:480859/90‑1 (MQ=255)                           cggcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCa                                                          <  1:696721/90‑1 (MQ=255)                             gcgcggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAac                                                        <  1:88280/90‑1 (MQ=255)                                cggcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACt                                                     <  2:817754/90‑1 (MQ=255)                                  gcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTg                                                   <  2:242619/90‑1 (MQ=255)                                  gcAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTg                                                   <  1:646885/90‑1 (MQ=255)                                       atGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAt                                              >  2:207376/1‑90 (MQ=255)                                        tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                             <  2:294365/90‑1 (MQ=255)                                        tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                             <  2:47185/90‑1 (MQ=255)                                        tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                             <  2:601180/90‑1 (MQ=255)                                        tgtgAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                             <  2:657611/90‑1 (MQ=255)                                         gtgGACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTAtt                                             <  2:813731/90‑1 (MQ=255)                                                     aTCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAg                                 >  1:132258/1‑90 (MQ=255)                                                     aTCAATTAGTTGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGGGGCAACACTTGAGAATTTGCTTAAGAAAGt                                >  2:735967/1‑90 (MQ=255)                                                               tGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTTCCGTaa                      >  2:230404/1‑90 (MQ=255)                                                               tGGCCAGGACCTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTaa                      <  1:287943/90‑1 (MQ=255)                                                                   cAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACa                   >  2:110976/1‑90 (MQ=255)                                                                         cTGCACCAAGAATTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAaga            <  2:800167/90‑1 (MQ=255)                                                                                   aaTTAAC‑AAAACCCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  <  1:680054/90‑1 (MQ=255)                                                                                   aaTTAAC‑AAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGa  <  1:728862/90‑1 (MQ=255)                                                                                          |                                                                                    GGCGTACTGGTAAAACTTACCACCGGCGGCGCGGCAATATGTGAACAATGCCATCAATTAGTTGGCCAGGACCTGCACCAAGAATTAACAAAAAACCTGACGGCGGACGAAGTGGCAACACTTGAGTATTTGCTTAAGAAAGTCCTGCCGTAAACAAAAAAGAGGTATGACGA  >  NZ_CP009273/1613659‑1613831

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