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%	240.4 / NA	64	D185H (GAC→CAC)	rpoB	DNA‑directed RNA polymerase subunit beta
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base C (42/22);  total (42/22)

CTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCC  >  NZ_CP009273/4171646‑4171809                                                                                |                                                                                     cTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAAtt                                                                            >  2:1317151/1‑90 (MQ=255)    cGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGa                                                                         >  2:1011398/1‑90 (MQ=255)         aaGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGa                                                                    <  2:1225174/90‑1 (MQ=255)                 aCCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAc                                                            <  2:410306/90‑1 (MQ=255)                     aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  1:643999/1‑90 (MQ=255)                     aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  2:1128963/1‑90 (MQ=255)                     aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  1:353503/1‑90 (MQ=255)                     aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  2:343694/1‑90 (MQ=255)                     aCTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGt                                                        >  2:1283704/1‑90 (MQ=255)                         ttCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgt                                                    <  1:1528878/90‑1 (MQ=255)                                aaaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           <  1:1311614/76‑1 (MQ=255)                                aaaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           >  2:1311614/1‑76 (MQ=255)                                aaaGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATc                                             >  2:1237240/1‑90 (MQ=255)                                      cTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       <  2:567257/90‑1 (MQ=255)                                      cTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       >  2:769058/1‑90 (MQ=255)                                       tGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTc                                      >  2:869105/1‑90 (MQ=255)                                       tGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTc                                      >  1:990237/1‑90 (MQ=255)                                       tGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTc                                      >  1:1115837/1‑90 (MQ=255)                                        gTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCg                                     >  2:1390967/1‑90 (MQ=255)                                        gTATAACGCGCGTATCATCCCTTAACGTGGTTCCTGGCTGCACTGCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCg                                     >  2:64139/1‑90 (MQ=255)                                            aaCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGt                                 >  1:542928/1‑90 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           <  1:910052/63‑1 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           >  2:99914/1‑63 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           <  1:99914/63‑1 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           >  2:910052/1‑63 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgta                                                   >  2:1473367/1‑71 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTcgta                                                   <  1:1473367/71‑1 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  2:152321/1‑90 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  2:232374/1‑90 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  2:399970/1‑90 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  2:1063839/1‑90 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  1:1362082/1‑90 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTa                                >  1:309787/1‑90 (MQ=255)                                             aCGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCATACGTATCGACCGTCGCCGTa                                >  2:583292/1‑90 (MQ=255)                                              cgcgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTAGTACGTATCGACCGTCGCCGTaa                               >  2:1146565/1‑90 (MQ=255)                                               gcgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       >  2:623273/1‑81 (MQ=255)                                               gcgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGt                                       <  1:623273/81‑1 (MQ=255)                                               gcgcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTaaa                              >  1:1507955/1‑90 (MQ=255)                                                cgcgTATCATCCCTTACCGTGGTTCCTGGCTTCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAAc                             >  1:1462034/1‑90 (MQ=255)                                                 gcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAAct                            >  1:1352606/1‑90 (MQ=255)                                                 gcgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAAct                            <  2:275970/90‑1 (MQ=255)                                                  cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           <  1:41495/58‑1 (MQ=255)                                                  cgTATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAAcc                                                           >  2:41495/1‑58 (MQ=255)                                                    tATCATCCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAActgcc                         <  1:941017/90‑1 (MQ=255)                                                        atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCg                     <  1:341066/90‑1 (MQ=255)                                                        atcCCTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCg                     <  2:1141655/90‑1 (MQ=255)                                                          cccTTCCCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGAc                   <  1:610112/90‑1 (MQ=255)                                                          cccTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGAc                   <  1:387847/90‑1 (MQ=255)                                                            cTTACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACca                 >  1:1064797/1‑90 (MQ=255)                                                            cTTACCGTGGTTCCTAGCTGCACTTCGAATTCGAGCCGAAGGACAACCTGTTCATACGTATCGACCGTCGCCGTAAACTCCCTGCGAcga                 >  2:776104/1‑88 (MQ=255)                                                             ttACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACcat                >  2:1483237/1‑90 (MQ=255)                                                             ttACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACcat                >  1:1486484/1‑90 (MQ=255)                                                              tACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               >  1:440188/1‑42 (MQ=255)                                                              tACCGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGAc                                                               <  2:440188/42‑1 (MQ=255)                                                                 cGTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCAtt            >  1:420358/1‑90 (MQ=255)                                                                  gTGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTc           >  1:695212/1‑90 (MQ=255)                                                                   tGGTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCCCCGTAAACTGCCTGCGACCATCATTCt          >  2:1577286/1‑90 (MQ=255)                                                                    ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTg         >  2:128467/1‑90 (MQ=255)                                                                    ggTTCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTg         >  1:1035001/1‑90 (MQ=255)                                                                      ttCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcg       <  2:420358/90‑1 (MQ=255)                                                                      ttCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcg       <  2:1240449/90‑1 (MQ=255)                                                                       tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcgc      <  1:424424/90‑1 (MQ=255)                                                                       tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcgc      <  1:769058/90‑1 (MQ=255)                                                                       tCCTGGCTGCACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTgcgc      <  2:1035001/90‑1 (MQ=255)                                                                           ggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGccc  <  2:292971/90‑1 (MQ=255)                                                                           ggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGccc  <  1:669413/90‑1 (MQ=255)                                                                           ggctgcACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGccc  <  1:1237240/90‑1 (MQ=255)                                                                                |                                                                                     CTCCGACAAAGGTAAAACCCACTCTTCGGGTAAAGTGCTGTATAACGCGCGTATCATCCCTTACCGTGGTTCCTGGCTGGACTTCGAATTCGATCCGAAGGACAACCTGTTCGTACGTATCGACCGTCGCCGTAAACTGCCTGCGACCATCATTCTGCGCGCCC  >  NZ_CP009273/4171646‑4171809

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