Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NZ_CP009273	4,171,725	G→C	D185H (GAC→CAC)	rpoB →	DNA‑directed RNA polymerase subunit beta

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NZ_CP009273	4,171,725	0	G	C	100.0%	128.5 / NA	40	D185H (GAC→CAC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (22/18);  total (22/18)

AGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACCATCATTCTGCGCGCC  >  NZ_CP009273/4171655‑4171808                                                                       |                                                                                     aGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGaa                                                                   <  1:76216/90‑1 (MQ=255)         cccACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           <  1:173705/90‑1 (MQ=255)           cACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTg                                                         >  2:492813/1‑90 (MQ=255)            aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               <  2:519834/83‑1 (MQ=255)            aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               >  1:519834/1‑83 (MQ=255)            aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  1:67576/1‑90 (MQ=255)            aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        <  1:252955/90‑1 (MQ=255)                               gTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCg                                     >  2:462939/1‑90 (MQ=255)                               gTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCg                                     >  2:411626/1‑90 (MQ=255)                                    aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           <  2:139387/63‑1 (MQ=255)                                    aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           >  1:139387/1‑63 (MQ=255)                                    aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTa                               >  1:217132/1‑90 (MQ=255)                                    aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTa                               >  1:581847/1‑90 (MQ=255)                                    aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTa                               >  2:513364/1‑90 (MQ=255)                                      gcgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTaaa                             >  1:370805/1‑90 (MQ=255)                                       cgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAAc                            <  1:411626/90‑1 (MQ=255)                                       cgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAAc                            <  1:462939/90‑1 (MQ=255)                                       cgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAAc                            <  1:55378/90‑1 (MQ=255)                                         cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    >  2:144374/1‑65 (MQ=255)                                         cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    <  1:144374/65‑1 (MQ=255)                                         cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAActg                          >  2:197184/1‑90 (MQ=255)                                         cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAActg                          >  1:582058/1‑90 (MQ=255)                                         cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAActg                          >  1:16918/1‑90 (MQ=255)                                         cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAActg                          >  2:542603/1‑90 (MQ=255)                                            atcatcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAActgcct                       <  2:243544/90‑1 (MQ=255)                                             tcatcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAActgcctg                      <  2:67576/90‑1 (MQ=255)                                               atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCg                    >  2:233877/1‑90 (MQ=255)                                               atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCg                    <  1:43487/90‑1 (MQ=255)                                                   cTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACca                >  1:115506/1‑90 (MQ=255)                                                    ttACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACcat               <  1:233877/90‑1 (MQ=255)                                                     tACCGTGGTTCCTGGCTGCACTTCTAATTCGATCCGAAGGACAACCTGTTcgt                                                    >  1:559550/1‑53 (MQ=255)                                                     tACCGTGGTTCCTGGCTGCACTTCTAATTCGATCCGAAGGACAACCTGTTcgt                                                    <  2:559550/53‑1 (MQ=255)                                                     tACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACcatc              >  2:52014/1‑90 (MQ=255)                                                         gTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCCGTAAACTGCCTGCGACCATCAtt           >  1:216556/1‑90 (MQ=255)                                                         gTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACCATCATTc          >  1:111423/1‑90 (MQ=255)                                                             ttCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACCATCATTCTgcg      <  1:11851/90‑1 (MQ=255)                                                             ttCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACCATCATTCTgcg      <  2:582058/90‑1 (MQ=255)                                                              tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACCATCATTCTgcgc     <  1:542603/90‑1 (MQ=255)                                                              tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACCATCATTCTgcgc     >  2:269373/1‑90 (MQ=255)                                                              tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACCATCATTCTgcgc     <  2:360266/90‑1 (MQ=255)                                                                 tggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACCATCATTCTGCGCGcc  >  1:181353/1‑90 (MQ=255)                                                                       |                                                                                     AGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCG‑CCGTAAACTGCCTGCGACCATCATTCTGCGCGCC  >  NZ_CP009273/4171655‑4171808

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