Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	3,576,634	A→G	R79R (CGA→CGG)	BW25113_RS17825 →	YrhA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	3,576,634	0	A	G	100.0%	61.9 / NA	19	R79R (CGA→CGG)	BW25113_RS17825	YrhA family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (15/4);  total (15/4)

AATGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAA  >  NZ_CP009273/3576548‑3576719                                                                                       |                                                                                      aaTGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGgtta                                                                                    >  1:1177187/1‑90 (MQ=255)  aTGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGgttat                                                                                   >  2:941124/1‑90 (MQ=255)  aTGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGgttat                                                                                   >  1:313078/1‑90 (MQ=255)    gAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTAttt                                                                                 >  2:1237270/1‑90 (MQ=255)       gTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGc                                                                              >  1:570768/1‑90 (MQ=255)         tAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCct                                                                            >  2:797681/1‑90 (MQ=255)           gTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCctct                                                                          >  1:45229/1‑90 (MQ=255)                cTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCtctc                                                                         <  1:782444/86‑1 (MQ=255)                cTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCtctc                                                                         >  2:782444/1‑86 (MQ=255)                 tttttAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTc                                                                    >  1:85913/1‑90 (MQ=255)                                            cTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTg                                         >  2:691759/1‑90 (MQ=255)                                                                      gAAGTTAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAaataat               >  2:756892/1‑90 (MQ=255)                                                                           tAATGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAAtgatga          >  2:1209969/1‑90 (MQ=255)                                                                             aTGGATTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGAtt        <  2:1213129/90‑1 (MQ=255)                                                                                 aTTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCAt    >  1:287043/1‑90 (MQ=255)                                                                                 aTTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCAt    >  2:542324/1‑90 (MQ=255)                                                                                 aTTACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCAt    >  2:7466/1‑90 (MQ=255)                                                                                  ttACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATa   <  1:64329/90‑1 (MQ=255)                                                                                   tACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATaa  >  1:270546/1‑90 (MQ=255)                                                                                   tACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATaa  <  2:292259/90‑1 (MQ=255)                                                                                   tACGGTTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATaa  >  1:57846/1‑90 (MQ=255)                                                                                   tACGGTTATTGAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATaa  >  1:1099823/1‑90 (MQ=255)                                                                                       |                                                                                      AATGAAGTTAGTATTCTTTTTAAGAATCAACCTGATTATCTTACTTTTTTAAGAGCAATGGATGGATTCGAAGTTAATGGATTACGATTATTTAGCCTCTCGATTCCAGAACCTTCAGTTAAAAACCTTTTTGCCGTAAATGAATTTTATAGAAATAATGATGATTTCATAA  >  NZ_CP009273/3576548‑3576719

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